Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “gradient 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 1,369 records · Page 76
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

An error-in-constitutive equations strategy for topology optimization for frequency-domain dynamics

This paper presents a topology optimization formulation for frequency-domain dynamics to reduce solution dependence upon initial guess and considered loading conditions. Due to resonance phenomena in undamped steady-state dynamics, objectives measuring dynamic response possess many local minima that may represent poor solutions to a design problem, an issue exacerbated for design with respect to multiple frequencies. In this work, we propose an extension of the modified error-in-constitutive-equations (MECE) method, used previously in material identification inverse problems, as a new approach for frequency-domain dynamics topology optimization to mitigate these issues. The main idea of the proposed framework is to incorporate an additional penalty-like term in the objective function that measures the discrepancy in the constitutive relations between stresses and strains and between inertial forces and displacements. Then, the design problem is cast within a PDE-constrained optimization formulation in which we seek displacements, stresses, inertial forces, and a density-field solution that minimize our new objective subject to conservation of linear momentum plus some additional constraints. We show that this approach yields superior designs to conventional gradient-based optimization approaches that solely use a functional of displacements as the objective, while strictly enforcing the constitutive equations. The MECE strategy integrates into a density-based topology optimization scheme for void–solid or two-phase material structural design. We highlight the merits of our approach in a variety of scenarios for direct frequency response design, considering multiple frequency load cases and structural objectives.

42 ENGINEERING↗
DOE OSTI2026

Avoidance of disruptions on KSTAR due to vertical displacement events via novel real-time stability assessment

Disruption avoidance via the DECAF approach has been achieved on KSTAR using a novel real-time vertical stability assessment and a multiactuator feedback control strategy. The development of disruption avoidance strategies with reactor-relevant reliability is an urgent activity, enabling future fusion power plants. The stability metric employed is based on a new formulation of a vertical force gradient balance metric evaluated across the poloidal cross section of the plasma, with parameters tuned using historical data. Evaluation of this metric on a validation set of 400 recent KSTAR shots indicates >82% of Vertical displacement events can be avoided via feedback control. Essential to its calculation is the two-dimensional toroidal current density distribution in the plasma. Measurement of this profile faster than fully-converged equilibrium reconstructions can deliver is found to improve forecaster performance and is achieved with a surrogate model that takes as input magnetic diagnostic measurements and outputs the current profile on a basis comprising the top principal components of historical current profiles (from past equilibrium reconstructions). This method solves the non-uniqueness problem typically faced when reconstructing current profiles directly from diagnostics, while improving computational time and accuracy. On average, profiles produced by this model reach coefficients of determination of >0.99 with respect to those from equilibrium reconstructions. The avoidance actuators employed include poloidal field coils and an electron cyclotron current drive system. The multiactuator approach, as shown in this first demonstration, allows disruption avoidance while minimizing impact to operational performance. This ability, along with its flexibility and speed, makes this new approach an attractive option for avoiding these types of disruptions in reactors.

Tobin, Matthew [Columbia Univ., New York, NY (Unit↗
DOE OSTI2023

DSPS: Differentiable stellar population synthesis

ABSTRACT Models of stellar population synthesis (SPS) are the fundamental tool that relates the physical properties of a galaxy to its spectral energy distribution (SED). In this paper, we present DSPS: a python package for SPS. All of the functionality in DSPS is implemented natively in the JAX library for automatic differentiation, and so our predictions for galaxy photometry are fully differentiable, and directly inherit the performance benefits of JAX, including portability onto GPUs. DSPS also implements several novel features, such as i) a flexible empirical model for stellar metallicity that incorporates correlations with stellar age, ii) support for the Diffstar model that provides a physically-motivated connection between the star formation history of a galaxy (SFH) and the mass assembly of its underlying dark matter halo. We detail a set of theoretical techniques for using autodiff to calculate gradients of predictions for galaxy SEDs with respect to SPS parameters that control a range of physical effects, including SFH, stellar metallicity, nebular emission, and dust attenuation. When forward modelling the colours of a synthetic galaxy population, we find that DSPS can provide a factor of 5 speed-up over standard SPS codes on a CPU, and a factor of 300-400 on a modern GPU. When coupled with gradient-based techniques for optimization and inference, DSPS makes it practical to conduct expansive likelihood analyses of simulation-based models of the galaxy–halo connection that fully forward model galaxy spectra and photometry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

High fidelity model of directed energy deposition: Laser-powder-melt pool interaction and effect of laser beam profile on solidification microstructure

Metal additive manufacturing technologies keep receiving a great deal of interest as well as strong requests to develop methods to link the process science to printed parts performance and understand how to overcome inherent limitations. Here, a high-fidelity model based on the multiphysics ALE3D code was developed to reproduce the directed energy deposition process down to the powder scale. This includes resolving the laser-powder-melt pool interactions (powder impingement and incorporation into melt pool, hydrodynamics flow condition and laser absorption inefficiencies) as well as the resulting solidification microstructure. This micrometer scale digital twin captured the effect of powder incorporation process and powder flow rate on porosity. Furthermore, it was used to explore how a ring laser beam profile instead of the standard Gaussian laser profile could decrease the thermal gradient along the solidification front in the melt pool, which in turn can increase propensity for more desirable equiaxed grains.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Cosmogenic 10Be constraints on deglacial snowline rise in the Southern Alps, New Zealand

We report geochronological dating of glacial landforms, such as terminal and lateral moraines, is useful for determining the extent and timing of past glaciation and for reconstructing the magnitude and rate of past climate changes. In the Southern Alps of New Zealand, well-dated glacial geomorphological records constrain the last glacial cycle across much of the Waitaki River basin (e.g. Ōhau, Pukaki, Tekapo) but its southern sector such as the Ahuriri River valley remains comparatively unconstrained. Recently, there has been debate on the scale and rapidity of mountain glacier retreat during the last glacial termination, particularly the 20–17 ka period in New Zealand. Missing from this debate is well-constrained equilibrium-line altitude (ELA) and associated temperature reconstructions, particularly over the period around 17 ka, which can help us to develop a more complete picture of how past temperature changes drove glacier retreat. Here we report the first glacial chronology dataset from the Last Glacial Maximum (LGM) and subsequent deglaciation from the Ahuriri River valley, Southern Alps, New Zealand (44°23'54"S, 169°39'48"E) based on 38 beryllium-10 ( 10 Be) surface-exposure ages from terminal moraine systems and glaciated bedrock situated at the lower and middle sections of the valley. Our results show that the former Ahuriri Glacier reached its maximum extent at 19.8 ± 0.3 ka, which coincides with the global Last Glacial Maximum. By 16.7 ± 0.3 ka, the glacier had retreat ~18 km up-valley suggesting at least ~43% glacier-length loss relative to its full LGM extent. This deglaciation was accompanied by the formation of a shallow proglacial lake. Using the accumulation area ratio (AAR) method, we estimate that the ELA was lower than present by ~880 m (~1120 m a.s.l.) at 19.8 ± 0.3 ka, and ~770 m lower (~1230 m a.s.l.) at 16.7 ± 0.3 ka. Applying an estimate for temperature lapse rate, this ELA anomaly implies that local air temperature was 5 ± 1 °C colder than present (1981–2010) at 19.8 ± 0.3 ka, while it was 4.4 ± 0.9 °C colder at 16.7 ± 0.3 ka, assuming no change in precipitation. The substantial glacier retreat in response to a relatively small accompanying increases in ELA (110 m) and temperature (0.6 °C) may have been a result of the high glacier-length sensitivity of this glacier system due to its low gradient of former ice surface. Our low warming estimate differs markedly from other deglaciation studies, specifically from Rakaia River valley, which reports a much larger temperature increase at the onset of the last deglaciation. This precisely-dated moraine record along with reconstructed ELA as proxies for atmospheric conditions, provides new insight into post LGM glacier behaviour and climate conditions in New Zealand.

58 GEOSCIENCES↗
DOE OSTI2024

Photon classification with Gradient Boosted Trees at CLAS12

Dihadron semi-inclusive deep inelastic scattering (SIDIS) of 10.6 GeV longitudinally polarized electrons off the proton has been measured using the CLAS12 detector at Jefferson Lab. Two separate channels, π + π 0 and π - π 0 , were analyzed, requiring the reconstruction of diphoton pairs. Here, in this analysis, we addressed the problem of false neutral particles being reconstructed by CLAS12's event builder, polluting the otherwise physical combinatorial background underneath the π 0 peak. A photon classifier using a Gradient Boosted Trees (GBTs) architecture was trained with Monte Carlo simulations to reduce the amount of background π 0 's. We show that the nearest-neighbor features learned by the model lead to a substantial increase in signal vs. background discrimination compared to previous CLAS12 π^0 analyses. The machine learning approach recovers several times more dihadron statistics for the dataset.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

A machine learning approach for efficient multi-dimensional integration

Many physics problems involve integration in multi-dimensional space whose analytic solution is not available. The integrals can be evaluated using numerical integration methods, but it requires a large computational cost in some cases, so an efficient algorithm plays an important role in solving the physics problems. We propose a novel numerical multi-dimensional integration algorithm using machine learning (ML). After training a ML regression model to mimic a target integrand, the regression model is used to evaluate an approximation of the integral. Then, the difference between the approximation and the true answer is calculated to correct the bias in the approximation of the integral induced by ML prediction errors. Because of the bias correction, the final estimate of the integral is unbiased and has a statistically correct error estimation. Three ML models of multi-layer perceptron, gradient boosting decision tree, and Gaussian process regression algorithms are investigated. The performance of the proposed algorithm is demonstrated on six different families of integrands that typically appear in physics problems at various dimensions and integrand difficulties. The results show that, for the same total number of integrand evaluations, the new algorithm provides integral estimates with more than an order of magnitude smaller uncertainties than those of the VEGAS algorithm in most of the test cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Design of a Thomson scattering diagnostic for the SMall Aspect Ratio Tokamak (SMART)

Here, we describe the design of a Thomson scattering (TS) diagnostic to be used on the SMall Aspect Ratio Tokamak (SMART). SMART is a spherical tokamak being commissioned in Spain that aims to explore positive triangularity and negative triangularity plasma scenarios at a low aspect ratio. The SMART TS diagnostic is designed to operate at high spatial resolution, 6 mm scattering length in the low-field side and 9 mm in the high-field side regions, and a wide dynamic range, electron temperature from 1 eV to 1 keV and density from 5 x 10 18 m -3 to 1 x 10 20 m -3 ⁠, to resolve large gradients formed at the plasma edge and in the scrape-off layer (SOL) under different triangularities and low aspect ratios. A 2 J @1064 nm laser will be used that is capable of operating in the burst mode at 1, 2, and 4 kHz to investigate fast phenomena and at 30 Hz to study 1 s (or more) long discharges. The scattered light will be collected over an angular range of 60° – 120° from 28 spatial points in the midplane covering the entire plasma width and the outer midplane SOL. Each scattering signal will be spectrally resolved on five wavelength channels of a polychromator to obtain the electron temperature measurement. We will also present a method to monitor in situ laser alignment in the core during calibrations and plasma operations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

AIVT: Inference of turbulent thermal convection from measured 3D velocity data by physics-informed Kolmogorov-Arnold networks

We propose the artificial intelligence velocimetry-thermometry (AIVT) method to reconstruct a continuous and differentiable representation of the temperature and velocity in turbulent convection from measured three-dimensional (3D) velocity data. AIVT is based on physics-informed Kolmogorov-Arnold networks and trained by optimizing a loss function that minimizes residuals of the velocity data, boundary conditions, and governing equations. We apply AIVT to a set of simultaneously measured 3D temperature and velocity data of Rayleigh-Bénard convection, obtained by combining particle image thermometry and Lagrangian particle tracking. This enables us to directly compare machine learning results to true volumetric, simultaneous temperature and velocity measurements. We demonstrate that AIVT can reconstruct and infer continuous, instantaneous velocity and temperature fields and their gradients from sparse experimental data at a high resolution, providing an additional approach for understanding thermal turbulence.

Science & Technology - Other Topics↗
DOE OSTI2025

Arctic shrub and Eriophorum leaf and root decomposition, northern Alaska, 2017-2018

This data package contains litter decomposition data collected from 170 plots of rapidly expanding shrub genera (Alnus, Betula, and Salix) and a widespread sedge (Eriophorum vaginatum) along a latitudinal and temperature gradient in northern Alaska. These data were produced from a litter bag experiment that took place from July 2017 to July 2018 and include mass loss and nitrogen loss decomposition metrics for both leaf and root litters. These raw data support a submitted manuscript that examines the variability in decomposition between shrub and graminoid leaf and root litters across a 1-year experiment across the graminoid-dominated Arctic tundra and reveals how deciduous shrub expansion affects litter decomposition in tundra ecosystems. Data are presented by site (n=5) and patch (shrub or sedge plot) in csv files. The site and plot location data and environmental measurement data are provided in Fraterrigo and Chen (2020). Additional methods regarding plot distribution and environmental measurements are in Chen et al. (2020) and Fraterrigo et al. (2024).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Ice nucleating particle concentrations measured by the PNNL-CFDC during the AGINSGP field experiment

Real-time, immersion mode ambient ice nucleating particle (INP) concentrations were measured at ARM’s Southern Great Plains (SGP) site during the AGINSGP field experiment using a continuous flow diffusion chamber (CFDC) from Pacific Northwest National Laboratory (PNNL). The PNNL-CFDC is an ice-thermal gradient diffusion chamber that optically detects the freezing of single aerosol particles (Kulkarni et al. 2020). The instrument was located in the guest instrument facility (GIF) and measured INP concentrations from April 8th to April 29th. Air was drawn into the GIF through custom aluminum stack inlets (6” inner diameter) attached vertically to the outer GIF wall and the GIF outside platform. A blower pulled air through the stacks at a velocity of 1 m/s. Aerosol instruments subsampled via wall ports, through either 3/8” or 1/4" copper lines, depending on sample flow of the instrument. Kulkarni, G., Hiranuma, N., Möhler, O., Höhler, K., China, S., Cziczo, D. J., & DeMott, P. J. (2020). A new method for operating a continuous-flow diffusion chamber to investigate immersion freezing: Assessment and performance study. Atmospheric Measurement Techniques, 13(12), 6631–6643. https://doi.org/10.5194/amt-13-6631-2020.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Unlocking bulk and surface oxygen transport properties of mixed oxide-ion and electron conducting membranes with combined oxygen permeation cell and oxygen probe method

Surface exchange kinetics and bulk diffusion of oxygen are of paramount importance to the activity of oxygen electrocatalysis and performance of electrochemical devices such as fuel cell, metal-air batteries, and oxygen separation membranes. Conventional approaches to obtaining these transport properties are often limited to single property under a specific non-operation related condition. Here we use a combined oxygen permeation cell and oxygen probe methodology to simultaneously attain rates of oxygen surface exchange and bulk conductivity/chemical diffusivity of three representative mixed oxide-ion and electron conductors, namely SrCo 0.9 Ta 0.1 O 3-δ (SCT), La 0.6 Sr 0.4 CoO 3-δ (LSC) and La 0.6 Sr 0.4 FeO 3-δ (LSF), operated under a steady-state oxygen flux. The results explicitly show that SCT exhibit the highest oxide-ion conductivity/chemical diffusivity, fastest rates of surface oxygen exchange kinetics, thus promising to be the best oxygen electrocatalyst. Here, we have also mapped out the distribution of oxygen chemical potential gradient across the membranes and applied B-transport number concept to illustrate the rate-limiting steps in the overall oxygen permeation process.

42 ENGINEERING↗
DOE OSTI2024

Geometric origin of the intrinsic transverse spin transport in a canted-antiferromagnet/heavy-metal heterostructure

We theoretically study the conditions under which an intrinsic spin Nernst effect–a transverse spin current induced by an applied temperature gradient–can occur in a canted-antiferromagnet insulator, such as LaFeO 3 and other materials of the same family. The spin Nernst effect may provide a microscopic mechanism for an experimentally observed anomalous thermovoltage in LaFeO 3 /Pt heterostructures, where spin is transferred across the insulator/metal interface when a temperature gradient is applied to LaFeO 3 parallel to the interface. We find that LaFeO 3 exhibits an intrinsic spin Nernst effect when inversion symmetry is broken on the axes parallel to both the applied temperature gradient and the direction of spin transport, which can result in a spin injection across the insulator/metal interface. Furthermore, our paper provides a general derivation of a symmetry-breaking-induced spin Nernst effect, which may open a path to engineering a finite spin Nernst effect in systems where it would otherwise not arise.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Arctic shrub root traits, northern Alaska, summer 2017 (Version 2.0)

This data package contains root trait data collected from 170 plots of rapidly expanding shrub genera (Alnus, Betula, and Salix) and a widespread sedge (Eriophorum vaginatum) along a latitudinal and temperature gradient in northern Alaska. The trait data were collected in July 2017 and include root architecture (root diameter and branching patterns), mycorrhizal colonization (%), nitrogen concentration (%), delta 15N (per mil), and vertical root biomass. These raw data support a submitted manuscript that examines the distribution and interspecific variations of absorptive root traits of shrubs and graminoids across the graminoid-dominated nutrient-poor arctic tundra and reveals how deciduous shrub expansion affects plant nutrient acquisition strategies in tundra ecosystems. Data are presented by site (n=5) and patch (shrub or sedge plot) in separate csv files. The location data are provided in the “plot coordinates” file; all other files contain the data in the file title. Detailed methods are in Chen et al. (2020).2023/07/20 Update: The latest version of this dataset publication is version 2.0. The latest version of the data package was updated to include the alder nodule biomass dataset (alder nodule biomass.csv) and associated metadata (metadata_alder nodule biomass.csv). The name of the previous metadata file was updated (metadata_ root traits and biomass.csv ) to distinguish it from the new metadata file.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Robustness of Deep Learning Classification to Adversarial Input on GPUs: Asynchronous Parallel Accumulation Is a Source of Vulnerability

The ability of machine learning (ML) classification models to resist small, targeted input perturbations—known as adversarial attacks—is a key measure of their safety and reliability. We show that floating-point non associativity (FPNA) coupled with asynchronous parallel programming on GPUs is sufficient to result in misclassification, without any perturbation to the input. Additionally, we show that this misclassification is particularly significant for inputs close to the decision boundary and that standard adversarial robustness results may be overestimated up to 4.6 when not considering machine-level details. We first study a linear classifier, before focusing on standard Graph Neural Network (GNN) architectures and datasets used in robustness assessments. We develop a novel black-box attack using Bayesian optimization to discover external workloads that can change the instruction scheduling which bias the output of reductions on GPUs and reliably lead to misclassification. Motivated by these results, we present a new learnable permutation (LP) gradient-based approach to learning floating-point operation orderings that lead to misclassifications. The LP approach provides a worst-case estimate in a computationally efficient manner, avoiding the need to run identical experiments tens of thousands of times over a potentially large set of possible GPU states or architectures. Finally, using instrumentation-based testing, we investigate parallel reduction ordering across different GPU architectures under external background workloads, when utilizing multi-GPU virtualization, and when applying power capping. Our results demonstrate that parallel reduction ordering varies significantly across architectures under the first two conditions, substantially increasing the search space required to fully test the effects of this parallel scheduler-based vulnerability. These results and the methods developed here can help to include machine-level considerations into adversarial robustness assessments, which can make a difference in safety and mission critical applications.

Shanmugavelu, Sanjif [Maxeler Technologies, a Groq↗
DOE OSTI2020

Final Technical Report for the SciDAC Center for the Study of Plasma Microturbulence

This final technical report for the SciDAC Center for the Study of Plasma Microturbulence documents contributions from MIT to the project over the period of the award (February 15, 2011 to February 14, 2020). This project focused on the development of neoclassical and gyrokinetic simulation capabilities relevant for the description of particle, momentum and energy transport in tokamak fusion plasmas. Extensive validation against experimental results from Alcator C-Mod and DIII-D experiments was performed. Overall areas of emphasis included gyrokinetic simulations, experiments, and observation of density gradient driven Trapped Electron Mode (TEM) turbulence relevant to burning plasmas; theory and gyrokinetic simulations of the parametric variation of the nonlinear upshift of the TEM critical density gradient; identification of a new mechanism for the isotope effect due to density gradient driven TEM turbulence; new synthetic diagnostics based on full-wave simulations which enabled unprecedented agreement between gyrokinetic simulations and measured fluctuation spectra with and without strong electron heating; development of new continuum pedestal neoclassical and turbulence codes; formulation and implementation of the first gyrokinetic exact Fokker-Planck collision operator in gyrokinetic turbulence simulations; development of new spectral approaches in velocity space for highly efficient treatments of collisions; new model collision operators implemented in the discontinuous Galerkin gyrokinetic edge turbulence code GKEYLL; and a new model and spectral code to study the interaction of ion- and electron-scale turbulence, which will serve as an efficient test-bed for multi-rate and multi-scale methods. These fundamental advances contribute to the development of a validated, predictive capability for turbulent and collisional particle, momentum and energy transport in both the core and edge of magnetically confined fusion plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Molecule-surface interaction from van der Waals-corrected semilocal density functionals: The example of thiophene on transition-metal surfaces

Semilocal density functional approximations are widely used. None of them can capture the long-range van der Waals (vdW) attraction between separated subsystems, but they differ remarkably in the extent to which they capture intermediate-range vdW effects responsible for equilibrium bonds between neighboring small closed-shell subsystems. The local density approximation (LDA) often overestimates this effect, while the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) underestimates it. The strongly constrained and appropriately normed (SCAN) meta-GGA often estimates it well. All of these semilocal functionals require an additive nonlocal correction such as the revised Vydrov-Van Voorhis 2010 (rVV10) to capture the long-range part. This work reports adsorption energies and the corresponding geometry of aromatic thiophene (C 4 H 4 S) bound to transition-metal surfaces. The adsorption process requires a genuine interplay of covalent and weak binding and requires a simultaneously accurate description of surface and adsorption energies with the correct prediction of the adsorption site. All these quantities must come from well-balanced short- and long-range correlation effects for a universally applicable method for weak interactions with chemical accuracy. Here, our methods indicate that the correct interplay is not present in any combination of recent meta-GGA's and rVV10. The simple short-range damping of the vdW correction scheme that is practically successful in the combination of GGA's and vdW approximations is less transferable in SCAN+rVV10 or in the revised version, revSCAN+rVV10. In addition, we present accurate random-phase-approximation-quality adsorption energies from a model based on the one of Zaremba and Kohn.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Collaborative Research: Natural Organic Matter and Microbial Controls on Mobilization/Immobilization of I and Pu in Soils and Waters Affected by Radionuclide Releases in USA and Japan

In this project, the relationship between natural organic matter (NOM) and two radioactive elements that are relevant to nuclear waste disposal were studied: Plutonium (Pu) and Iodine. The human and environmental risks associated with Pu stem mainly from the very long half-lives of several of its isotopes ( 238 Pu, 88 yr; 239 Pu, 24,100 yr; 240 Pu, 6560 yr) and its radiotoxicity. Understanding Pu biogeochemical behavior in both near-field (>10-11M) and far-field scales (<10-11M) is imperative to the development of approaches for reprocessing Pu, remediation of Pu contamination and accurate assessment of risks posed by disposal practices for Pu-bearing wastes. The environmental mobility of Pu can be affected by redox potential, pH, adsorption, precipitation, complexation, colloid formation, and microbial activity, of which the first characteristic has the most profound influence. Numerous studies have shown high affinity of Pu towards NOM, as well as to mineral phases. NOM is ubiquitous in the environment, e.g., both fulvic and humic acids are able to reduce Pu(V,VI) to Pu(IV) and the redox potential of NOM is positively related to the abundance of phenolic/acidic OH groups. NOM can either facilitate or limit actinide migration, depending on specific biogeochemical conditions including pH, mineral and organic matter characteristics, etc. The other radionuclide of interest is radioiodine ( 129 I). 129I is a major by-product of nuclear fission and of serious concern to the Department of Energy (DOE) as it is among the top risk drivers at existing and potential radiowaste-contaminated sites. The risk of 129 I stems largely from its high bioconcentration factor (90% of the body’s iodine is accumulated in the thyroid), a high inventory at source terms, a very long-half life (16M years), and rapid mobility in the subsurface environment. As a consequence, 129 I has the lowest drinking water standard (1 pCi/L) among all radionuclides in the Federal Register. With a novel and sensitive gas chromatography-mass spectrometry (GC-MS) method developed in our lab, it is possible to quickly and simultaneously determine the distribution of 129 I and stable 127 I forms in environments, as low as 2 pCi/L for 129 I. This method was subsequently validated using accelerator mass spectroscopy, AMS. IO 3 - and organo-I were determined as major species in the groundwater of SRS and the Hanford Site, contrary to thermodynamic predictions that I- should be the dominant species at these sites. Mobility of 129 I was also demonstrated to depend greatly on the I species and its concentration, sediment pH, and redox state, with times to achieve equilibrium taking up to 12 weeks. Along the groundwater pathway in the F-Area of SRS, 129 I- supplied from the seepage basins was transformed to 129 IO 3 - and organo- 129 I with increasing iodine sediment sorption, causing the lower total 127 I and 129 I concentrations along the gradient transect of the waste plume. By contrast, groundwater 129 I concentrations in the wetlands (as high as 1617.3 pCi/L) were greatly elevated with respect to the source term (159.3 pCi/L). While the NOM promoted the uptake of 129 I to the wetland sediment, it also promoted the formation of soluble organic fraction. A small fraction of NOM that is bound to iodine can behave as a mobile organo-I source. Iodide was enzymatically incorporated into NOM, whereas both iodide and iodate were abiotically bound to NOM, under certain conditions. Iodate removal from the mobile aqueous phase can also occur through incorporation into carbonate (e.g., at the Hanford Site, USA). Thus immobilization and re-mobilization of iodine species were influenced by pH, Eh and the presence of NOM and metal oxides, which adds to the complexity of site remediation action. A ground-breaking result was to elucidate the products (i.e. organo-iodine moieties formed via enzymatic and non-enzymatic processes) at the molecular level by nuclear magnetic resonance (NMR) and electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICRMS). We found that iodine-NOM interactions may be influenced by NOM hydrophobic aliphatic moieties. From the perspective of ESI-FTICRMS, organo-iodine formulas were ascribed to the groups of unsaturated hydrocarbons, lignins and proteins. Iodate is likely abiotically reduced to reactive iodine species by lignin- and tannin-like compounds or carboxylic-rich alicyclic molecules (CRAM). We also investigated microbial mechanisms in iodine incorporation into NOM. We established that soil bacteria isolated from F-Area of SRS did not accumulate significant amounts of I- (0.2-2%). Intracellular uptake of I- decreases with increasing pH when pH ranged from 4 to 6. In contrast, 44 out of 84 strains isolated from the F-Area of SRS can transform I- to IO 3 - and organo-iodine. In some cases, oxidation was facilitated in the presence of H 2 O 2 . Microbes can also excrete organic acids that enhance I- oxidation by lowering the ambient pH and reacting with H 2 O 2 to form peroxy carboxylic acids. At lower pH values (≤5), H 2 O 2 hydrolysis was the driving force for iodide-oxidation; whereas, at pH ≥ 6, spontaneous decomposition of peroxy carboxylic acids, originating from H 2 O 2 and organic acids were the primary cause of iodide oxidation. Lastly, it was determined that microbial processes involved in Mn (II) are capable of directly oxidizing I- via enzymatic catalysis (i.e., multicopper oxidases), or indirectly through the formation of reactive oxygen species (ROS) and/or biogenic manganese oxides. ROS-mediated oxidation of I- was found to predominate at pH >5, whereas the enzymatic and Mn oxide pathways were more active at pH < 5. Together, this project has resulted in 9 publications in high-impact journals, and the training of 1 Ph.D and 4 undergraduate students.

54 ENVIRONMENTAL SCIENCES↗