Plan
Comptes Rendus

Estimating the stabilities of actinide aqueous species. Influence of sulfoxy-anions on uranium(IV) geochemistry and discussion of Pa(V) first hydrolysis
Comptes Rendus. Chimie, Volume 10 (2007) no. 10-11, pp. 978-993.

Résumés

Qualitative chemical information is used as a guideline for correlations between equilibrium constants or between equilibrium constants and atomic charges (deduced from quantum mechanics calculations). Pa(V) and Nb(V) hydrolysis constants are also recalculated from experimental data. logK1°(AnIV/RO22)=6.59±0.55(S2O32),10.06±0.88(SO32),11.97±1.07(CO32), and10.05±0.88(HPO42) are estimated based on the trend of affinity for An cations in the series CO32>HPO42SO32>SO42S2O32. These ideas and values are used to discuss U(IV) chemistry in S-containing ground-waters.

Des connaissances qualitatives ont été concrétisées sous la forme de corrélations empiriques entre constantes d'équilibres, voire avec les charges atomiques (issues de calculs quantiques) dans la série CO32>HPO42SO32>SO42S2O32, pour, par exemple, estimer logK1°(AnIV/RO22)=6.59±0.55(S2O32),10.06±0.88(SO32),11.97±1.07(CO32), et10.05±0.88(HPO42). Ces valeurs sont utilisées pour prévoir l'influence éventuelle d'anions soufrés sur la chimie de U(IV) dans des eaux souterraines.

Métadonnées
Reçu le :
Accepté le :
Publié le :
DOI : 10.1016/j.crci.2007.04.015
Keywords: Actinides, Hydrolysis, Sulphate, Sulfite, Thiosulphate, Carbonate, Groundwaters
Mots clés : Actinides, Hydrolyse, Sulfate, Sulfite, Thiosulfate, Carbonate, Eaux souterraines

Pierre Vitorge 1, 2 ; Vannapha Phrommavanh 3 ; Bertrand Siboulet 4 ; Dominique You 1 ; Thomas Vercouter 1 ; Michael Descostes 2, 3 ; Colin J. Marsden 5 ; Catherine Beaucaire 3 ; Jean-Paul Gaudet 6

1 Laboratoire de spéciation des radionucléides et des molécules (LSRM), CEA Saclay, DEN/DANS/DPC/SECR, 91191 Gif-sur-Yvette cedex, France
2 UMR 8587, CEA Saclay, DEN/DANS/DPC/SECR, 91191 Gif-sur-Yvette cedex, France
3 Laboratoire de mesures et modélisation de la migration des radionucléides (L3MR), CEA Saclay, DEN/DANS/DPC/SECR, 91191 Gif-sur-Yvette cedex, France
4 CEA Marcoule, DEN/DRCP/SCPS, 30207 Bagnols-sur-Cèze cedex, France
5 Laboratoire de physique quantique, CNRS–UMR 5626, Université Paul-Sabatier, 118, route de Narbonne, 31062 Toulouse cedex 4, France
6 Laboratoire d'étude des transferts en hydrologie et environnement (LTHE), UMR 5564, CNRS/INPG/IRD/UJF, BP 53, 38041 Grenoble cedex 9, France
@article{CRCHIM_2007__10_10-11_978_0,
     author = {Pierre Vitorge and Vannapha Phrommavanh and Bertrand Siboulet and Dominique You and Thomas Vercouter and Michael Descostes and Colin J. Marsden and Catherine Beaucaire and Jean-Paul Gaudet},
     title = {Estimating the stabilities of actinide aqueous species. {Influence} of sulfoxy-anions on {uranium(IV)} geochemistry and discussion of {Pa(V)} first hydrolysis},
     journal = {Comptes Rendus. Chimie},
     pages = {978--993},
     publisher = {Elsevier},
     volume = {10},
     number = {10-11},
     year = {2007},
     doi = {10.1016/j.crci.2007.04.015},
     language = {en},
}
TY  - JOUR
AU  - Pierre Vitorge
AU  - Vannapha Phrommavanh
AU  - Bertrand Siboulet
AU  - Dominique You
AU  - Thomas Vercouter
AU  - Michael Descostes
AU  - Colin J. Marsden
AU  - Catherine Beaucaire
AU  - Jean-Paul Gaudet
TI  - Estimating the stabilities of actinide aqueous species. Influence of sulfoxy-anions on uranium(IV) geochemistry and discussion of Pa(V) first hydrolysis
JO  - Comptes Rendus. Chimie
PY  - 2007
SP  - 978
EP  - 993
VL  - 10
IS  - 10-11
PB  - Elsevier
DO  - 10.1016/j.crci.2007.04.015
LA  - en
ID  - CRCHIM_2007__10_10-11_978_0
ER  - 
%0 Journal Article
%A Pierre Vitorge
%A Vannapha Phrommavanh
%A Bertrand Siboulet
%A Dominique You
%A Thomas Vercouter
%A Michael Descostes
%A Colin J. Marsden
%A Catherine Beaucaire
%A Jean-Paul Gaudet
%T Estimating the stabilities of actinide aqueous species. Influence of sulfoxy-anions on uranium(IV) geochemistry and discussion of Pa(V) first hydrolysis
%J Comptes Rendus. Chimie
%D 2007
%P 978-993
%V 10
%N 10-11
%I Elsevier
%R 10.1016/j.crci.2007.04.015
%G en
%F CRCHIM_2007__10_10-11_978_0
Pierre Vitorge; Vannapha Phrommavanh; Bertrand Siboulet; Dominique You; Thomas Vercouter; Michael Descostes; Colin J. Marsden; Catherine Beaucaire; Jean-Paul Gaudet. Estimating the stabilities of actinide aqueous species. Influence of sulfoxy-anions on uranium(IV) geochemistry and discussion of Pa(V) first hydrolysis. Comptes Rendus. Chimie, Volume 10 (2007) no. 10-11, pp. 978-993. doi : 10.1016/j.crci.2007.04.015. https://comptes-rendus.academie-sciences.fr/chimie/articles/10.1016/j.crci.2007.04.015/

Version originale du texte intégral

1 Introduction

The geochemical behaviour of actinides have been extensively studied for understanding uranium and thorium ore deposits, and more recently for assessing the environmental impact of the possible disposal of wastes that contain Pa, Np, Pu, Am or Cm. Starting with uranium [1], the Thermochemical Data Base project of the Nuclear Energy Agency (NEA-TDB) organised the reviewing of published experimental data relevant for modelling aqueous chemistries and solubilities of the most important radionuclides. The results of the NEA-TDB project are now well accepted as a reference critical review, essentially for aqueous chemistry and solubility at room temperature, but these reviews proposed data only when convincing experimental validations were published. There is therefore a gap between this restricted set of quantitative validated thermochemical data and qualitative chemical knowledge.

However, besides selected numerical values, the NEA-TDB reviews also provide qualitative information [2,3], which we used together with analogies for estimating the hydrolysis constants and standard potentials [4] needed for drawing Pourbaix' diagrams of actinides [5]. The present paper aims at testing such rough estimates for complexation.

Rules of thumb are currently used by chemists for checking the possible formation of hypothetical chemical species in specific chemical conditions, when these chemical species are not in databases. This might typically be the case for the environmental aqueous chemistry of An(IV), the actinide elements at the +4 oxidation state (An = Th, U, Np and Pu) in the presence of S-containing inorganic ligands; for this reason it is also an aim of the present paper to estimate the stabilities of An(IV) complexes with sulfoxy-anions.

For storing radioactive wastes, several projects are looking for geological sites that are well isolated from surface waters. These often correspond to anoxic conditions, where U, Np and Pu are expected to be stable in the +4 oxidation state [2,4]. Interestingly, chemical analogues are Ce(IV) and Th, and probably Zr and Hf. Selected NEA-TDB equilibrium constants and redox potentials [2] are adequate for reliable modelling of uranium chemistry in most equilibrated groundwaters. Uranium is predicted to be stable in anoxic waters in the form of the U(OH)4(aq) aqueous species in equilibrium with uraninite, UO2(s), a compound of low solubility [2,5,6]. Similar behaviour is expected for Np and Pu, even though Pu3+ might also be stable [4–7].

In natural under-ground-waters CO32, the carbonate anion is often the dominating ligand among the inorganic ligands (for actinide cations). However, carbonate complexes are predicted to be of little important for An(IV) [2,4,7]. Nevertheless, many An(CO3)i(OH)j42ij complexes have been proposed, but no reliable values could be validated for most of the corresponding formation constants [2,4], for which maximum possible values have been estimated from experimental observations, that also confirmed the similar behaviour of Th, U(IV), Np(IV) and Pu(IV) in carbonate/bicarbonate aqueous solutions [4,5,8]. Values have recently been proposed for the formation constants of several Th(CO3)i(OH)j42ij complexes in an attempt to interpret a solubility study of ThO2(s) [8,9], but the Th(CO3)44 species were not included in the interpretation, even though we shall see it should not be completely negligible according to the NEA-TDB data [4]. For probing such competition between the HO and CO32 ligands, the relative stabilities of the corresponding 1:1 complexes of An4+ will be compared. In this framework, the stability of AnCO32+ will be estimated: its existence has never been demonstrated, since it is always hidden by hydrolysed species. This over-stabilisation of hydrolysis is specific to the +4 oxidation state: the 1:1 carbonate complexes are well known, and their stabilities were well established for An(III), An(V) and An(VI) (it is not known for Pa(V), which is known to have a very different aqueous chemical behaviour from the other An(V)).

Several other inorganic hard anions are quite reactive toward actinide cations, specially those of high charge (PO43), of small size (F), or polydentade (HPO42): depending on their content in groundwaters they might form complexes with actinide cations. For example, it was recently proposed that the pore-waters of Callovo–Oxfordian clay minerals, where an underground research laboratory is being built to study the feasibility of a deep geological repository for radioactive wastes in France [10], contain quite large amounts of SO42. The pore-water composition is approximately at the SO42/HS/H2S frontier point (Fig. 1). Although the most stable aqueous sulfur species - those bolded on Pourbaix' diagrams - are sulfide (HS) and sulfate (SO42), some other sulfoxy-anions are often detected in natural environments, typically as thiosulfate (S2O32) and sulfite (SO32) ions. For example, about 25% of the total S content has been reported to be S2O32 in a reducing groundwater [11], and even in more oxidizing conditions [12]. Sulfoxy-anions are also suspected to form in the course of the oxidative dissolution of pyrite (FeS2) [13], a mineral often associated with the redox regulation of ground-waters. For this reason, we focus on sulfoxy-anion ligands, specially on S2O32.

Fig. 1

Pourbaix diagram of sulfur. (a) The names of the major species (HS, H2S, SO42,HSO4 and solid S) are in bold, and their predominance domains are drawn with thick lines. The domains (thin dashed lines) of the other (minor sulfur) species are obtained by suppressing SO42,HSO4 and H2SO4. (b and c) The speciation (namelylog[HiSjOkz]/[S]total) of each species (HiSjOkz) is represented for two kinetic assumptions (see text) in redox conditions corresponding to line B of Fig. 6.

Grenthe et al. have selected complexing constants for U(VI) complexes with these anions [2], but they wrote in their review that confirmation is needed: “The only quantitative information about aqueous uranium thiosulfate complexes is the study by Melton and Amis […] This review tentatively accepts [their] value, (although confirmation of the results from another study would be useful)”. Furthermore: “The solid formed seems to be a mixture indicating decomposition of thiosulfate into sulfite and elemental sulfur. This review finds no reliable evidence for the formation of solid uranium thiosulfate compounds.” This decomposition might very well be the result of redox reactions (disproportionation), since both uranium and sulfur have a wide range of possible oxidation states, and their stability domains are not well established in mixtures of uranium and sulfur. This might in fact be a problem for other uranium compounds and complexes with sulfur-containing ligands. Unfortunately the NEA-TDB reviews could not validate such data for the Th and Zr analogues [14,15]. Since the NEA-TDB review has selected data for a thiosulfate complex of U(VI) [2], one might very well expect thiosulfate complexes of U(IV) in more reducing conditions – unless thiosulfate is strongly reduced, when U(IV) is formed – because the U4+ hard cation is usually more reactive (than UO22+) toward oxygen-donor ligands. Based on the same hardness rule, thiosulfate should bind to hard cations via the O rather than the S atom of the sulfoxy-anion ligands. The same problem holds for the U/SO32 system: the NEA-TDB review has selected data for sulfite complexes of U(VI); but not for U(IV): “Formation of aqueous uranium(IV) sulfite complexes was reported in a qualitative study by Rosenheim and Kelmy. However, no experimental chemical thermodynamic data on these species are available.” As a probe for ligand competition (between S2O32orSO32 and typically OH or CO32), we shall estimate the formation constants of the corresponding 1:1 An(IV) complexes, namely AnS2O32+andAnSO32+.

Various methods are commonly used for estimating equilibrium constants, typically as empirical correlations with physical (or phenomenological) parameters (atomic radii and charges, solvent interactions…). They can also be obtained from molecular modelling methods. Indeed we recently estimated an uncertainty of 10 kJ mol−1 on ΔrG for Pa(V) hydrolysis [16], which is about 20 times higher than the uncertainty of the experimental determinations of equilibrium constants and standard redox potentials in aqueous solutions. Furthermore, such methods need caution, when using calculated energies [17]. After others [18,19] (with a comment in Ref. [7]) we also tested empirical correlations, which appear to work surprisingly well for some of them: they actually help in putting numbers for quite encyclopedic qualitative knowledge. It is also a way to check such knowledge and corresponding chemical intuition. We use such correlations here, hoping this special issue will help such rules of thumb used by actinide chemistry specialists to become less mysterious. The rough estimates are only guidelines, which need experimental confirmation. They often originate in geometrical and electrostatics reasoning. Indeed, the chemical stabilities of hard cations are often correlated with the charge2/radius ratios of the reactants. Nevertheless, we shall see that the best correlations are not specially with the atomic charges; correlations between measured equilibrium constants will often appear to fit better: several physical contributions probably cancel out in those cases.

The present paper is organised as follows. First, definitions are given, together with features and explanations of methodologies from the NEA-TDB reviews. Results are first reported for correlations between U(VI) (or analogous An(VI)) 1:1 (1 cation with 1 ligand) complexes and protonation of the corresponding ligands, two types of reactions for which many published data are available. In a next step, such correlations are extended to An(III) and An(V). There are virtually no reliable published 1:1 complexing constants of An(IV), since An4+ hydrolysis usually overcomes complexation. Consequently we shall estimate An(IV) 1:1 complexing constants. We shall also consider the position of Pa(V) in our correlations, since its chemical behaviour is an exception as compared to that of the other An(V): this comparison will be based only on hydrolysis data, since there are very few other published equilibrium constants for Pa(V) aqueous complexes. Prior to this comparison, we shall examine the impact of our estimated complexing constants on the geochemical behaviour of U(IV).

2 Methods

2.1 Equilibrium constants

For consistency, we used reaction data – i.e. ΔrG° (equivalently standard equilibrium constants and potentials of redox couples) – they were preferred to formation data (ΔfG°) [2].

K1°=|MLzM+zL||MzM||LzL|(1)
where |A| is the activity of species A.
Ka°=|H+||LzL||HL1+zL|(2)
pKa°=logKa° is pH1/2, the pH at the half-point reaction (where |LzL|=|HL1+zL|). pKa° appears to be the ionic product of water, when LzL=OH and when using |H2O| = 1; but when comparing the complexing strengths of various ligands we used the concentration of liquid water: |H2O|=CH2O (55.34 mol kg−1). Superscript ° stands for infinite dilution (ionic strength, I = 0), the standard conditions (see Section 2.2).

logK1° values were plotted as a function of pKa°, for example when K1° is the formation constant of AmCO32+, 1/Ka° is the protonation constant of CO32. We observed

logK1,x°=ax+bxpKa°(3)
linear correlations, where we fitted the ax and bx parameters for a series of actinides (or analogues) of the same oxidation number X: An(X)=MzM=An3+,An4+,AnO2+orAnO22+. 10a=K1°(Ka°)b=|MLzM+zL||H+|b/(|MzM||LzL|1b|HL1+zL|b) is obtained from Eqs. (1)–(3): it is actually an equilibrium constant. When b = 1, it simplifies into the 10a(1)=K1°Ka°=|MLzM+zL||H+|/(|MzM||HL1+zL|), the
MzM+HL1+zLMLzM+zL+H+(4)
exchange equilibrium constant: |LzL| cancels out in K1°Ka° (=10a(1)), which is now the same for all the reactions; it also only depends on a, not on each Ka° (or equivalently K1°). In that case, the half-point reaction definition is [|H+|/|MzM|]1/2=K1°Ka°: a solution of activity |H+| has the same reactivity (for the LzL ligand) as a solution of activity |MzM|10a(1). Conversely, when the (b) slope is not 1, this simple equivalent solution definition is no longer relevant: the general half-point reaction definition is 10a=K1°(Ka°)b=[(|H+||LzL|)b/(|MzM||LzL|)]1/2, whose interpretation is less intuitive. Furthermore the constant of the exchange equilibrium (Eq. (4)) is now specific for each reaction: K1°Ka°=10a(Ka°)1b=10a/b(K1°)11/b depends on each Ka° (or equivalently K1°), not only on a and b.

Since bx did not seem to strongly depend on X, we finally used the same bx (=b) value for all the oxidation states; in that case, for comparing An(X) complexes with An(VI) ones, we used:

10aXaVI=K1,X°K1,VI°=KX/VI°=|MLzM+zL||AnO22+||AnO2L2+zL||MzM|(5)
the constant of the
AnO2L2+zL+MzMMLzM+zL+AnO22+(6)
exchange equilibrium. This AnO22+/MzM exchange equilibrium is similar to the H+/MzM one (Eq. (4)), where now the half-point reaction definition is [|AnO22+|/|MzM|]1/2=Kx/VI°: a solution of activity |AnO22+| has the same reactivity (for the LzL ligand) as a solution of activity |MzM|Kx/VI°. Similarly, for comparing hydrolysis equilibria, we used the
AnO2OH++MzMMOHzM1+AnO22+(7)
hydrolysis competition equilibrium of constant
K1,x°K1,VI°=KX/VI°=|MOHzM1||AnO22+||AnO2OH+||MzM|(8)
where
Ki,X°=|M(OH)izMi||H+||M(OH)i1zMi+1|(9)
is a classical stepwise standard equilibrium constant. logKi,x° is pH1/2.

2.2 Activity coefficients

For ionic strength corrections, we used γi, the molal activity coefficient of ion i calculated with the “Specific Interaction Theory”, the SIT formula [2,4]. The corresponding ɛij empirical (pair interaction) coefficients are taken from the NEA-TDB reviews [4]:

logγi=zi2D+jɛijmj(10)
In most cases, the summation could here be restricted to the ClO4 and Na+ (dominating) counter-ions (= j). mj is j (molal) concentration (mol per kg of pure water). Molar (M = mol L−1) to molal conversion coefficients are tabulated in Handbooks (including the cited NEA-TDB books):
D=AIm1+BrIm(11)
is a Debye–Hückel term, where Im is molal I, A=0.509[298ɛ298/Tɛr,T]1.5[d298/dT]0.5kg1/2mol1/2, and B = 3.28 × 109 [298ɛ298dT/(Tɛr,Td298)]0.5kg1/2mol1/2m1 are calculated from physical constants, and ɛr,T the relative dielectric constant of the solvent (water) and dT the density at absolute temperature T. r accounts for geometric exclusion about ion i. r is assumed to be constant, B r = 1.5 kg1/2 mol−1/2 at 25 °C, an approximation that enables many measured mean activity coefficients of strong inorganic aqueous electrolytes to be fitted [2]; it corresponds to r = 4.57 × 10−10 m, which appears to be of the correct order of magnitude for most inorganic hydrated ions including their first hydration sphere (by definition if a counter-ion stays in the first hydration sphere, it is treated as complex formation, not non-ideality). Furthermore, BI=1/lD, where lD is the Debye distance, the distance between an ion and its counter-ion atmosphere: 3.05 × 10−10 m (at 25 °C and I = 1 M) in the SIT approximation. The Debye–Hückel formula is valid for large lD namely for lD >>r: B r = 1.5 kg1/2 mol−1/2 corresponds to r/lD=1.5Im<<1. Indeed the Debye–Hückel formula (D term alone) is only valid at ionic strengths less than 10 or 1 mM, while the SIT formula is usually a good approximation for aqueous solutions with ionic strength up to 4 mol kg−1 [2]. The use of ɛij empirical coefficients certainly partly compensates systematic errors (of D in the SIT formula), which might very well explain why the numerical values for ɛij are correlated with the size and the charges of the ions (Fig. 2).

Fig. 2

Empirical correlations for estimating SIT ion pair coefficients.

Unknown ɛij numerical values can be estimated by analogy with ions of same charge and similar sizes [3,4]. In that case it was proposed to increase the ɛij uncertainties by ±0.05 kg mol−1 [4]. Moreover, we observed a reasonable linear correlation between interaction coefficients and the charge/radius ratios (Fig. 2(a)). For this correlation we used the formal charge of the cation (or the complex), which is indeed the charge seen by the (ClO4) counter-anion at large distance. Unfortunately this assumption does not hold for high ionic strengths, the only conditions where the ɛij mj term cannot be neglected. Indeed some data for AnO2X4 cations are not in the middle of the correlation cloud, indicating that the relevant phenomenological charge might be higher (than the formal one). Conversely, the corresponding points are moved to the other side of this correlation plot (not represented on the figure), when using the atomic charge (Section 2.3) of An (in AnO2X4, for this reason, we also plotted other correlations without using the charges (Fig. 2(b) and (c))).

2.3 Quantum calculations

In our hydrolysis correlation study, we used the formal charges of An3+ (3) and An4+ (4), while we used the atomic charge of An(X) in AnO2X4 (X = 5 or 6). This latter was deduced from quantum (DFT) calculations performed at the same calculation levels (ECP and basis sets) as in our recent previous works [17,20], from which we extracted NPA atomic charges [21,22]. Note that in Gaussian98 and 03, these NPA charges are calculated with the NBO software, which is known to consider the UO22+ 6d orbitals as Rydberg orbitals, despite the final result gives the 7s5f6d electronic configuration. This overestimates the U atomic charge by 0.50 electron. For this reason we recalculated the NPA charge assuming the 6d to be valence orbitals for all the oxy actinide cations.

For the ligands (alone) closed-shell ab initio calculations were performed at the MP2/6-311+g(2df,2p) level; open-shell calculations are not needed, as checked at the B3LYP/6-31+g(d,p) level. All the quantum calculations were done with the Gaussian98 and 03 suites of programs [23,24].

ΔrE, the ab initio energy, was calculated for the protonation reaction corresponding to the pKa equilibrium (Fig. 3). pKa represents pH1/2; similarly, when adding the entropic contribution to ΔrE, it represents PH1/2+, the H+ partial pressure: the hydration energy of H+ is related to the slope in the ΔrE plot as a function of pKa (Fig. 3), while the intercept rather reflects the hydration energy change between the reactants and products. Since the slope is much greater than 1, the E scale had to be contracted, which means that the final correlation is quite poor. Slope ≠ 1 also means it is not equivalent to use pKa or ΔrE for our estimates of complexing constants: pKa appeared to be better correlated with measured complexing constants. Furthermore pKa represents a marginal value of ΔrE, which already is a marginal part of the DFT calculated electronic energy: the experimental pKa values are more reliable for our purpose and more accurate than the ΔrE ones obtained from quantum calculations in the gas phase.

Fig. 3

Protonation reactions of RO22 ligands in the gas phase and in aqueous solutions. The pKa (in pure liquid water at 25 °C) are converted to kJ mol−1rG = −RTpKa) and compared to the energy of the corresponding reaction ab initio calculated (with neither matrix nor temperature correction), where the protonation is on an O (black points) or S (white points) atom.

3 Results and discussions

3.1 Comparing the affinities of hard anions for H+ and AnO22+

We first examined U(VI) aqueous complexation and hydrolysis data, since a sufficiently large set of complexing constants is available. Positive logK1°vspKa° correlations are observed, and even linear correlations appear to fit the data reasonably well. The plot appeared to be less scattered when restricted to the consistent set of data selected by the NEA-TDB reviews (Fig. 4): we finally used these (NEA-TDB) consistent sets of data for the figures and numerical correlations given in the present paper with the following exceptions. For An(VI) we did not use the NEA-TDB PuO2CO3(aq) formation constant based on only (logK1°=13.80.6+0.8 [26] and 9.3 ± 0.5 [27]) two experimental determinations that are not consistent (within uncertainties). The latter value (that we published ourselves) is closer to the U(VI) and Np(VI) ones (Table 1) [2,4], but it is rather a maximum possible value corresponding to the detection limit of our solubility measurements. For this reason, we also do not rely on the 9.5 ± 0.5 similar value more recently updated by the NEA-TDB review [16].

Fig. 4

Proton and other cation affinities for ligands. K1° is the standard formation constant of the MLzM+zL complex, and Ka° is the protonation constant of the LzL ligand written on the figure (Table 1). Since all the (logK1,X°=aX+bXpKa°) lines are virtually parallel for all the (X) oxidation states (a). We also used the same (bx = 0.62 ± 0.16) slope in logK1,XaX°=((0.62±0.16)pKa°) regressions for the RO22 ligands, and shifted the curves by ax (b): aXaVI=(logK1,X°logK1,VI°)=logKX/VI° (Eq. (5)), where KX/VI° is the constant for the AnO2RO2/MRO2zM2 exchange equilibrium (Eq. (6a)) for MzM = An3+, An4+ or AnO2+ For An at the oxidation state X,aX(=logK1,X°0.62pKa°) can be used as a definition of a quantitative scale for the (up to now qualitative) An4+>AnO22+An3+>AnO2+ series (y axis of (c)). It is compared (black filled symbols) with the hydrolysis constant (*K1) scale on the x axis (of (c)), and zAn, the An atomic charge (blue open symbols corresponding to the top scale). zO, the atomic charge of O in RO22 is not specially correlated with its pKa (d).The x-axes are the same for (a), (b) and (d): the names of the ligands are only written on (a).

Table 1

Formation constants used to draw Fig. 4: pKa° of RO22 ligands and corresponding K1,X° formation constantsc with cations of f-block elements [2–4,16]

RO22pKablogK1,III° clogK1,IV°clogK1,V°clogK1,VI°c
S2O321.596.6 ± 0.8a2.8(U)
6.58(U)3.15(U)
SO421.983.85(Am)6.85(Np)0.44(Np)3.28(Np)
3.91(Pu)6.89(Pu)3.38(Pu)
HPO427.2129.5 ± 2.3a2.95(Np)7.24(U)
6.2(Np)
SO327.229.5 ± 2.3a6.6(U)
CO3210.3297.7 ± 0.3(Am)11.1 ± 3.2a4.962(Np)9.67(U)
7.8 ± 0.2(Eu)5.12(Pu)9.32(Np)
5.1(Am)
ad2.25.8 ± 0.4−0.72.0 ± 0.3
bd0.550.51 ± 0.270.540.665 ± 0.085
axe2.253 ± 0.0755.444 ± 0.249−1.236 ± 0.4892.408 ± 0.453
logKX/VI°−0.16 ± 0.53.04 ± 0.5−3.64 ± 0.50
logKX/VI°−1.74.81−6.20

a Value estimated in the present work.

b pKa of the RO22 ligand.

c K1,X° is the standard constant of equilibrium Mz++RO22MRO2(z2)+ for Mz+ = Anx+ (X = 3 or 4) or AnO2(X4)+ (X = 5 or 6) from published data (see text) [2–4,7].

d a and b, the coefficients of the logK1,X°=(a+bpKa°) linear regression are fitted for a given oxidation state, X.

e while ax is the (fitted) intercept of a similar regression, but with the same (0.62 ± 0.16) slope for all the oxidation states: logK1,X°ax=((0.62±0.16)pKa°).

The PO43 tri-anion is out of the correlation, while the F mono anion could probably be included into the correlation, but HO cannot (Fig. 4). When restricting the correlations to LzL=RO22 potentially bidentate oxygen-donor ligands, virtually the same logK1° value is observed for a given ligand with AnO22+ for An = U, Np, Pu or Am: the variations along this series are only slightly higher than experimental uncertainties. logK1,VI°=((2.0±0.3)+(0.665±0.085)pKa°) is obtained for these (AnO22+) cations. Adding not critically reviewed (logK1,VI°) values gives similar regression coefficients with increased uncertainties, namely (2.6 ± 0.8) instead of (2.0 ± 0.3), and (0.58 ± 0.08) instead of (0.665 ± 0.085).

3.2 AnO2RO2/MRO2zM2 exchanges for MzM = An3+, An4+ or AnO2+

Similar logK1°vspKa° correlations and observations (as for AnO22+, Section 3.1) are made for the An3+ and AnO2+, cations, while there are too few data for M4+. For the M3+/CO32 systems, we used our own complexing constants (Table 1) [28], AmHCO3+ is an outlier, in part due to the stabilisation of H2CO3 as CO2(aq): the measured pKa (of the HCO3/CO2(aq) couple) is not the relevant parameter for our correlations. For LzL=RO22 we found K1° values in the order An4+>AnO22+An3+>AnO2+, a classical order for the reactivity of actinide cations toward hard anions. However, AnO22+An3+ is often written (instead of AnO22+An3+). Here the available data for An3+ appear to be within the correlation lines of AnO22+. The atomic charge of U in UO22+ is 2.8 electron, which compares with the charge of Am3+. Nevertheless, using atomic charges in the correlations gave poorer results (see below and Fig. 4(c)). We obtained logK1,III°2.2+0.55pKa° for An3+ and logK1,v°0.7+0.54pKa° for AnO22+.

There are too few data for a statistical evaluation of all the uncertainties. Nevertheless, it seems that bx, the slopes (of every logK1°vs pKa° correlations) are the same within uncertainties. Indeed a reasonable fit is obtained when fixing bx = bVI (=0.665), namely logK1,X°=(aX+0.665pKa°) for all the oxidation states. We finally fitted the slope (and ax) on all the data, and obtained logK1,X°=(aX+(0.616±0.055)pKa°) with aIII = 2.253 ± 0.075, aIV = 5.444 ± 0.249, aV = −1.236 ± 0.489 and aVI = 2.408 ± 0.453 fitted values, where uncertainties are 1.96σ (the maximum errors on logK1,X° was +0.9 for UO2CO3 and −0.8 for AnCO3+). Since there are not enough data for a meaningful statistical analysis, we finally increased the uncertainties: logK1,X°aX=(0.62±0.16)pKa° (Fig. 4(b)). The standard deviation of the fit is 0.5 (it represents less than 3 kJ mol−1), a value only a little higher than that of many direct experimental determinations in aqueous solutions, and of the order of magnitude of the variations of the logK1° values among the series of the analogous cations considered here. However, when a slope is fitted for each oxidation state, a small systematic deviation can be inferred, namely the slopes seem slightly different for each type of central cation (AnO22+, An3+ and AnO2+). On the other hand, this difference is small as compared to the scatter in the available data, and it is not specially correlated with the charge of the cation.

Since the (0.62 ± 0.16) slope does not depend on X, (aXaVI)=logKX,VI° (Eq. (5)), the shift between the lines in Fig. 4(a) is related to the following equation:

AnO2RO2+MzMMRO2zM2+AnO22+(6a)
AnO2RO2/MRO2zM2 exchange equilibrium (Eq. (5)) whose equilibrium constant is KX/VI°. For such anions the order of their reactivity toward H+ and actinide ions is found to be
CO32>HPO42SO32>SO42S2O32
also corresponding to their pKa° values, but not specially to zO, the atomic charge of O (Fig. 4(d)).

We also attempted to force b = 1 in the fits, but without success: the correlations cannot simply be interpreted with Eq. (4).

Similarly, hydrolysis equilibria were compared using KX/VI° (Eq. (8)), the constant of a hydrolysis competition equilibrium (Eq. (7)) between actinide ions at oxidation states X and +6: logKIII/VI°=1.7, logKIV/VI°=4.81 and logKV/VI°=6.2 for Np [4]. These results give nearly the same An4+>>AnO22+>An3+>>AnO2+ qualitative scale (Fig. 4(c)), and the logKX/VI°=(0.36+0.60)logKX/VI° linear correlation.

3.3 Estimating the stabilities of An(IV)RO22+ complexes

We now want to estimate missing complexing constants for U(IV), namely with the SO32 and S2O32 anions. For this we could not use the same type of correlations as those observed for UO22+ and the other cations (Section 3.2), because reliable formation constant has been published for only one type of 1:1 An(IV) complexes, namely for the MSO42+ complexes: logK1,IV°=6.58(USO42+) [2], 6.85 (NpSO42+) [4], 6.89 (PuSO42+) [4] and 7.04 (ZrSO42+) [15]. When using this single datum and the (0.62 ± 0.16) slope value estimated above, the (logK1,IV°=((5.6±0.2)+(0.62±0.16)pKa°) line is drawn, from which logK1,IV°=(12.0±1.9), (10.1 ± 1.4), (6.8 ± 0.5) and (6.6 ± 0.5) are calculated for the An(IV) complexes of CO32, SO32,SO42 and S2O32, respectively. The value of aIV = 5.6 ± 0.2 was chosen to fit the 6.58 (USO42+), 6.85 (NpSO42+) and 6.89 (PuSO42+) data by giving less weight to USO42+, because the +4 oxidation state is more difficult to stabilise for uranium, even if more experimental studies of U(IV) are available. However, the logK1,IV° (=6.6 ± 0.5) MS2O32+ value is close to the MSO42+ fixed point (of the correlation). Therefore it does not depend sensitively on the value estimated for the slope (of the logK1,IV°vspKa° linear correlation). Consequently the uncertainties are relatively small for (the logK1,IV° estimate of) MS2O32+. Conversely the biggest uncertainty is for MCO32+ the most stable complex, since CO32 appeared to be the most reactive RO22 ligand we studied: it is one of the endpoints of the correlation lines. For this reason, we estimated an upper bound for the value of logK1,IV° from experimental data: we re-interpreted available published experimental solubilities of actinides(IV) by using the same methodology as in Refs. [4,5,29]. We obtained formation constants consistent with the original interpretation (by the authors of Ref. [9]). However, for consistency we added the known stabilities of An(CO3)i4−2i (i = 4 and 5), and we tested many possible complexes for which we also estimated maximum possible stabilities (Table 2) for sensibility analysis; our purpose was essentially to estimate a value for AnCO32+. As expected from pH vs log[CO32] predominance diagrams [7,26,28], the most restrictive conditions were found for published solubilities of An(IV) measured at low pH and high CO2 partial pressure, namely in Refs. [9,30]. In both cases we obtained virtually the same values: logK1,IV°11.1or11.5 for the ThCO32+ complex. Using this value (and the sulfate data), logK1,IV°=(5.8+0.51pKa°) is calculated, from which logK1,IV°=11.1, 9.5, 6.8 and 6.6 are calculated for the An(IV) complexes CO32, SO32, SO42 and S2O32, respectively. These values are within the uncertainties of the previous estimates above. The AnCO32+ value of 11.1 is identical as, or slightly smaller (by 0.9 ± 1.9) than the central value (12.0 ± 1.9) of the above estimate, and the same trend was observed for AmCO3+: it was overestimated by 0.8 log10 unit when using the same slope (0.62) for all oxidation states. This correlation predicts under-stabilisation of the AnCO3X2 complexes (X = 3 or 4) as compared to the AnO2CO3X6 systems (X = 5 or 6): it might be attributed to the planar structure of CO32, which offers a better fit to the geometry of the coordinating equatorial plane of the AnO2X4 actinyl cations. Unfortunately, such an explanation would also hold for the SO32 (quasi-) planar ligand (Fig. 5), which is not specially confirmed (Fig. 4). For consistency with our estimate, logK1,IV°11.1, we prefer the last correlation logK1,IV°=(5.8+0.51pKa°), where we increased the uncertainties to encompass the previous logK1,IV°=((5.6±0.2)+(0.62±0.16)pKa°) correlation: logK1,IV°=((5.8±0.4)+(0.51±0.27)pKa°) from which we obtain logK1,IV°11.1±3.2 for AnCO32+,=9.5±2.3 for AnSO32+, 9.5 ± 2.3 for AnHPO42+ and 6.6 ± 0.8 for AnS2O32+, where uncertainties are increased take into account the lack of experimental data.

Table 2

Stabilities of An(IV) carbonato complexes

An(CO3)i(OH)j42ijlogKi,j°(Pu)a, [5]logKi,j°(Th) a, [9]logKi,j°(Th) a,b
AnOH3+13.2
AnCO32+≤11.1
AnCO3OH+<21.1
An(CO3)2<20.8
AnCO3(OH)2<<4227.0<30.1
An(OH)4<<47.938.5
An(CO3)2OH<<40.5<29.4
AnCO3(OH)3<<47.734.8≤38.5
An(CO3)32<<37.6<27.5
An(CO3)2(OH)22<<46.233.3≤36.8
AnCO3(OH)42<<51.837.4≤39.9
An(CO3)3OH3<<42≤34
An(CO3)2(OH)33<<50.5<39.2
An(CO3)443729.9
An(CO3)3(OH)24<41<37.6
An(CO3)3(OH)35<40.5<37.9
An(CO3)4OH5<3934.435.4
An(CO3)5635.628.4
An(CO3)4(OH)26<37<36.4
An(CO3)3(OH)46<38.5<39.3
An(CO3)68<36

a Ki,j° is the standard constant of equilibrium An4++iCO32+jHOAn(CO3)i(OH)j42ij.

b Maximum possible values from the experimental data of Ref. [9]; these estimations are consistent with the original interpretation (by the authors of Ref. [9]), but the known stabilities of An(CO3)i42i are here added in the fits for i = 4 and 5; our purpose was essentially to estimate a value for AnCO32+ (see text) and to outline sensitivity analysis.

Fig. 5

Geometry of RO22 ligands optimized in the gas phase. R–O bond distances (Å), angles (°) and zO, the atomic charge of O are written on the figure.

logK1,IV°=6.6±0.8, the value estimated for AnS2O32+ is quite similar to those for the AnSO42+ complexes: 6.58 (USO42+), 6.85 (NpSO42+) and 6.89 (PuSO42+). This is consistent with the molecular structures of these ligands: SO42 and S2O32 both have a tetrahedral structure (Fig. 5) – an O atom (of SO42) being replaced by an S atom (in S2O32) and similar pKa values. SO32 has a different structure, a higher pKa, and a higher estimated value of the complexation constant. Unfortunately, it does not seem easy to draw any simple explanation just from zO, the atomic charge of the free ligand in vacuum (Fig. 4(d)).

3.4 Uranium geochemistry

In non-complexing aqueous solutions, the solubility of uranium is controlled by uraninite (UO2(s)) in reducing conditions, and by schoepite (UO3(s)) in oxidizing and neutral conditions, as illustrated in Fig. 6(a) for 1 μM [U]t. Besides those minerals, aqueous U(VI) species are predominant in a large E–pH domain, while aqueous U(IV) is stable only in reducing conditions, and is mostly hydrolysed.

Fig. 6

Pourbaix' diagrams of uranium. The predominance domains of the major soluble species are shown as a function of pH and E, the redox potential of the solution in (a) non-complexing media, (b) adding the influences of carbonate ions and (c) sulfur species for a typical ([CO3]t = 11.6 mM and [S]t = 4 mM) composition of underground waters [11]. pH is taken into account for carbonate speciation (CO32/HCO3/CO2(aq)), but its reduction -typically into CH4(g)- is not. The mixed dashed line (A) represents redox conditions in Fig. 7. The long dashed line (B) represents redox conditions for Fig. 1(c).

On adding carbonate ions at a typical concentration of underground waters, aqueous U(VI) carbonate complexes prevail between pH 4 and 12 (Fig. 6(b)), UO3(s) is totally dissolved, and the UO2(s) stability domain is reduced. We have ignored the reduction of the carbonate ions, since this reaction is usually very slow. Nevertheless, no carbonate complex of U(IV) appears on the Pourbaix diagrams. These complexes would predominate only at higher carbonate concentrations than those studied here. Rai and Ryan have already proposed that carbonate complexation of actinide(IV) ions in environmental conditions can be neglected [8].

Among the main sulfur species, SO42 prevails over a large (E–pH) domain in oxidizing to slightly reducing equilibrium conditions. For high pH values, this domain even extends to reducing conditions (Fig. 1). H2S, HS and S2− prevail in reducing conditions, but their complexing properties are not significant for the (hard) cations of the f-block elements. UO2SO4(aq) is the only predominating sulfur species in our typical underground water conditions (Fig. 6(c)), even when introducing the 1:1 complexing constants estimated above (Section 3.3, Table 1) for US2O32+ and USO32+. The RO22 ligands are protonated in acidic conditions, which decreases their concentrations: the corresponding 1:1 complexing equilibria are actually the

MzM+HRO2MRO2zM2+H+(4a)
exchange equilibria. Despite this concurrency between protonation and complexation, several complexes with the 1:1 stoichiometry are stable. Conversely, the existence of AnCO32+ complexes has never been demonstrated, since the competition (between carbonate complexation and hydrolysis) is less favourable for An(IV) as compared to the actinides in the other oxidation states (Fig. 4(a)). The RO22 concentration can also dramatically decrease as a result of redox reactions: CO32 and SO42(=RO22) are reduced in equilibrium chemical conditions, where U(IV) is stable. However, these reduction reactions are very slow: we have neglected them for the carbonate ions (Fig. 6). Similarly, the coexistence of S(VI) and S(–II) species is often observed, indicating that equilibrium conditions are not always achieved for S natural aqueous systems. The achievement of the U(IV)/U(VI) equilibrium can also last for a few weeks. Note that all these slow reduction reactions can be explained by the need to break strong (covalent) bonds of O(–II) with C, S(VI) or U(VI). Such kinetics is often handled by assuming that slow reactions are blocked, while equilibrium conditions are achieved for other equilibria. Therefore U(IV) might be complexed by sulfate and by reduced intermediary S species as well: we plotted two simulation diagrams corresponding to different kinetic assumptions that would stabilise sulfoxy-anion complexes (Fig. 7). In both cases, E, the redox potential, is at the limit of the HiS2O3i2 and HiSO3i2 domains as plotted in Fig. 1(a) and line A in Fig. 6(c). In such redox conditions uranium aqueous species are stable at the +6 oxidation state. For this reason, the sulfate complexes of U(IV) are not seen on the Pourbaix' diagram (Fig. 6(c)): the sulfate ions are reduced in equilibrium conditions, where U(IV) aqueous species are dominating – i.e. the S(VI)/S(–II) frontier (line B in Fig. 6(c)) is at higher potentials than the U(VI)/U(IV) frontier in most pH conditions. Nevertheless we only focus on U(IV) aqueous species; therefore they would be minor uranium species, when the U(VI)/U(IV) equilibria are achieved.

Fig. 7

Distribution of aqueous U(IV) species. lg(|A|/|U(IV)|) is plotted as a function of pH, where |A| is the activity of an aqueous U(IV) species. The conditions are the same as in Fig. 6(c) except for E, the redox potential of the solution: HiSO4i2 (a) or HiS2O3i2 (b) are assumed to be the major S species. [HiS2O3i2]/[HiSO3i2]10 (mixed dashed line A in Fig. 6(c)). [HiSO4i2]/[HiS2O3i2]20(a) or1(b) corresponding to different kinetic assumptions.

In a first simulation we assumed that the major S species are in the +6 oxidation state (namely HiSO4i2), and that [HiSO4i2]/[HiS2O3i2]20, which corresponds to much higher HiS2O3i2 concentrations than in equilibrium conditions (for kinetic reasons). The USO42+ and U(SO4)2 sulfate complexes still appear to be U(IV) major species in acidic conditions (pH < 3.2), while nearly 1% of US2O32+ can be formed (Fig. 7(a)). Yet US2O32+ might be a kinetic intermediate, since uraninite is often associated with pyrite (FeS2); indeed S2O32 is an intermediary product of its oxidative dissolution [14]. S2O32 and SO42 have similar reactivities (Fig. 4), but (in our hypothesis) the SO42 concentration is higher (than the S2O32 one), which explains why the sulfate complexes dominate. Although the SO32 ions are more reactive, the USO32 complex is negligible in that simulation, since the reduction of SO32 (Fig. 1) decreases its concentration and consequently its complexing ability.

In a second simulation we assumed that SO42 is not formed at all (still for kinetic reasons). Therefore S2O32 is now the dominating S aqueous species (Fig. 1). US2O32+ is the major U(IV) species, and up to a few percent of USO32+ is formed at pH less than 3.0 (Fig. 7(b)).

Even if sulfoxy-anions complexes of U(IV) are certainly not stable in equilibrium conditions, these simulations indicate that they might be formed as kinetic intermediates typically in the course of uraninite and pyrite oxidative dissolutions, or the interaction of aqueous uranium (including U(VI)) with pyrite surfaces [31,32]. Of course, this statement needs experimental confirmation, specially in less acidic to basic conditions, where higher and mixed An(RO2)i(OH)j4-2i-j complexes may form. Furthermore, at pH > 3, U(OH)4(aq) is the dominant species in equilibrium reducing conditions, and the formation of S-containing complexes is less favourable, specially in the usual pH conditions of equilibrated natural under-ground-waters.

3.5 Hydrolysis of actinide ions

We now consider where protactinium should be placed in the above correlations. There is little experimental published information on Pa aqueous chemistry. We first used the available information on its aqueous species in non-complexing (acidic) solutions and its hydrolysis. Although Pa is known to be an f-block element [17], Pa(V) aqueous chemistry is closer to that of d elements (in the same – the 5th- column of the periodic table) than to the chemistry of AnO2+ trans-protactinian actinide cations. Here is a strong indication that PaO2+ is not the dominating aqueous Pa species.

Jaussaud et al. recently reviewed previous measurements of the Pa(V) hydrolysis constants [33–35]. They are essentially based on liquid–liquid extraction measurements, for which they pointed out many experimental difficulties: side reactions of the organic chemicals, sorption of Pa(V) on vessels and unwanted side reactions. Surprisingly most of the published studies are only based on measurements at trace concentrations of Pa(V), while sorption reactions are classically made negligible by using (non-radioactive) chemical analogues at macro-concentrations. The hydrolysis behaviour of Nb(V) seems similar to that proposed for Pa(V): we calculated the hydrolysis constants of Nb(V) by interpreting the aqueous solubility of Nb2O5 (Table 3) [25]. We re-interpreted the (Pa) experimental measurements at 25 °C by giving more weight to the data for which systematic errors seemed the lowest according to Jaussaud's observations and comments: the measurements at high ionic strength and low (TTA = thenoyltrifluoroacetone) extractant concentrations. This essentially confirmed the original interpretation given by the authors, but increased uncertainties and suppressed an inconsistency (Table 3). Namely, classical slope analyses of the raw experimental data imply a Pa species of charge +1 in 0.1 M NaClO4 aqueous solutions, and +2 at higher I (3 M), consistent with the PaOOH(OH)+ and PaOOH2+ stoichiometries [17] in nearly the whole pH range studied (0 < pH < 4). Indeed, high I usually stabilizes the species of higher charges. There is no clear evidence for neutral and tri-cationic species, due to the scattering of the data (Fig. 8). A neutral Pa species is certainly formed, but it is clearly stabilized by increasing [TTA]t, the total TTA concentration (results not reproduced in Fig. 8 for clarity), which suggests the corresponding aqueous Pa(V) species might include ionised deprotonated TTA ligands, which would hide the formation of PaOOH(OH)2 (equivalently written PaO(OH)3 or Pa(OH)5). The formation of PaO3+ is not proposed in the original interpretation of the authors, and indeed needs confirmation. Finally, the most reliable published Pa(V) standard hydrolysis constant is certainly logK2°=(1.26±0.15) [34], a value consistent with our re-interpretation (−1.44 ± 0.71) (Fig. 8 and Table 3), and with the Nb(V) value of −(1.65 ± 0.2). We also propose logK1°=(0.04±0.36) and logK3°3.6, giving no credit to the (−(7.15 ± 0.4) [34] and −(7.03 ± 0.15) [33,35]) published interpretations from the same experimental data, because the corresponding hydrolysis species could not be detected in the experimental conditions used in this work (pH ≤ 4 corresponding to a maximum value of about 10−7+4 = 0.1% for the concentration of the hydrolysis species to be compared with the 100.71 uncertainty in the measurements of logK2°), even though they might provide a reasonable value: the Nb(V) value of −4.95 ± 0.2 might also be used as a rough estimate. These Ki° constants are actually the hydrolysis constant of PaO3+, but the existence of this species has not clearly been demonstrated, and for this reason only the second stepwise experimental determination (K2°) seems reliable, while the first one (K1°) needs confirmation. The logK2° value is only a little larger than that of logK1°: the two first hydrolyses of Pa(V) would be nearly simultaneous. This is also observed for the other actinide aqueous species.

Table 3

Pa(V) hydrolysis constants

I (M)logK1,m*elogK2,m*elogK3,m*e
3−2.0 ± 0.15−5.8 ± 0.3b[33–35]
3−0.35 ± 0.29−1.75 ± 0.91<−3.5[33]a
1−1.7 ± 0.2−6.9 ± 0.6b[33–35]
1≥−0.38c−1.50 ± 0.56<−0.98b[33]a
0.5−1.6 ± 0.2−6.9 ± 0.6b[33,35]
0.5−1.49c[33]a
0.1−1.5 ± 0.2−7.0 ± 0.6b[33–35]
0.1−1.80 ± 0.34[33]a
0.1≥−0.75c−1.65 ± 0.2−4.95 ± 0.2Nb(V)d [36]a
0−1.24 ± 0.02−7.03 ± 0.15b[33,35]
0−1.26 ± 0.15−7.15 ± 0.4b[34]
0−0.04 ± 0.36c−1.44 ± 0.71<−3.6c[33]a

a Graphically interpreted in the present work, and (last line) extrapolated to I = 0. See Fig. 8.

b Inconsistent value: since the measurements were performed at pH < 4, it is not possible to fit (from them) a log Ki value less than about −4.

c There is no clear experimental evidence of the corresponding reaction.

d Nb(V) value fitted from solubility data [36], and here tabulated for comparison.

e Ki,m is the constant of equilibrium PaO.

Fig. 8

Predominance diagram of Pa(V) hydrolysis species. The hydrolysis constants (logK1°=pH1/2,i correspond to the pH value at the half-point reaction) have been re-interpreted (Table 3) from the experimental data (and extrapolated to I = 0 with the SIT formula) of Ref. [30] in 0.1 (open symbols) to 3 (filled symbols of the same colour) mol L−1 NaClO4 aqueous solutions for the successive formations of species of charges 2 (dark blue diamonds), 1 (brown circles) and 0 (green squares). The intermediate colours are for 1 and 0.5 mol L−1. For higher [TTA] values (not represented), a systematic deviation is observed. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)

A linear regression correlates reasonably well logK1°, the logarithm of the hydrolysis constant and z, the charge of the Mz+ cation for Mz+ = Ra2+, Am3+, Pu4+ (brown dotted line in Fig. 9(a)): logKi°(25.246.16z). Adding other An3+ and An4+ cations gives virtually the same correlation (brown dotted line in Fig. 9(b)): logKi°(25.676.28z). The actinyl ions such as UO22+, and NpO2+ are below this correlation line. This stabilisation of their hydrolysed species might be originated in an An–Oyl (intramolecular) charge transfer induced on approaching equatorial water molecules to UO22+ [21], an inductive effects of Owater, which is expected to be more important for HO equatorial ligands (where Oyl is an O atom of an AnO2(X4)+ actinyl cation, while Owater is an O atom of a water equatorial ligand of the actinyl). Nevertheless, the actinyl cations are not far from the correlation lines: including UO22+ and NpO2+ but excluding Ra2+ (black lines in Fig. 9); logKi°(24.436.01z), and adding other AnO22+ and AnO2+: logKi°(22.735.62z), while extrapolating (blue dashed lines in Fig. 9) the AnO22+ and AnO2+ data give (blue dashed) lines (in Fig. 9) that are between the PaO3+ and PaOOH2+ data. Namely, PaO3+(aq) would be intermediate between the bare actinide hard cations (Anx+, X = 3 or 4), and the usual actinyl cations (AnO2(x−4)+, X = 5 or 6, An = U, Np, Pu or Am), whose polarization by equatorial ligands essentially results in the intramolecular charge transfer between An and Oyl. Note that despite this polarization AnO2(X4)+ is still a hard cation for the equatorial ligands – i.e. negligible charge transfer from the equatorial ligands – we may wonder whether assuming this concept of (equatorial) hardness/(intra-actinyl) softness is still relevant for such a behaviour. However, the AnO2+/AnO22+/PaO3+ set can as well be considered as a series of actinide oxo-cations on their own line (approximately the blue dashed one in Fig. 9). As a result of their polarizability, these oxo-cations are below the (brown dotted) line of the bare Mz+ cations. This polarizibility probably decreases zAn, the An atomic charge (in the oxo-cation) by inductive effect. This inductive effect would increase with zAn: this can explain the quite surprisingly low calculated atomic charges of U (in UO2+) and Pa (in PaO3+), which would therefore be at the origin of the more negative slope for the actinide oxo-cation line (as compared to that of the bare cations). However, this proposition needs confirmation, since the Pa(V) experimental results are neither very accurate, nor validated, and since physical explanations cannot be proven by only such empirical correlations: it is a limit of such correlations, rather than problems in the (very few) experimental data. According to this interpretation, two correlation lines are expected, one for the bare Anx+ cations, and the other one for the AnO2(X4)+ di-oxo-cations, while the PaO3+ mono-oxo-cation would fall between these lines, as we actually observed.

Fig. 9

Hydrolysis constants of actinide aqueous ions: logK1° (filled symbols) and logK2° (open symbols). The lines represent the linear regressions for (black line) all the points plotted on the figures excepted Pa(V) ones, only for the (Ra2+, Am3+, Pu4+) set (brown dotted line), or only for the (NpO2+ and UO22+) set (blue dashed line). The (atomic NPA) charges of the bare An cations (written in Fig. 9(a)) were obtained from DFT quantum calculations in gas phase. Data and corresponding lines are plotted by adding similar cations in Fig. 9(b). Two possible stoichiometries of “the” Pa(V) aqueous species are represented. (For interpretation of the references to colour in this figure legend, the reader is referred to the web version of this article.)

The experimental hydrolysis data already indicated that PaOOH2+ should be only a little less reactive than Pu4+. Finally, it can still be under debate whether it is better to consider PaO3+ or PaOOH2+, or both, as the Pa(V) aqueous species. Using atomic charges here gave unexpected correlations, as compared to the other chemical studies reviewed in the present paper. Nevertheless, it is not a convincing proof that PaO3+ can be an important aqueous species (i.e. it is not a confirmation of the above estimate of its logK1° value).

Our DFT calculations give pictures confirming that PaOOH2+(aq) is a logical species, similar to UO22+(aq), while PaO2+(aq) is easily protonated. PaOOH2+ is merely protonated PaO2+ [17]. However, our continuing DFT calculations for interpreting published EXAFS measurements clearly indicate that the PaO3+ geometry exists in the Pa(V) species formed in concentrated H2SO4 aqueous solutions, since none of the other tested model geometries reproduced the shortest experimental distance, including PaOOH2+: it thus dissociates in complexing media, to give aqueous sulfate complexes of PaO3+.

4 Concluding remarks on analogies

(logK1°vspKa°) correlations fit experimental data surprisingly well (Fig. 4(a)) for An(III to VI) complexes, with the exception of Pa(V), with 5 anionic potentially bidentate and oxygen-donor RO22 ligands: CO32>HPO42SO32>SO42S2O32 in the order of their reactivities for the actinide cations. Furthermore, the correlations were found to be linear for each oxidation state, and the (0.62) slope of these H+/An(X) correlations are approximately the same for all the oxidation states (X). This, of course, means that the An(X)/An(VI) correlations deduced are linear with slope 1: ax, the intercept, is interpreted in terms of the half-reaction point, namely 10aXaVI=K1,X°/K1,VI°=KX/VI° (Eq. (5)) is the constant of the exchange equilibrium (Eq. (6a)), which takes (approximately) the same value (Table 1) for all the RO22 ligands and for a given oxidation state. This provides ax-based numerical values for the (up to now qualitative) scale An4+>AnO22+An3+>AnO2+ (Fig. 4(c)). The ax values are also qualitatively correlated with z, the atomic charges of the An cations; but using z instead of ax gives a poorer correlation (Fig. 4(c)). Note that measured energies of reactions, namely the hydrolysis constants of Pa(V), show clearly that it is not an analogue of the other An(V); quantum calculations give a chemical explanation of the destabilization of PaO2+ by hydration, typically resulting in clear apical H-bonded water molecules [17]. The two approaches are complementary: the experimental energies of reactions are more accurate, while quantum calculations provide geometries and other physical data that can be interpreted in terms of usual chemical concepts (atomic charges, bounds and their covalency…) for explaining the observed chemical reactivity.

All the correlations we used are totally empirical. They are probably the result of various effects, some of them more or less cancel out in the correlations, and they are probably all linked to z, the atomic charge of the cation. However, using z gave poorer linear correlations than comparing only experimental equilibrium constants. For this reason it is certainly better to interpret the correlations with chemical concepts than with any unproven physical explanation. Nevertheless, intramolecular charge transfers are deduced from quantum calculations for (possibly protonated) actinyl cations. These charge transfers can be related to the slightly different trends observed between the hydrolysis behaviours for these three types of cations (Anz+, PaO3+ and AnO2X4). This effect would be a little less important for PaO3+ than for the AnO2X4 cations where there are more covalent bounds (which promote intramolecular charge transfer). Charge transfers might as well be at the origin of the over-stabilisation of the aqueous An(IV) hydroxides as compared to complexation. These predictions are restricted to similar complexes: the linear regression numbers cannot be extended to all types of complexes without validations. Indeed such approaches exist in literature, overestimating the stabilities of several hypothetical chemical species, whose existences have never been confirmed; such numbers are not considered in the NEA-TDB reviews [2–4].

We have in fact used several levels of analogy giving different rules of thumb (although some have similar mathematical equations), which might indicate that it is hopeless to develop more general empirical formulae. The strongest analogy is typically for M(X), the series of An and other cations with the same charge, geometry and oxidation state, X: they form soluble complexes and hydroxides with the same stoichiometries, and with stability constant values that hardly differ by more than the experimental uncertainties. For different (X) oxidation states, linear correlations were found with slope 1, which defines a second type of analogy: considering the AnO22+/MzM exchange, a solution of activity |AnO22+| has the same reactivity (for the LzL ligand) as a solution of activity of |MzM| KX/VI°.

In a third type of analogy, the (b) slope of the linear correlation is no longer 1 (b ≠ 1), as found here for H+/MzM exchanges (b is not correlated with zM, so it is better not to consider the (zM H+)/MzM exchange); K1°Ka°, as its equilibrium constant is no longer a relevant parameter for the analogy: analogue solutions should instead be based on the parameter 10a=K1°(Ka°)b=[(|H+||LzL|)b/(|MzM||LzL|)]1/2, whose interpretation is less intuitive.

Similarly when comparing protonation energies in aqueous and gas phases, the slope of the energy correlation is determined by the ratio of the H+ activity scales in both phases. It is far from 1 (Fig. 3), corresponding to differences in the orders of magnitude (of the energies of reactions) in each phase. The intercept is related to the balance of the hydration energies (it can be shifted by changing the reference state).

Finally the rules of thumb are often characterized by the value of the (fitted) slopes in linear correlations of equilibrium constants, or equivalently by the value of the corresponding exponent in ratios of activities, the ideal concentrations.


Bibliographie

[1] V. Phrommavanh, M. Descostes, P. Vitorge, C. Beaucaire, J.-P. Gaudet, Estimating the Stabilities of Aqueous Actinide Complexes with Sulfoxyanions. Poster PA3-9 at the 10th International Conference on Chemistry and Migration Behaviour of Actinides and Fission Products in the Geosphere, MIGRATION'05, September 18–23, 2005. Avignon, France.

[2] I. Grenthe; J. Fuger; R.J.M. Konings; R. Lemire; A.B. Muller; C. Nguyen-Trung; H. Wanner Chemical Thermodynamics of Uranium, Paris OECD/NEA and Elsevier, 1992 http://www.nea.fr/html/dbtdb/pubs/uranium.pdf (download:)

[3] R. Silva; G. Bidoglio; M. Rand; P. Robouch; H. Wanner; I. Puigdomenech Chemical Thermodynamics of Americium, Paris OECD/NEA and Elsevier, 1995 http://www.nea.fr/html/dbtdb/pubs/americium.pdf (download:)

[4] R. Lemire; J. Fuger; H. Nitsche; M. Rand; K. Spahiu; J. Sullivan; W. Ullman; P. Vitorge Chemical Thermodynamics of Neptunium and Plutonium, Paris OECD/NEA and Elsevier, 2001

[5] P. Vitorge; H. Capdevila Radiochim. Acta, 91 (2003), p. 623

[6] P. Vitorge; H. Capdevila; S. Maillard; M.-H. Fauré; T. Vercouter J. Nuclear Sci. Technol., 3 (2002), p. 713

[7] P. Vitorge, Chimie des actinides. Article and Formulaire form B, Techniques de l'Ingénieur, vol. 3520, 1999.

[8] D. Rai; J.L. Ryan Inorg. Chem., 24 (1985), p. 247

[9] M. Altmaier, V. Neck, T. Fanghänel, Solubility and Ternary Hydroxo-Carbonate Complexes of Thorium, NRC6 2004-03, Sixth International Conference on Nuclear and Radiochemistry, August 29th to September 3rd 2004, Aachen, Germany. http://www.fz-juelich.de/NRC6/.

[10] M. Altmaier; V. Neck; R. Muller; T. Fanghanel Radiochim. Acta, 93 (2005) no. 2, p. 83

[11] E. Gaucher; C. Robelin; J.-M. Matray; G. Negral; Y. Gros; J.-F. Heitz; A. Vinsot; H. Rebours; A. Cassagnabere; A. Bouchet Phys. Chem. Earth, 29 (2004) no. 1, p. 55

[12] C. Beaucaire; P. Toulhoat Appl. Geochem., 2 (1987), p. 417

[13] M. Descostes; C. Beaucaire; F. Mercier; S. Savoye; J. Sow; P. Zuddas Bull. Soc. Géol. Fr., 173 (2002), p. 265

[14] M. Descostes; P. Vitorge; C. Beaucaire Geochim. Cosmochim. Acta, 68 (2004) no. 22, p. 4559

[15] P.L. Brown; E. Curti; B. Grambow Chemical Thermodynamics of Zirconium, Paris OECD/NEA and Elsevier, 2005

[16] R. Guillaumont; T. Fanghänel; J. Fuger; I. Grenthe; V. Neck; D.A. Palmer; M.H. Rand Update on the Chemical Thermodynamics of Uranium, Neptunium, Plutonium, Americium and Technetium, Paris OECD/NEA and Elsevier, 2003

[17] B. Siboulet, C.J. Marsden, P. Vitorge, Submitted for publication.

[18] R.L. Martin; P.J. Hay; L.R. Pratt J. Phys. Chem. A, 102 (1998) no. 20, p. 3565

[19] G. Choppin Principles and Practices of Solvent Extraction (J. Rydberg; C. Musikas; G. Choppin, eds.), Marcel Dekker, New York, USA, 1992, p. 71

[20] B. Allard; J. Rydberg; C. Musikas; G. Choppin Principles and Practices of Solvent Extraction (J. Rydberg; C. Musikas; G.R. Choppin, eds.), Marcel Dekker, New York, USA, 1992, p. 209

[21] B. Siboulet, C.J. Marsden, P. Vitorge, Chem. Phys., in press.

[22] A.E. Reed; L.A. Curtiss; F. Weinhold; W. Aas Chem. Rev., 88 (1988), p. 899

[23] E.D. Glendening, A.E. Reed, J.E. Carpenter, F. Weinhold, Natural Population Analysis in Gaussian. NBO Version 3.1, 1998.

[24] M.J. Frisch; G.W. Trucks; H.B. Schlegel; G.E. Scuseria; M.A. Robb; J.R. Cheeseman; V.G. Zakrzewski; J.A. Montgomery; R.E. Stratmann; J.C. Burant; S. Dapprich; J.M. Millam; A.D. Daniels; K.N. Kudin; M.C. Strain; O. Farkas; J. Tomasi; V. Barone; M. Cossi; R. Cammi; B. Mennucci; C. Pomelli; C. Adamo; S. Clifford; J. Ochterski; G.A. Petersson; P.Y. Ayala; Q. Cui; K. Morokuma; D.K. Malick; A.D. Rabuck; K. Raghavachari; J.B. Foresman; J. Cioslowski; J.V. Ortiz; A.G. Baboul; B.B. Stefanov; G. Liu; A. Liashenko; P. Piskorz; I. Komaromi; R. Gomperts; R.L. Martin; D.J. Fox; T. Keith; M.A. Al-Laham; C.Y. Peng; A. Nanayakkara; C. Gonzalez; M. Challacombe; P.M.W. Gill; B.G. Johnson; W. Chen; M.W. Wong; J.L. Andres; M. Head-Gordon; E.S. Replogle; J.A. Pople Gaussian 98 (Revision A.9), Gaussian, Inc., Pittsburgh PA, 1998

[25] M.J. Frisch; G.W. Trucks; H.B. Schlegel; G.E. Scuseria; M.A. Robb; J.R. Cheeseman; J.A. Montgomery; T. Vreven; K.N. Kudin; J.C. Burant; J.M. Millam; S.S. Iyengar; J. Tomasi; V. Barone; B. Mennucci; M. Cossi; G. Scalmani; N. Rega; G.A. Petersson; H. Nakatsuji; M. Hada; M. Ehara; K. Toyota; R. Fukuda; J. Hasegawa; M. Ishida; T. Nakajima; Y. Honda; O. Kitao; H. Nakai; M. Klene; X. Li; J.E. Knox; H.P. Hratchian; J.B. Cross; V. Bakken; C. Adamo; J. Jaramillo; R. Gomperts; R.E. Stratmann; O. Yazyev; A.J. Austin; R. Cammi; C. Pomelli; J.W. Ochterski; P.Y. Ayala; K. Morokuma; G.A. Voth; P. Salvador; J.J. Dannenberg; V.G. Zakrzewski; S. Dapprich; A.D. Daniels; M.C. Strain; O. Farkas; D.K. Malick; A.D. Rabuck; K. Raghavachari; J.B. Foresman; J.V. Ortiz; Q. Cui; A.G. Baboul; S. Clifford; J. Cioslowski; B.B. Stefanov; G. Liu; A. Liashenko; P. Piskorz; I. Komaromi; R.L. Martin; D.J. Fox; T. Keith; M.A. Al-Laham; C.Y. Peng; A. Nanayakkara; M. Challacombe; P.M.W. Gill; B. Johnson; W. Chen; M.W. Wong; C. Gonzalez; J.A. Pople Gaussian 03, Revision C.02, Gaussian, Inc., Wallingford CT, 2004

[26] P. Robouch, Contribution à la prévision du comportement de l'américium, du plutonium et du neptunium dans la géosphère ; données géochimiques, Thèse n°1987STR13230-SUDOC n° 043589863, Université Louis-Pasteur, Strasbourg, France, 13/11/1987. CEA-R-5473, 1989.

[27] P. Robouch; P. Vitorge Inorg. Chim. Acta, 140 (1987) no. 1/2, p. 239

[28] T. Vercouter, Complexes aqueux de lanthanides (III) et actinides (III) avec les ions carbonate et sulfate. Étude thermodynamique par spectrofluorimétrie laser résolue en temps et spectrométrie de masse à ionisation électrospray. Thèse n°2005EVRY0003-SUDOC n° 09483699X, Université d'Évry, France, 2005.

[29] P. Vitorge, H. Capdevila, CEA-R-5793, 1998.

[30] E. Osthols; J. Bruno; I. Grenthe Geochim. Cosmochim. Acta, 58 (1994) no. 2, p. 613

[31] N. Eglizaud, M. Descostes, M. Schlegel, E. Simoni, 10es Journées nationales de radiochimie et de chimie nucléaire, Avignon, France, 7–8/09/2006.

[32] N. Eglizaud, F. Miserque, E. Simoni, M. Schlegel, M. Descostes, MIGRATION'05, Avignon, France, 18–23/09/2005.

[33] C. Jaussaud, Contribution à l'étude thermodynamique de l'hydrolyse de Pa(V) à l'échelle des traces par la technique d'extraction liquide–liquide avec la TTA, Thèse, Université Paris-Sud, Orsay, 2003.

[34] D. Trubert; C. Le Naour; C. Jaussaud J. Solution Chem., 31 (2002) no. 4, p. 261

[35] D. Trubert; C. Le Naour; C. Jaussaud; O. Mrad J. Solution Chem., 32 (2003) no. 6, p. 505

[36] C. Peiffert, Solubilité et hydrolyse du niobium en solution aqueuse à 25 °C et 0.1 Mpa, Poster, Journées scientifiques ANDRA, Nancy (France), 7–9 décembre 1999.


Partially presented at the Migration 05 conference [1], and part of the Ph.D. thesis of V. Phrommavanh.

Commentaires - Politique