Diffusion and sorption of neptunium(V) in compacted montmorillonite: effects of carbonate and salinity

Exploring the kinetics of actinyl-EDTA reduction by ferrous iron using quantum-mechanical calculations

The reduction of An(vi) (An = U, Np, and Pu) to An(iv) significantly decreases its solubility and mobility. This reaction can be hindered by complexation with inorganic (e.g., carbonate) or organic ligands. Ethylenediaminetetraacetic acid (EDTA) is one such organic ligand that forms stable complexes with actinides. Therefore, it may enhance the mobility of actinides. However, the redox kinetics and mechanisms of actinyl (An(v/vi)O2+/2+)-EDTA are not well characterized yet and are thus studied here using quantum-mechanical calculations. The principle is to approach the actinyl-EDTA and Fe2+ (reductant) in small incremental steps and calculate the system energy at each distance. The overall reaction is then delineated into sub-processes (encounter frequency in bulk solution, formation of outer-sphere complex, transition from outer- to inner-sphere complex, and electron transfer), and reaction rates are determined for each sub-process. The formation of outer-sphere complexes occurs rapidly in microseconds to seconds over a wide range of actinyl concentrations (pM to μM); in contrast, the transition to the inner-sphere complex is relatively slow (milliseconds to a few seconds). Immediate electron transfer to form the pentavalent actinide is observed along the reaction path for Np(vi) and Pu(vi), but not for U(vi). Surprisingly, in acidic conditions, one of the carboxylic groups gets protonated in EDTA of [UO2(edta)]2- rather than one of the amino groups. This process-based series of calculations can be applied to any redox reaction and allows the prediction of changes to the rate law and rate-limiting step in a more fundamental way for different environments.

Decreasing Kd uncertainties through the application of thermodynamic sorption models

Radionuclide retardation processes during transport are expected to play an important role in the safety assessment of subsurface disposal facilities for radioactive waste. The linear distribution coefficient (Kd) is often used to represent radionuclide retention, because analytical solutions to the classic advection-diffusion-retardation equation under simple boundary conditions are readily obtainable, and because numerical implementation of this approach is relatively straightforward. For these reasons, the Kd approach lends itself to probabilistic calculations required by Performance Assessment (PA) calculations. However, it is widely recognised that Kd values derived from laboratory experiments generally have a narrow field of validity, and that the uncertainty of the Kd outside this field increases significantly. Mechanistic multicomponent geochemical simulators can be used to calculate Kd values under a wide range of conditions. This approach is powerful and flexible, but requires expert knowledge on the part of the user. The work presented in this paper aims to develop a simplified approach of estimating Kd values whose level of accuracy would be comparable with those obtained by fully-fledged geochemical simulators. The proposed approach consists of deriving simplified algebraic expressions by combining relevant mass action equations. This approach was applied to three distinct geochemical systems involving surface complexation and ion-exchange processes. Within bounds imposed by model simplifications, the presented approach allows radionuclide Kd values to be estimated as a function of key system-controlling parameters, such as the pH and mineralogy. This approach could be used by PA professionals to assess the impact of key geochemical parameters on the variability of radionuclide Kd values. Moreover, the presented approach could be relatively easily implemented in existing codes to represent the influence of temporal and spatial changes in geochemistry on Kd values.

Matrix diffusion and sorption of Cs+, Na+, I- and HTO in granodiorite: Laboratory-scale results and their extrapolation to the in situ condition

Matrix diffusion and sorption are important processes controlling radionuclide transport in crystalline rocks. Such processes are typically studied in the laboratory using borehole core samples however there is still much uncertainty on the changes to rock transport properties during coring and decompression. It is therefore important to show how such laboratory-based results compare with in situ conditions. This paper focuses on laboratory-scale mechanistic understanding and how this can be extrapolated to in situ conditions as part of the Long Term Diffusion (LTD) project at the Grimsel Test Site, Switzerland. Diffusion and sorption of (137)Cs(+), (22)Na(+), (125)I(-) and tritiated water (HTO) in Grimsel granodiorite were studied using through-diffusion and batch sorption experiments. Effective diffusivities (De) of these tracers showed typical cation excess and anion exclusion effects and their salinity dependence, although the extent of these effects varied due to the heterogeneous pore networks in the crystalline rock samples. Rock capacity factors (α) and distribution coefficients (Kd) for Cs(+) and Na(+) were found to be sensitive to porewater salinity. Through-diffusion experiments indicated dual depth profiles for Cs(+) and Na(+) which could be explained by a near-surface Kd increment. A microscopic analysis indicated that this is caused by high porosity and sorption capacities in disturbed biotite minerals on the surface of the samples. The Kd values derived from the dual profiles are likely to correspond to Kd dependence on the grain sizes of crushed samples in the batch sorption experiments. The results of the in situ LTD experiments were interpreted reasonably well by using transport parameters derived from laboratory data and extrapolating them to in situ conditions. These comparative experimental and modelling studies provided a way to extrapolate from laboratory scale to in situ condition. It is well known that the difference in porosity between laboratory and in situ conditions is a key factor to scale laboratory-derived De to in situ conditions. We also show that cation excess diffusion is likely to be a key mechanism in crystalline rocks and that high Kd in the disturbed surfaces is critically important to evaluate transport in both laboratory and in situ tests.

Spectroscopic identification of binary and ternary surface complexes of Np(V) on gibbsite

For the first time, detailed molecular information on the Np(V) sorption species on amorphous Al(OH)3 and crystalline gibbsite was obtained by in situ time-resolved Attenuated Total Reflection Fourier-Transform Infrared (ATR FT-IR) and Extended X-ray Absorption Fine Structure (EXAFS) spectroscopy. The results consistently demonstrate the formation of mononuclear inner sphere complexes of the NpO2(+) ion irrespective of the prevailing atmospheric condition. The impact of the presence of atmospheric equivalent added carbonate on the speciation in solution and on the surfaces becomes evident from vibrational data. While the 1:1 aqueous carbonato species (NpO2CO3(-)) was found to become predominant in the circumneutral pH range, it is most likely that this species is sorbed onto the gibbsite surface as a ternary inner sphere surface complex where the NpO2(+) moiety is directly coordinated to the functional groups of the gibbsite's surface. These findings are corroborated by results obtained from EXAFS spectroscopy providing further evidence for a bidentate coordination of the Np(V) ion on amorphous Al(OH)3. The identification of the Np(V) surface species on gibbsite constitutes a basic finding for a comprehensive description of the dissemination of neptunium in groundwater systems.

Influence of humic acid on neptunium(V) sorption and diffusion in Opalinus Clay

The influence of14C-labeled M42 humic acid (HA) on the interaction between neptunium(V) and natural clay rock (Opalinus Clay (OPA), Switzerland) has been investigated in batch sorption and diffusion experiments under ambient air conditions. The effect of 10 mg/L HA on the diffusion of 8 μM Np(V) in OPA has been investigated in synthetic OPA pore water (pH 7.6, I = 0.4M) for the first time. Batch sorption experiments as a function of solid-to-liquid ratio (4-20 g/L) were performed under same experimental conditions to compare distribution coefficients obtained from both diffusion and sorption experiments. These experiments showed only a slight influence of HA on Np(V) uptake by OPA in both cases and provided comparable distribution coefficients (presence of HA: Kd= 22-32 L/kg, absence of HA: Kd= 30-46 L/kg). As it is known that the interactions of humic substances with actinides depend on various experimental parameters, the effect of HA on Np(V) sorption on 15 g/L OPA was also investigated as a function of pH (6-10) and initial Np concentration (8 μM and 7 pM). A saturated calcite solution was used as a background electrolyte in this case to prevent any dissolution of calcite contained in OPA at low pH. The results showed that the presence of M42 HA increases Np(V) sorption at pH < 7 while lower sorption was obtained above pH 8. Higher sorption at acidic pH can be attributed to strong sorption of HA which increases the number of sorption sites, while sorption decreases at higher pH probably due to formation of soluble ternary neptunium humate carbonate species.