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

Abstract

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)O₂^+/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 Fe²⁺ (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 [UO₂(edta)]²⁻ 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.