Search NASA⌕ Search

SEARCH · Search NASA

Results for “Potentials”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Comparison of DeePMD, MTP, GAP, ACE and MACE Machine‐Learned Potentials for Radiation‐Damage Simulations: A User Perspective

Accurate and efficient interatomic potentials are essential for molecular dynamics (MD) simulations of radiation damage, gas diffusion, and phase stability in complex ceramics such as LiAlO 2 , especially under extreme conditions relevant to tritium production. Here, we evaluate the performance of six machine-learned interatomic potentials (MLIPs), moment tensor potential (MTP), Gaussian approximation potential, deep potential (DeePMD), atomic cluster expansion (ACE), message-passing ACE (multilayer atomic cluster expansion (MACE) pretrained) and MACE (trained from-scratch), all trained on the same density functional theory dataset with inclusion of tritium. The MLIPs are benchmarked against traditional Buckingham and ReaxFF potentials in terms of energy accuracy, density predictions, thermal equilibration behavior, threshold displacement energy (E d ), tritium diffusivity, and computational cost. Among the models, MTP shows the best overall balance between efficiency and accuracy, with low force and energy errors and realistic E d values for Li and Al. The ACE and MACE (pretrained and trained from scratch) models exhibit high E d (>200 eV) and unphysical pair interactions. DeePMD underestimates Ed due to overly repulsive behavior even at equilibrium distances. All models over-estimate tritium diffusion but the pretrained MACE model behaves well during tritium-diffusion simulations up to 500 K, maintaining diffusivities in the physically consistent 10 −11 m 2 /s range. Finally, we quantify the computational cost of each potential in large-scale atomic/molecular massively parallel simulator, finding that only MTP is more efficient than traditional empirical potentials, while others are significantly more expensive. These findings explain the trade-offs between accuracy and computational cost in MLIP development and provide essential guidance for use in high-throughput radiation damage and gas diffusion simulations in nuclear ceramics.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Review and construction of interatomic potentials for molecular dynamics studies of hydrogen embrittlement in Fe─C based steels

Abstract Reducing hydrogen embrittlement in the low‐cost Fe─C based steels have the potential to significantly impact the development of hydrogen energy technologies. Molecular dynamics studies of hydrogen interactions with Fe─C steels provide fundamental information about the behavior of hydrogen at microstructural length scales, although such studies have not been performed due to the lack of an Fe─C─H ternary interatomic potential. In this work, the literature on interatomic potentials related to the Fe─C─H systems are reviewed with the aim of constructing an Fe─C─H potential from the published binary potentials. We found that Fe─C, Fe─H, and C─H bond order potentials exist and can be combined to construct an Fe─C─H ternary potential. Therefore, we constructed two such Fe─C─H potentials and demonstrate that these ternary potentials can reasonably capture hydrogen effects on deformation characteristics and deformation mechanisms for a variety of microstructural variations of the Fe─C steels, including martensite that results from γ to α phase transformation, and pearlite that results from the eutectic formation of the Fe 3 C cementite compound.

Zhou, Xiaowang↗
DOE OSTI2025

Insight Into the High-Potential Branch of the Alternate Nad+-Dependent NADPH:Ferredoxin Oxidoreductase II (NfnII) from Pyrococcus Furiosus. Evidence for the Gating Step in Electron Bifurcation

We have investigated the rapid-reaction kinetics of the NAD+-dependent NADPH:ferredoxin oxidoreductase II (NfnII) from Pyrococcus furiosus, permitting a comparison with recent work done with the paralog NfnI from the same organism. The half-potentials of the electron-bifurcating L-FAD are highly crossed in both NfnI, meaning the potential for the quinone/semiquinone couple is significantly lower than that for the semiquinone/hydroquinone couple so that the semiquinone oxidation state is thermodynamically unstable. The same appears to be the case with NfnII on the basis of its similar behavior in transient absorption spectroscopy experiments and the absence of any evidence for FAD*- accumulation in the course of reductive titrations (which would be manifested as a transient increase in absorbance at ~380 nm)1. Reductive titrations with both the one-electron donor sodium dithionite and the obligate two-electron donor NADPH demonstrate that, in contrast to NfnI, little FADH* accumulates in the high-potential pathway of NfnII in the course of reduction, as reflected in the absence of a transient increase in absorbance in the 500-600 nm region. This indicates that the half-potentials of the S-FAD are crossed in NfnII by a minimum of 120 mV. Rapid-reaction experiments mixing oxidized NfnII with NADPH also show no evidence of S-FADH* accumulation. Furthermore, the enzyme is only partially reduced at the end of the reaction with NADPH, indicating that there is little electron transfer into the high-potential pathway of NfnII. When the reaction is carried out in the presence of the one-electron carrier ferredoxin, only approximately one equivalent of ferredoxin becomes reduced, in contrast to the maximum of three equivalents seen with NfnI. This is consistent with only a single electron bifurcation event taking place under these conditions with NfnII, with only a single high-potential electron passing into the [2Fe-2S] cluster of the high-potential pathway of NfnII but not progressing further to the S-FAD. This accounts for the observation that NfnII has greatly reduced bifurcating activity compared to NfnI. The lack of electron transfer into the S-FAD of NfnII due to its crossed half-potentials prevents successive bifurcation activity and indicates that electron transfer within the high-potential branch gates bifurcation activity in NfnI.

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

Shadow energy functionals and potentials in Born–Oppenheimer molecular dynamics

In Born–Oppenheimer molecular dynamics (BOMD) simulations based on the density functional theory (DFT), the potential energy and the interatomic forces are calculated from an electronic ground state density that is determined by an iterative self-consistent field optimization procedure, which, in practice, never is fully converged. The calculated energies and forces are, therefore, only approximate, which may lead to an unphysical energy drift and instabilities. Here, we discuss an alternative shadow BOMD approach that is based on backward error analysis. Instead of calculating approximate solutions for an underlying exact regular Born–Oppenheimer potential, we do the opposite. Instead, we calculate the exact electron density, energies, and forces, but for an underlying approximate shadow Born–Oppenheimer potential energy surface. In this way, the calculated forces are conservative with respect to the approximate shadow potential and generate accurate molecular trajectories with long-term energy stabilities. We show how such shadow Born–Oppenheimer potentials can be constructed at different levels of accuracy as a function of the integration time step, δt, from the constrained minimization of a sequence of systematically improvable, but approximate, shadow energy density functionals. For each energy functional, there is a corresponding ground state Born–Oppenheimer potential. These pairs of shadow energy functionals and potentials are higher-level generalizations of the original “zeroth-level” shadow energy functionals and potentials used in extended Lagrangian BOMD. The proposed shadow energy functionals and potentials are useful only within this extended dynamical framework, where also the electronic degrees of freedom are propagated as dynamical field variables together with the atomic positions and velocities. The theory is quite general and can be applied to MD simulations using approximate DFT, Hartree–Fock, or semi-empirical methods, as well as to coarse-grained flexible charge models.

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

Thermophysical properties of FLiBe using moment tensor potentials

Fluoride salts are prospective materials for applications in some next-generation nuclear reactors and their thermophysical properties at various conditions are of interest. Experimental measurement of the properties of these salts is often difficult and, in some cases, unfeasible due to challenges from high temperatures, impurity control, and corrosivity. Therefore, accurate theoretical methods are needed for fluoride salt property prediction. In this work, we used moment tensor potentials (MTP) to approximate the potential energy surface of eutectic FLiBe (66.6% LiF – 33.3% BeF2) predicted by the ab initio (DFT-D3) method. Here, we then used the developed potential and molecular dynamics to obtain several thermophysical properties of FLiBe, including radial distribution functions, density, self-diffusion coefficients, thermal expansion, specific heat capacity, bulk modulus, viscosity, and thermal conductivity. Our results show that the MTP potential approximates the potential energy surface accurately and the overall approach yields very good agreement with experimental values. The converged fitting can be obtained with less than 600 configurations generated from DFT calculations, which data can be generated in just 1200 core hours on today's typical processors. The MTP potential is faster than many machine learning potentials and about one order of magnitude slower than widely used empirical molten salt potentials such as Tosi/Fumi.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Development of an interatomic potential for the W–Ta system

A physics-inspired, data-driven interatomic potential framework for multi-element system is presented. This potential is based on the generalization of the embedded atom method potential and systematically incorporates two-, three- and many-body effects. Different atomic environments are described by atom-centered Gaussian basis sets which provides sufficient flexibility to capture diverse atomic environments and allows for easy optimization of the free parameters. An interatomic potential model for the tungsten–tantalum (W–Ta) system is developed using this framework by training on data from ab initio density functional theory (DFT) calculations. The potential is thoroughly tested at various compositions by comparing bulk and defect properties from experimental and DFT data. It is shown that the potential predicts the elastic constants, defect properties, such as vacancy and interstitial formation energies, core structures of dislocations, and melting points for both pure and tantalum and tungsten with an accuracy comparable to other single element machine learning based potentials. In conclusion, the potential model is used to investigate the formation energies of ordered alloys and vacancy formation energies for chemically random W–Ta alloys.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Nanometer-scale electrical potential imaging on absorber of CdSeTe solar cells

Here, we report on nm-scale electrical potential imaging throughout As-doped and Cu-doped CdSeTe absorbers using Kelvin probe force microscopy (KPFM). The potential imaging was conducted both laterally and vertically on beveled films using ion milling at small glancing angle. KPFM images electrical potential on the beveled surface and assesses defect charging in the subsurface region within a screening length from the beveled surface. We found that the grain boundaries were positively charged and that there were significant potential fluctuations in both grain boundary versus grain interior and intragrain. We further found that these potential fluctuations decreased significantly toward the front interface. Time of flight secondary ion mass spectrometry imaging shows that Se content increased toward the front interface, consistent with Se passivation of defects. The potential fluctuation was induced by defect charging, and the results elucidate different details of the defect configurations and grain structures of the films with different CdCl 2 treatment temperatures in the As-doped CdSeTe. The defect configurations in the region near the front interface can be a main factor contributing to the device performance difference. Our potential imaging provides insights about the defects throughout the absorber films, and shows that the potential fluctuation has a direct correlation to the V oc deficit.

14 SOLAR ENERGY↗
DOE OSTI2021

Continental-Scale Geographic Trends in Barometric-Pumping Efficiency Potential: A North American Case Study

Barometric pumping is a gas transport mechanism that has important implications for many applications involving subsurface gas seepage processes. Here, this study provides the first continental-scale analysis of barometric-pumping efficiency potential based on meteorology. We quantified the barometric-pumping efficiency potential at 1,257 locations across the continental US and Canada. The results provide continental-scale geographic dependencies of barometric-pumping efficiency potential, indicating a significant correlation with latitude and a nonlinear dependence on longitude. The analysis also indicates that variability in barometric-pumping efficiency potential decreases with distance from the coast and as elevation increases. Locations far from the coastline are more likely to have upper mid-range potentials, while higher elevation locations are more likely to have low potentials. The highest barometric-pumping efficiency potentials are mostly found around the Gulf of St. Lawrence around 50°N. Locations along the Atlantic coast exhibit large-scale variations in potentials with a clear increasing trend with latitude.

58 GEOSCIENCES↗
DOE OSTI2023

Hybrid interatomic potential for Sn

To design materials for extreme applications, it is important to understand and predict phase transitions and their influence on material properties under high pressures and temperatures. Atomistic modeling can be a useful tool to assess these behaviors. However, this can be difficult due to the lack of fidelity of the interatomic potentials in reproducing this high pressure and temperature extreme behavior. Here, in this work, a hybrid EAM-R—which is the combination of embedded atom method (EAM) and rapid artificial neural network potential—for Tin (Sn) is described which is capable of accurately modeling the complex sequence of phase transitions between different metallic polymorphs as a function of pressure. This hybrid approach ensures that a basic empirical potential like EAM is used as a lower energy bound. By using the final activation function, the neural network contribution to energy must be positive, assuring stability over the whole configuration space. This implementation has the capacity to reproduce density functional theory results at 6 orders of magnitude slower than a pair potential for molecular dynamics simulation, including elastic and plastic characteristics and relative energies of each phase. Using calculations of the Gibbs free energy, it is demonstrated that the potential precisely predicts the experimentally observed phase changes at temperatures and pressures across the whole phase diagram. At 10.2 GPa, the present potential predicts a first-order phase transition between body-centered tetragonal (BCT) β-Sn and another polymorph of BCT-Sn. This structure transforms into body-centered cubic near the experimentally reported value at 33 GPa. Thus, the Sn potential developed in this paper can be used to study complex deformation mechanisms under extreme conditions of high pressure and strain rates unlike existing potentials. Moreover, the framework developed in this paper can be extended for different material systems with complex phase diagrams.

36 MATERIALS SCIENCE↗
DOE OSTI2020

In Operando Calorimetric Measurements for Activated Carbon Electrodes in Ionic Liquid Electrolytes under Large Potential Windows

This study aims to investigate the effect of the potential window on heat generation in carbon-based electrical double layer capacitors (EDLCs) with ionic-liquid (IL)-based electrolytes using in operando calorimetry. The EDLCs consisted of two identical activated-carbon electrodes with either neat 1-butyl-1-methylpyrrolidinium bis(trifluoromethane-sulfonyl)imide ([Pyr 14 ][TFSI]) electrolyte or 1.0 m [Pyr 14 ][TFSI] in propylene carbonate (PC) as electrolyte. The instantaneous heat generation rate at each electrode was measured under galvanostatic cycling for different potential windows ranging from 1 to 4 V. First, the heat generation rates at the positive and negative electrodes differed significantly in neat IL owing to the differences in the ion sizes and diffusion coefficients. However, these differences were minimized when the IL was diluted in PC. Second, for EDLC in neat [Pyr 14 ][TFSI] at high potential window (4 V), a pronounced endothermic peak was observed at the beginning of the charging step at the positive electrode owing to TFSI - intercalation in the activated carbon. On the other hand, for EDLC in 1.0 m [Pyr 14 ][TFSI] in PC at potential window above 3 V, an endothermic peak was observed only at the negative electrode owing to the decomposition of PC. Third, for both neat and diluted [Pyr 14 ][TFSI] electrolytes, the irreversible heat generation rate increased with increasing potential window and exceeded Joule heating. This was attributed to the effect of potential-dependent charge redistribution resistance. Additionally, a further increase in the irreversible heat generation rate was observed for the largest potential windows owing to the degradation of the PC solvent. Finally, for both types of electrolyte, the reversible heat generation rate increased with increasing potential window because of the increase in the amount of ion adsorbed/desorbed at the electrode/electrolyte interface.

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

Long-term carbon intensity reduction potential of K-12 school buildings in the United States

School buildings have a great potential for carbon emission reduction since their annual emission is about 72 million metric tons. Currently, more than 30% of school buildings were built before 1960 and are underperforming. To effectively reduce carbon emissions via school building retrofits, it is critical for policymakers to understand the carbon intensity reduction potential of retrofitting school buildings in different regions. Hence, this study develops a method to comprehensively assess the long-term carbon intensity reduction potential of aggregated commercial buildings on a county-by-county basis in the continental U.S. We apply this method to the K-12 school buildings including primary and secondary school buildings. Here, this paper predicts the carbon intensity reduction potential of K-12 school buildings with eight building retrofit measures from 2022 to 2050 in the continental U.S. The results reveal several interesting findings: (1) In the approximately 3,000 counties of the U.S. from 2022 to 2050, the carbon intensity reduction potential of retrofitting K-12 school buildings in each county ranges from 0.41 kg/m 2 to 40.00 kg/m 2 . (2) Even in the same climate zone, the trends of carbon intensity reduction potential from 2022 to 2050 are different depending on their electricity sources. For example, in a hot and humid climate zone, the carbon intensity reduction potential in Florida will decrease from 2044 to 2048. However, in Mississippi, the carbon intensity reduction potential from 2044 to 2046 will increase due to the termination of the nuclear energy usage. (3) Generally, reducing lighting power density leads to more carbon intensity reduction in most states, but it might not be applicable for states with high clean energy penetration, such as Washington.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2024

Prediction of Redox Potentials for Ac, Th, and Pa in Aqueous Solution

Density functional theory in conjunction with small core pseudopotentials and the associated basis sets was used to calculate potentials for multiple redox couples, covering a range of oxidation states for Ac (0 to III), Th (0 to IV), and Pa (0 to V) in aqueous solution. Solvation effects were incorporated using a supermolecule-continuum approach, with 30 water molecules representing two solvation shells, and the COSMO and SMD implicit solvation models. The calculated geometries for Ac(III), Th(IV), and Pa(V) were in reasonable agreement with the available experimental data. Using the COSMO model with the B3LYP functional, the calculated redox potentials were within ± 0.2 V from experiment for most redox couples. Several pathways were explored for the Pa(V/IV) redox couple for different forms of Pa(V) and Pa(IV). Most Pa(V/IV) redox couples have very similar potentials, ranging from 0 to -0.4 V up to a pH of 1.4. At pH = 1.4, the potentials shift to values that are more negative than -0.7 V, reflecting the growing unfavorable nature of the redox process at higher pH levels. The calculated values for An(III/II) potentials were consistent with prior estimates and the available experimental data. The predicted redox potentials for An(II/I) were highly negative, as expected. For An(I/0) potentials, Th and Pa exhibited positive values, contrasting with the negative values calculated for Ac. Furthermore, the An +m /An(0) potentials agreed better with the experimental data when using the COSMO solvation model as compared to the SMD model.

Chemical calculations↗
DOE OSTI2022

Plasma parameters and the reduction potential at a plasma-liquid interface

We report nonthermal plasmas in contact with liquids have been shown to generate a variety of reactive species capable of initiating reduction–oxidation (redox) reactions at the electrochemically active plasma–liquid interface. In conventional electrochemical cells, selective redox chemistry is achieved by controlling the reduction potential at the solid electrode–electrolyte interface by applying a bias via an external circuit. In the case of plasma–liquid systems, an analogous means of tuning the reduction potential near the interface has not clearly been identified. When treated as a floating surface, the liquid is expected to adopt a net negative charge to balance the flux of hot electrons and relatively cold positive ions. The reduction potential near the plasma–liquid interface is hypothesized to be proportional to the floating potential, which can be approximated using an analytical model provided the plasma parameters are known. Herein, we present a framework for correlating the electron density and electron temperature of a noble gas plasma jet to the reduction potential near the plasma–liquid interface. The plasma parameters were acquired for an argon atmospheric plasma jet in contact with an aqueous solution by means of laser Thomson scattering. The reduction potential was determined using identical reference electrodes to measure the potential difference between the plasma–liquid interface and bulk solution. Interestingly, the measured reduction potentials near the plasma–liquid interface were found to be in good agreement with the model-predicted values determined using the plasma parameters obtained from the Thomson scattering experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Limits to forests-based mitigation in integrated assessment modelling: global potentials and impacts under constraining factors

Forests-based measures such as afforestation/reforestation (A/R) and reducing deforestation (RDF) are considered promising options to mitigate climate change, yet their mitigation potentials are limited by economic and biophysical factors that are largely uncertain. The range of mitigation potential estimates from integrated assessment models raises concerns about the capacity of land systems to provide realistic, cost-effective and permanent land-based mitigation. We use the Global Change Analysis Model to quantify the economic mitigation potential of forests-based measures by simulating a climate policy including a tax on greenhouse gas emissions from agriculture, forestry, and other land uses. In addition, we assess how constraining unused arable land (UAL) availability, forestland expansion rates, and global bioenergy demand may influence the forests-based mitigation potential by simulating scenarios with alternative combinations of constraints. Results show that the average forests-based mitigation potential in 2020–2050 increases from 738 MtCO 2 .yr -1 through a forestland increase of 86 Mha in the fully constrained scenario to 1394 MtCO 2 .yr -1 through a forestland increase of 146 Mha when all constraints are relaxed. Regional potentials in terms of A/R and RDF differ strongly between scenarios: unconstrained forest expansion rates mostly increase A/R potentials in northern regions (e.g., +120 MtCO 2 .yr -1 in North America); while unconstrained UAL conversion and low bioenergy demand mostly increase RDF potentials in tropical regions (e.g., +76 and +68 MtCO 2 .yr -1 in Southeast Asia, respectively). This study shows that forests-based mitigation is limited by many factors that constrain the rates of land use change across regions. These factors, often overlooked in modelling exercises, should be carefully addressed for understanding the role of forests in global climate mitigation and defining pledges towards the Paris Agreement.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Deep potential molecular dynamics simulations of low-temperature plasma-surface interactions

Machine learning approaches to potential generation for molecular dynamics (MD) simulations of low-temperature plasma-surface interactions could greatly extend the range of chemical systems that can be modeled. Empirical potentials are difficult to generalize to complex combinations of multiple elements with interactions that might include covalent, ionic, and metallic bonds. This work demonstrates that a specific machine learning approach, Deep Potential Molecular Dynamics (DeepMD), can generate potentials that provide a good model of plasma etching in the Si-Cl-Ar system. Comparisons are made between MD results using DeepMD models and empirical potentials, as well as experimental measurements. Pure Si properties predicted by the DeepMD model are in reasonable agreement with experimental results. Simulations of Si bombardment by Ar + ions demonstrate the ability of the DeepMD method to predict sputtering yields as well as the depth of the amorphous-crystalline interface. Etch yields as a function of flux ratio and ion energy for simultaneous Cl 2 and Ar + impacts are in good agreement with previous simulation results and experiment. Predictions of etch yields and etch products during plasma-assisted atomic layer etching of Si-Cl 2 -Ar are shown to be in good agreement with MD predictions using empirical potentials and with experiment. Finally, good agreement was also seen with measurements for the spontaneous etching of Si by Cl atoms at 300 K. Further, the demonstration that DeepMD can reproduce results from MD simulations using empirical potentials is a necessary condition to future efforts to extend the method to a much wider range of systems for which empirical potentials may be difficult or impossible to obtain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL Radiological Response Training Range

This assessment uses screening-level models to calculate potential environmental impacts from proposed tests at the Idaho National Laboratory (INL) Radiological Response Training Range (RRTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K2O, LaBr3, KBr, Cu, Zr, F, Ga, Ga2O3, NaNO2, Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface-soil impacts were calculated with the computer program mixing-cell model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated assuming 12 tests are conducted annually using all 11 material types for a period of 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer-code input, output, executable files, batch files, scripts, and spreadsheet files are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS). It is highly unlikely the test scenarios evaluated in this ECAR will adversely impact human health based on comparisons of calculated dose and concentration against regulatory standards and risk-based screening levels. Conservative estimates of dose to workers and the public from atmospheric transport of possible radionuclide releases are far below federal radiation protection standards. Conservative estimates of potential contaminant concentrations in groundwater are less than federal drinking water standards or screening levels. Predicted radionuclide concentrations in surface soils are below risk-based screening levels, except for Ge-68 (material Ga-68) for the worker. The Ge-68 soil concentration can be made less than the worker PRG, if the number of annual tests using Ga-68 is reduced from 12 to 6. However, the sum of ratios still exceeds one because of the high K-40 ratio. If the EF of the worker (number of days the worker is in the contaminated testing area) is reduced from 225 days/yr (default value for full time worker) to 112 days/yr, the Ge-68 ratio is less than one and the sum of ratios is less than one. Actual radiation doses and groundwater and surface-soil concentrations are likely to be much less than those calculated because of the conservative assumptions and parameters employed in the modeling. For example, atmospheric-transport calculations assume the entire inventory of each material type is readily released to the atmosphere and no plume deposition, depletion, or radioactive decay occurs during transport. The calculations also assume the same meteorological conditions (e.g., wind velocity, wind direction, stability class) that produce the maximum 95th percentile concentration (i.e., concentration representing the 95th percentile of a distribution of concentrations derived from 3 years of hourly meteorological data) at each receptor location are the same for all 12 tests during the year, and each receptor is assumed to be present during all 12 tests. The surface-soil assessment assumes the entire inventory of each test is deposited in the top 5 cm of soil. No atmospheric dispersal is assumed, and the radionuclides are subject only to leaching and radioactive decay. The groundwater-pathway modeling is conservative in that it is one-dimensional in the unsaturated zone (no lateral spreading/dilution) and assumes the entire inventory of contaminants infiltrates into the ground at the same location for every test. This is especially conservative for particulate radionuclides because they would have to dissolve or corrode first and some would be dispersed into the atmosphere. The groundwater receptor is also assumed to consume water directly from a hypothetical well positioned in the location of maximum concentration. In addition, conservative degradation rates were used, and volatilization was not considered for the nonradioactive chemical

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL National Security Test Range

This assessment uses screening-level models to calculate potential environmental impacts from proposed tests at two locations at the Idaho National Laboratory (INL) National Security Test Range (NSTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K 2 O, LaBr 3 , KBr, Cu, Zr, F, Ga, Ga 2 O 3 , NaNO 2 , Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface-soil impacts were calculated with the computer program mixing-cell model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated based on the assumption that 12 tests are conducted annually using all 11 material types for 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer-code input, output, executable files, batch files, scripts, and spreadsheet files, are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS). It is highly unlikely the test scenarios evaluated in this ECAR will adversely impact human health based on comparisons of calculated dose and concentration against regulatory standards and risk-based screening levels. Conservative estimates of dose to workers and the public from atmospheric transport of possible radionuclide releases are far below federal radiation protection standards. Conservative estimates of potential contaminant concentrations in groundwater are less than federal drinking water standards or screening levels. Predicted radionuclide concentrations in surface soils are below risk-based screening levels, except for Ge-68 (material Ga-68) for the worker. The Ge-68 soil concentration can be made less than the worker PRG, if the number of annual tests using Ga-68 is reduced from 12 to 6. However, the sum of ratios still exceeds one because of the high K-40 ratio. If the EF of the worker (number of days the worker is in the contaminated testing area) is reduced from 225 days/yr (default value for full time worker) to 112 days/yr, the Ge-68 ratio is less than one and the sum of ratios is less than one. Actual radiation doses and groundwater and surface-soil concentrations are likely to be much less than those calculated because of the conservative assumptions and parameters employed in the modeling. For example, atmospheric-transport calculations assume the entire inventory of each material type is readily released to the atmosphere and no plume deposition, depletion, or radioactive decay occurs during transport. The calculations also assume the same meteorological conditions (e.g., wind velocity, wind direction, stability class) that produce the maximum 95th percentile concentration (i.e., concentration representing the 95th percentile of a distribution of concentrations derived from 3 years of hourly meteorological data) at each receptor location are the same for all 12 tests during the year, and each receptor is assumed to be present during all 12 tests. The surface-soil assessment assumes the entire inventory of each test is deposited in the top 5 cm of soil. No atmospheric dispersal is assumed, and the radionuclides are subject only to leaching and radioactive decay. The groundwater-pathway modeling is conservative in that it is one-dimensional in the unsaturated zone (no lateral spreading/dilution) and assumes the entire inventory of contaminants infiltrates into the ground at the same location for every test. This is especially conservative for particulate radionuclides because they would have to dissolve or corrode first and some would be dispersed into the atmosphere. The groundwater receptor is also assumed to consume water directly from a hypothetical well positioned in the location of maximum concentration. In addition, conservative degradation rates were used, and volatilization was not considered for the nonradradioactive chemicals modeled. And finally, the calculations assume all 12 tests will be performed at the same place at both locations, and all 11 radioactive material types will be used for each test. This is conservative because it is anticipated that no more than two material types will be used per test.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2022

Quantifying the drivers of ecosystem fluxes and water potential across the soil-plant-atmosphere continuum in an arid woodland

Dryland ecosystems occupy a vast swath of the terrestrial land surface and exert a sizeable impact on the cycling of carbon and water globally. These biomes are characterized by tightly coupled carbon and water cycles that respond rapidly to transient pulses in water availability. However, there exist many mechanistic uncertainties regarding the environmental drivers of, and linkages between, plant and ecosystem processes. Thus, drylands are often poorly represented in many vegetation and land surface models. An enhanced understanding of dryland ecosystem function is limited by the lack of long-term, co-located, and frequent measurements of plant and ecosystem processes. At a piñon-juniper woodland in southeastern Utah, USA, we collected a continuous dataset of meteorological conditions, soil water potential from surface to bedrock, tree water potential, and ecosystem carbon and water fluxes from eddy covariance. We found that predawn and midday tree water potential and daily ecosystem fluxes were highly sensitive to fluctuations in soil water availability, particularly in shallower layers, and that daytime variability in atmospheric drivers only loosely controlled these processes. The strong connections between shallow soil water potential, tree water potential, and ecosystem fluxes occurred because of the dominant role of precipitation pulses in driving vegetation activity, as even small pulses of moisture stimulated shallow soil water potential, tree water potential, and evapotranspiration for between 1 and 2 weeks. Carbon fluxes (net ecosystem exchange and gross primary productivity) were sensitive to precipitation pulses for longer, up to 3 weeks. Our results highlight that improved monitoring and sensing of shallow soil moisture can greatly enhance our understanding of dryland ecosystem function. Here, a better mechanistic understanding of the impacts of precipitation pulses is also needed to improve vegetation modeling of dryland ecosystems.

54 ENVIRONMENTAL SCIENCES↗