1 Introduction
Quantitative information on chemical equilibria and the thermodynamics of actinide complex formation in aqueous solution is, in addition to its pure scientific interest, also important in many applications, such as separation technology and speciation in ground and surface water systems. Solution thermodynamic data provide information on the stoichiometry of the complexes formed, but little or no information on the structure of the first coordination sphere, i.e. the mode of coordination of ligands and the possible formation of isomers of various types. There is extensive information [1] on the chemical thermodynamics of actinide complex formation reactions, including the following uranyl(VI) sulfate reactions that will be discussed in this communication:
UO22+ + SO42− ⇌ UO2(SO4) | (1) |
UO22+ + 2SO42− ⇌ UO2(SO4)22− | (2) |
Single-crystal X-ray structures of a number of uranyl(VI) sulfate compounds reveal that the sulfate ion can be coordinated both in bidentate [2] and monodentate [3] fashion, but with predominance of the latter type. A characteristic feature of the bond distances in solid compounds that contain both mono- and bidentate sulfate coordination is that the U–Osulfate bond distance is about 0.15 Å shorter in the former than in the latter. Structure studies of uranyl(VI) sulfate complexes in solution have been made by Moll et al. using EXAFS [4] and by Neuefeind et al. using large-angle X-ray diffraction (LAXS) [5]. Both methods provide accurate information on bond distances, but much less precise information on coordination numbers. These two studies demonstrate that the sulfate ion also in solution can be coordinated both in a monodentate [5] and bidentate fashion [4]. The solution structure determinations are complicated by the fact that the test solutions may contain a mixture of different uranyl species, as will be discussed in the following.
The equilibrium constants for the formation of uranyl(VI) sulfate complexes are moderately large at zero ionic strength, log β01 = 3.15 and log β02 = 4.14, respectively [1] for the reactions (1) and (2). These equilibrium constants decrease substantially with increasing ionic strength [1], cf. Appendix, and from their magnitude it is not possible to make a reliable prediction of the mode of coordination of sulfate and the number of donor atoms in the equatorial plane of the UO22+ ion. We have therefore used quantum chemical methods to determine the structure and bond distances in various isomers of UO2(SO4) and UO2(SO4)22− and their relative electronic and Gibbs energies of formation. The bond distances can be directly compared with experimental data, while the relative energies are used to corroborate the conclusions drawn from the bond distances.
The quantum chemical calculations refer to a model where the complexes are dissolved in a model solvent where the water is described by a conductor-like polarizable continuum model (CPCM) with the macroscopic dielectric constant of water [6]. This differs in two ways from the experimental situation: hydrogen bonding is treated in a simple way and interactions between the dissolved complexes and the counter ions, Na+ in most experimental test solutions, are not taken into account, i.e. the calculation refers to a medium with zero ionic strength. For this reason we have tested how sensitive the geometry and the relative energy of the complexes are for different model assumptions.
The M–O–S bond angle in most solid-state structures of metal sulfates differs significantly from 180°, the configuration that is expected to minimize electrostatic repulsion. To rationalize this observation, we have probed the interaction between the uranyl unit and a unidentate sulfate group through the energy decomposition scheme [7] in ADF [8,9] (see Section 2.1) to determine whether the angular preference is driven by electronic or steric effects.
2 Quantum chemical calculations and results
2.1 Computational details
The small-core relativistic effective core potential of the Stuttgart–Bonn group was employed for uranium [10], along with the segmented (14s13p10d8f6g)/[10s9p5d4f3g] basis set [11]. The valence triple-ζ plus polarization basis sets were used for sulfur [12a], oxygen [12b], and hydrogen [12b] atoms. The correlation-consistent double-ζ plus polarization basis set (cc-pVDZ) from Woon et al. [13] was used for the sodium counterion. The structures of the various isomers of UO2(SO4) and UO2(SO4)22− species were optimized with the gradient-corrected hybrid functional B3PW91 [14], as this is known to result in geometries that are close to the experimental values. Single-point MP2 calculations were also performed at the DFT-optimized geometries, with the 1s core orbitals of oxygen, the 1s, 2s, and 2p orbitals of sulfur and 5s, 5p, 5d orbitals of uranium frozen. Thermochemical corrections to the electronic energies to obtain free energies were computed using the geometries and vibrational partition functions computed at the DFT level for a temperature of 25 °C and a pressure of 1 atm. Solvent effects on geometries, partition functions and energies were accounted for through the use of the conductor-like polarizable continuum model (CPCM) [6], with the United Atom Kohn Sham (UAKS) formula for the atomic radii [15]. All calculations were performed with the Gaussian 03 package [16].
2.2 Comment on the energy calculations
There is no significant difference in electronic energy of reaction calculated at the DFT and MP2 levels (Table 2) for reactions that do not involve a change in coordination number. This is not the case for reactions where the coordination number changes from 5 to 6, in which the DFT values are 10–20 kJ/mol higher than the MP2 values. The reason for this is known and related to the known tendency of current exchange correlation functions to underestimate the stability of high coordination numbers. These errors may compensate when comparing two isomers that differ by the number of bound water molecules and the mode of coordination of the sulfate group. Further details have been discussed by Rotzinger [17] for transition metals and by Wåhlin et al. [18] in the context of the uranyl aqua ion. It is thus important to compare DFT-based numbers for energetics to those obtained with wave-function-based correlated methods, such as MP2. The difference between the electronic and Gibbs energy of reaction is strongly dependent on the contribution of low-frequency vibration modes originating from shallow potential wells. In some cases it was not possible to obtain structures with only real frequencies; in other cases, such as [UO2(SO4-chelate)2(OH2)2]·(H2O), the difference between the electronic and Gibbs energy of reaction was fairly large, around 10 kJ/mol. We conclude that the calculated free energy of reaction may have an uncertainty of about the same magnitude.
Relative electronic and Gibbs free energies in kJ/mol between the different isomers of UO2(SO4) and UO2(SO4)22− computed in the CPCM solvent using both the DFT-B3PW91 and the MP2 method, in both cases using the DFT geometry with B3PW91 functional
Chemical model | Na | DFT-B3PW91 | MP2 | ||
ΔE | ΔG0 | ΔE | ΔG0 | ||
[UO2(SO4-chelate)(OH2)3],(H2O) | 5 | 0.0 | 0.0 | 0.0 | 0.0 |
UO2(SO4-mono)(OH2)4 | 5 | 12.3 | 13.0 | 16.3 | 17.0 |
[UO2(SO4-chelate)(OH2)4] | 6 | 22.2 | 32.4 | 12.0 | 22.2 |
[UO2(SO4-chelate)(OH2)3],(H2O),Na+ | 5 | 0.0 | –b | 0.0 | –b |
[UO2(SO4-mono)(OH2)4],Na+ | 5 | 12.5 | –b | 10.3 | –b |
[UO2(SO4-chelate)2(OH2)2]2−,(H2O) | 6 | 0.0 | 0.0 | 0.0 | 0.0 |
[UO2(SO4-chelate)2(OH2)]2−,(H2O)2 | 5 | 12.7 | −9.0 | 26.4 | 4.7 |
[UO2(SO4-mono)2(OH2)3]2− | 5 | −3.0 | −3.1 | 17.8 | 17.7 |
[UO2(SO4-chelate)(SO4-mono)(OH2)2]2−,(H2O) | 5 | −0.3 | –b | 18.7 | –b |
[UO2(SO4-chelate)(SO4-mono)(OH2)3]2− | 6 | 21.5 | 26.5 | 28.9 | 33.9 |
a N is the coordination number of uranyl(VI).
b The computation of the thermal contributions was not possible because of spurious imaginary frequencies obtained with the CPCM solvent model.
2.3 Details of the energy decomposition calculations
To analyze the interaction between uranyl and monodentate sulfate groups, we have computed the total interaction energy between the uranyl fragment and one monodentate sulfate group, using the energy decomposition scheme [7] implemented in the Amsterdam density functional (ADF) package [8,9]. For these calculations, we considered a simplified molecular system involving the uranyl unit and one monodentate sulfate group, without additional water molecules in the first coordination sphere. All-electron calculations were performed using the scalar relativistic zero-order regular approximation (ZORA) – DFT method [19] with the hybrid B3LYP functional [20]. All atoms were described with the uncontracted triplet-ζ Slater-type orbital (STO) augmented by two polarization functions (TZ2P) [21]. The various contributions to the total interaction energy are given in Table S3. The geometry was optimized at the B3LYP level without symmetry constraints. We then explored the evolution of the total bonding energy between the uranyl and the monosulfate groups for various U–O–S bending angles ranging from 120° to 180°, with steps of 10°.
2.4 Quantum chemical results
The bond distances for the various isomers of the uranyl monosulfate and uranyl bis-sulfate complexes calculated at the DFT-B3PW91 level in the CPCM solvent are given in Table 1. The total energy and coordinates of the different complexes are given as supporting data, Tables S1 and S2. The only other theoretical data available is that of Craw et al. [22]. They optimized the structure of uranyl monosulfate at the Hartree–Fock level, assuming a bidentate coordination and imposing C2v symmetry on the structure. They did not consider the lower-symmetry monodentate isomer, certainly because of lack of computational capacity at that time. Our calculated equatorial distances are generally shorter, by 0.20 Å, than those obtained by Craw et al. This is a result of relaxing the symmetry constrains and of the effect of electron correlation and solvent effects, as discussed for other uranyl structures by Wahlgren et al. and Vallet et al. [23]. In the experimental test solutions, there are counterions present (Na+ or H+), and we therefore explored if the Na+ counterions present in the test solutions might influence the quantum mechanical (QM) bond distances in the mono- and bidentate isomers; we found no significant changes. It is noteworthy that it was not possible to optimize the structure of the monodentate isomer in gas phase because of proton abstraction of one non-coordinated oxygen atom in the sulfate from a neighboring coordinated water molecule.
Geometries of the monodentate (mono) and bidentate (chelate) UO2(SO4) and UO2(SO4)22− complexes optimized in the solvent using the B3PW91 functional
Chemical model | Na | U–Oyl | U–Owat | U–Ochel | U–Omono | U–Oaver | U–Schel | U–Smono |
[UO2(OH2)5]2+·(H2O) | 5 | 1.746 | 2.43(1) | – | – | 2.43(1) | – | – |
EXAFS Ref. [23] | 5 | 1.77 | 2.41 | – | – | – | – | – |
UO2(SO4-mono)(OH2)4 | 5 | 1.759 | 2.48(3) | – | 2.22 | 2.43(8) | – | 3.61 |
[UO2(SO4-mono)(OH2)4]·Na+ | 5 | 1.758 | 2.46(1) | – | 2.25 | 2.42(7) | – | 3.64 |
[UO2(SO4-chel)(OH2)3]·(H2O) | 5 | 1.765 | 2.46(3) | 2.37(2) | – | 2.42(4) | 3.08 | – |
[UO2(SO4-chel)(OH2)3]·(H2O),Na+ | 5 | 1.760 | 2.43(2) | 2.39 | – | 2.42(3) | 3.07 | – |
[UO2(SO4-chel)(OH2)4] | 6 | 1.760 | 2.52(3) | 2.41 | – | 2.51(6) | 3.13 | – |
[UO2(SO4-chel)(OH2)3]b | 5 | 1.74 | 2.64 | 2.562 | – | – | 3.08 | – |
EXAFS Ref. [4] | 5 | 1.78 | – | – | – | 2.42 | 3.12 | – |
[UO2(SO4-mono)2(OH2)3]2− | 5 | 1.769 | 2.49(3) | – | 2.29 | 2.4(1) | – | 3.67 |
[UO2(SO4-chel)(SO4-mono)(OH2)2]2−·(H2O) | 5 | 1.769 | 2.46(1) | 2.42(1) | 2.27 | 2.41(6) | 3.11 | 3.67 |
[UO2(SO4-chel)(SO4-mono)(OH2)3]2− | 6 | 1.765 | 2.58(9) | 2.47(1) | 2.38 | 2.49(7) | 3.17 | 3.74 |
[UO2(SO4-chel)2(H2O)]2−·(OH2)2 | 5 | 1.768 | 2.43 | 2.40(1) | – | 2.41(1) | 3.09 | – |
[UO2(SO4-chel)2(H2O)2]2−·(OH2) | 6 | 1.769 | 2.54 | 2.46(1) | – | 2.49(3) | 3.16(1) | – |
a N is the coordination number of uranyl(VI).
b From Ref. [22]: structure optimized at the Hartree–Fock level in gas phase.
2.4.1 Five-coordinated UO2(SO4) and UO2(SO4)22−
Two different isomers, [UO2(SO4-mono)(OH2)4] and [UO2(SO4-chel)(OH2)3] (cf. Fig. 1) were identified for UO2(SO4). The U–S distances are 3.61 and 3.08 Å, respectively, the corresponding U–Osulfate distances 2.22 and 2.37(2) Å are also significantly different. The average U–OH2 distances in the two isomers are 2.48(3) and 2.46(3) Å, slightly longer than the experimental EXAFS distance, 2.43(1) Å in UO2(OH2)52+ [23a]. The negative Gibbs energy of reaction, ΔG0 = −17.0 kJ/mol (cf. Table 2) for:
UO2(SO4-mono)(OH2)4 → [UO2(SO4-chel)(OH2)3] | (3) |
![](/chimie/article/CRCHIM_2007__10_10-11_905_0/jpg/src/tex/figures/gr1.jpg)
The structure of the two five-coordinated isomers of UO2(SO4): (a) [UO2(SO4-mono)(OH2)4] and (b) [UO2(SO4-chelate)(OH2)3],(H2O) based on DFT geometry optimizations in the CPCM solvent. The dashed line indicates a hydrogen bond where the distance Hwater–Osulfate is about 1.60 Å and the angle Owater–H–Osulfate is 168°; the water outside the first coordination sphere in UO2(SO4-chelate) forms only very weak hydrogen bonds.
There are three five-coordinated isomers of UO2(SO4)22−, [UO2(SO4-mono)2(OH2)3]2−, [UO2(SO4-chel)2(OH2)]2−·(H2O)2 and [UO2(SO4-mono)(SO4-chel)(OH2)2]2−·(H2O), cf. Fig. 2. The U–S distances in the mono- and bidentate sulfate groups are close to those in the UO2(SO4) isomers and the same is true for the U–O distances. The quantum chemical bond distances for the isomer [UO2(SO4-chel)2(OH2)]2−·(H2O)2 are in good agreement with the experimental EXAFS data [4], indicating that this is the stable isomer. The calculated Gibbs energy of reaction for reaction (4):
[UO2(SO4-chel)2(OH2)]2−·(H2O)2 → [UO2(SO4-mono)2(OH2)]32− | (4) |
[UO2(SO4-chel)2(OH2)]2−·(H2O)2 → [UO2(SO4-chel)(SO4-mono)(OH2)2]2−·(H2O) | (5) |
![](/chimie/article/CRCHIM_2007__10_10-11_905_0/jpg/src/tex/figures/gr2.jpg)
The structures of the three different five-coordinated isomers of UO2(SO4)22−: (a) [UO2(SO4-mono)2(OH2)3]−2, (b) [UO2(SO4-chelate)2(OH2)2−]·(H2O)2 and (c) [UO2(SO4-mono)(SO4-chelate)(OH2)2(H2O)2−]·(H2O), based on DFT geometry optimizations in the CPCM solvent. The dashed line indicates a hydrogen bond where the distances Hwater–Omonosulfate and Hwater–Ochel-sulfate are about 1.60, 1.94 Å in structures (a), (b) and (c), respectively. The angles Owater–H–Omonosulfate and Owater–H–Ochel-sulfate are 168° and 177°, respectively.
2.4.2 Six-coordinated UO2(SO4) and UO2(SO4)22−
There is only one stable six-coordinated isomer of UO2(SO4), [UO2(SO4-chel)(OH2)4], cf. Fig. S1a; all bond distances are longer than in the corresponding five-coordinated isomer, the U–S distance is 3.13 Å, the U–Ochel, 2.41 Å and the U–OH2 2.52(3) Å. The average U–O distance in the first coordination sphere is 2.51(6) Å, larger than the value for the five-coordinated isomer, 2.42 Å, which is in good agreement with the experimental EXAFS distance, 2.40–2.44 Å [4].
There are two stable six-coordinated isomers of UO2(SO4)22−, [UO2(SO4-chel)2(OH2)2]2−·(H2O) and [UO2(SO4-chel)(SO4-mono)(OH2)3]2−, cf. Fig. S1b,c. All bond distances in these isomers are significantly longer than in the corresponding five-coordinated isomers and neither of these structures is consistent with the experimental EXAFS data. The positive Gibbs energy for reaction:
[UO2(SO4-chel)2(OH2)2]2−·(H2O) → [UO2(SO4-chel)(SO4-mono)(OH2)3]2−, ΔG0 = 33.9 kJ/mol, |
at the MP2 level, indicates that the isomer [UO2(SO4-chel)2(OH2)2]2−·(H2O) is the most stable one and confirms that bidentate sulfate coordination is preferred in six-coordinated isomers.
2.4.3 The relative stability of five- and six-coordinated isomers
Five-coordination is preferred for the UO2SO4 isomer as indicated by the negative Gibbs energy of reaction for:
[UO2(SO4-chel)(OH2)4] → [UO2(SO4-chel)(OH2)3]·(H2O), ΔG0 = −22.2 kJ/mol | (6) |
[UO2(SO4-chel)2(OH2)2]2−·(H2O) → [UO2(SO4-chel)2(OH2)]2−·(H2O)2 | (7) |
3 Discussion
The interpretation of experimental EXAFS and LAXS data rests both on assumptions of the chemical composition of the tests solutions and on the physical models used. The chemical composition of the test solutions, that is the stoichiometric concentrations of the various species present, can be estimated from published equilibrium constants. There is extensive information in the literature on the equilibrium constants formed in the uranyl(VI) sulfate system [1] and the species UO2(SO4) and UO2(SO4)22− are well established. The formation of UO2(SO4)34− has been suggested by Ahrland [24], but this proposal has not been accepted in the reviews [1a,1b]. The equilibrium constants for the complex formation reactions depend on the ionic strength and the ionic medium composition, as described in the reviews by Grenthe et al. [1a,25] and Guillaumont et al. [1b]. For the reactions:
UO22+ + HSO4− ⇌ UO2(SO4) + H+ | (8) |
UO22+ + 2HSO4− ⇌ UO2(SO4)22− + 2H+ | (9) |
The most straightforward experimental method to determine if sulfate is coordinated through one or two oxygen donors is to determine the U–S distance and this method has been used both by Moll et al. [4] and Neuefeind et al. [5]. Moll et al. find a U–S distance of 3.12 Å in test solutions A, B, and C, indicating that the sulfate is bidentate coordinated in all sulfate complexes present. The number of U–S distances in test solution C is close to 2, confirming the predominance of UO2(SO4)22−. Ahrland [24] has suggested that the complex UO2(SO4)34− is formed at high sulfate concentrations:
UO2(SO4)22− + SO42−⇌UO2(SO4)34− | (10) |
Neuefeind et al. [5] used LAXS data from two test solutions containing the same concentration of uranyl(VI), the first without and the second with sulfate present. In this way they could to a large extent eliminate the contributions of the solvent in the difference radial distribution functions, allowing a clear identification of the U–S distance and the difference in the average U–O distances in the two complexes. They find a U–S distance of 3.67 Å, in excellent agreement with the QM value, 3.64 Å, for a complex with monodentate sulfate coordination. In both studies [4,5] there is no evidence of a mixture of mono- and bidentate coordination; however, the experimental method is not sensitive enough to detect small amounts of the two isomers.
In the EXAFS study of Moll et al. [4], only one U–Oeq distance, 2.42 (2) Å, was found for test solutions A, B and C; Oeq denotes the oxygen atoms in the equatorial plane of the uranyl sulfate complex(es). This distance is equal to the average U–Oeq distance, 2.42 Å in the five-coordinated complexes obtained using QM, but significantly shorter than the corresponding average distance, 2.51 Å, computed for the six-coordinated complexes. The average U–OH2 distance in the five-coordinated complexes, 2.46 Å, is somewhat longer than in aqua ion UO2(OH2)52+, 2.43 Å, obtained in this study using the B3PW91 functional and 2.41 Å, obtained from the experimental EXAFS data [23].
The average difference between the U–O distances in the equatorial plane of UO22+(aq) ≡ UO2(OH2)52+ and in the sulfate test solutions used by Neuefeind et al. [5] is 0.017 Å. This value is very close to the difference between the average U–O distances in [UO2(OH2)5]2+·(H2O) and [UO2(SO4)(H2O)4]·(H2O), 0.008 Å, obtained using the QM optimized structures in the PCM solvent (see Table 1) where systematic errors in the distances will cancel to a large extent. It is not clear to us if the difference method used by Neuefeind et al. is precise enough to identify individual variations in the bond distances, in particular the appearance of one shorter U–O distance in the presence of four longer ones in the uranyl(VI) sulfate solution. The quantum chemical results show a significantly shorter U–O distance, 2.22 Å, in the monodentate sulfate isomer than in the chelate, 2.37 Å, for UO2(SO4). The trend is the same in the various isomers of UO2(SO4)22− and we find it unlikely that this difference is a computational artifact. In this context it is worth noticing that the U–Omono distance is significantly shorter, about 0.15 Å, than the U–Ochel distance in KNa5(UO2)(SO4)4·H2O [3a] and in K4(UO2)(SO4)3 [3b], which contain both mono- and bidentate coordination.
Both the average U–Oeq distances and the Gibbs energy of reaction suggest that the uranyl(VI) moiety has five donor atoms in the equatorial plane; this is consistent with the fact that there are no known crystal structures in the uranyl(VI) sulfate system with six-coordination [26].
The U–Osulfate–S angle is 141.5° and 143° in [UO2(SO4-mono)(OH2)4] and [UO2(SO4-mono)2(OH2)3], respectively. These values agree with that obtained by Neuefeind et al. [5], and is in the range of the values, 139° to 146°, with a grand average of 143°, found in crystal structures. Neuefeind et al. [5] suggested that the bent angle might be either an intrinsic property of uranyl–sulfate bonding, or due to hydrogen bond interactions. A screening of structure of crystals with unidentate sulfate and transition metal ions or alkali ions indicate that the M–O–S angle is always bent, around 130° or slightly larger; this property does not seem to be specific for the bonding to uranium. To rationalize this observation, we have probed the interaction between the uranyl unit and a unidentate sulfate group through the energy decomposition scheme [7] in ADF [8,9] (see Section 2.1) to determine whether the angular preference is driven by electronic or steric effects. The total bonding energy (EB) between the two fragments (the uranyl moiety and the sulfate group) was calculated for various bending angles ranging from 120° to 180° with steps of 10°. It is decomposed into a steric component (sum of the electrostatic and Pauli terms) and an electronic component (so-called orbital interaction). Fig. 3 shows the variation of the various contributions with the U–O–S angle. To facilitate the comparison, we have assigned the binding energy 0 kJ/mol at 180° to each of the contributions and calculated the energy at other angles relative to that reference. From the shape of the curve of the orbital interaction (electronic) component, it is clear that there is an electronic preference for linearity. However, steric interactions favor the bent structure, because the classical electrostatic interaction offsets the unfavorable Pauli repulsion upon bending. A plot of the total interaction energy indicates that there is an energy minimum around 125° and that the energy gain upon bending is about 16 kJ/mol. The bond angle is smaller than observed experimentally, but this might very well be due to the fact that the calculations were made using a model without water in the first coordination sphere. The key point is that the “driving force” for a bent U–O–S angle is the steric interaction energy, which overcomes the electronic preference for a linear U–O–S angle. Such steric interactions are most likely also important for the coordination of monodentate sulfate in other metal complexes and thus provide a rationalization for the general observation of bent M–O–S bonds.
![](/chimie/article/CRCHIM_2007__10_10-11_905_0/jpg/src/tex/figures/gr3.jpg)
Contributions to the total bonding energy of UO2(SO4-mono) for different U–O–S angles: (a) the steric contribution is the sum of the electrostatic interaction and Pauli repulsion; (b) the total bonding energy is the sum of the steric contributions and orbital interaction and has a minimum around 125°.
It is very clear from the experimental observations that the composition of the test solutions, (total concentration of uranium, ionic strength, and ionic medium) results in different modes of coordination of sulfate. This indicates that the difference in Gibbs energy between the two isomers is small, as also found in the quantum chemical calculation (Table 2) that refers to reactions (4)–(6). These reactions involve a transfer of water from the first to the second coordination sphere of U(VI), while in the test solutions the reactions describe the transfer of water from the first coordination sphere to the bulk solvent, as exemplified by:
UO2(SO4-mono)(OH2)4⇌UO2(SO4-chel)(OH2)3 + H2O | (11) |
Reaction (11) shows that the relative stability of the mono- and bidentate isomers depends on the water activity, where a low water activity will favor the chelate. As the water activity decreases with increasing ionic strength and temperature, it seems possible to test this by performing LAXS and EXAFS studies at higher temperatures. A very small change in Gibbs energy of reaction for reaction (11) will result in a substantial change in the relative amounts of the two isomers, as shown in Table 3. With the approximations used in the chemical and quantum chemical models, it seems difficult to claim a higher accuracy than 10 kJ/mol in the Gibbs energy and this uncertainty range spans the range of ΔG0 values in Table 3, for which there is a substantial change in the relative amounts of the isomers.
Ratio between the concentrations [UO2(SO4-chel)(OH2)3]/[UO2(SO4-mono)(OH2)4] at different values of the Gibbs energy of reaction
ΔG0 (kJ/mol) | [UO2(SO4-chel)(OH2)3]/[UO2(SO4-mono)(OH2)4] |
−5.0 | 7.5 |
0 | 1.00 |
+5.0 | 0.133 |
4 Conclusions
Standard solution chemical methods give information on the stoichiometry of the complexes with respect to the reactants, in this case UO22+ and SO42−, but not on the mode of coordination of the ligand and the number of coordinated water ligands. Structure chemical methods such as EXAFS provide this information, but the structure models are not always unique. In addition, EXAFS methods require information on the speciation, i.e. the concentration of the different species in the test solutions used. In systems where several different complexes occur in comparable concentrations, it is very difficult to deduce structure information. Quantum chemistry provides tools to obtain both the geometry and relative energy of different complexes and this information, combined with experimental solution chemical and EXAFS data, makes it possible to obtain a better understanding of the chemistry of aqueous complex formation than using the experimental or theoretical methods alone. In the present system we conclude that there are several possible isomers of the complexes UO2(SO4) and UO2(SO4)22− that have very near the same energy. These isomers consist both of species where the sulfate ligand is coordinated in bi- and monodentate fashion. However, the number of donor atoms in the equatorial plane of the UO2 moiety is five in all isomers, as supported both by experimental and calculated average bond distances and the relative energy of isomers with the same mode of coordination of sulfate, but a different number of coordinated water ligands. The relative stability of the different isomers is small and the mode of coordination of the sulfate ligand seems therefore to be dependent on the water activity, as also noted experimentally.
Acknowledgments
This study was supported by the ACTINET network of excellence and through a grant from the Trygger Foundation. Computational resources have been provided by the ‘Institut de développement et de ressources en informatique scientifique’ of the ‘Centre national de la recherche scientifique’, IDRIS (Project 61859), the ‘Centre informatique national de l'enseignement supérieur’ (CINES, Project phl2531), and the ‘Centre de ressources informatiques’ (CRI) of the ‘Université des sciences et technologies’ of Lille. CERLA is supported by the ‘Ministère chargé de la recherche’, the ‘Région Nord-Pas-de-Calais’, and the FEDER.
Appendix
The experimental EXAFS and LAXS investigations have been made in test solutions that contain UO22+(aq), UO2(SO4)(aq) and UO2(SO4)22−(aq), and it is important to know these concentrations for the analysis of the experimental data. Such calculations are based on the known total concentrations of the components and known equilibrium constants for the formation of the two sulfate complexes. The latter are tabulated in data bases [1] and usually refer to standard state conditions; they have therefore to be recalculated to the actual ionic medium of the test solutions. We have reanalyzed the composition of the test solutions described by Moll et al. [4] using the specific ion interaction theory (SIT), as described in the example given below. The theoretical basis for the SIT is described in Refs. [1,25].
The activity coefficient of an ion, i, is described by Eq. (A.1):
(A.1) |
The activity correction is particularly simple in sulfuric acid solutions, as shown by the following example for the reactions:
UO22+ + HSO4− ⇌ UO2(SO4) + H+ | (A.2) |
(A.3) |
(A.4) |
(A.5) |
(A.6) |
(A.7) |
The equilibrium constant for the reaction:
UO22+ + 2HSO4− ⇌ UO2(SO4)22− + 2H+ | (A.8) |
The specific ion interaction parameters used are taken from Ref. [1b] and are listed in Table A.1.
Using the equilibrium constants at zero ionic strength from Ref. [1] (log K0(A.2) = 1.17 and log K0(A.8) = 0.18) and the interaction coefficients from Table A.1, we obtain the following equilibrium constants valid in 1.0 M H2SO4: log K(A.2) = 0.71 and log K(A.8) = 0.56; the uncertainty in these quantities is at least 0.10 log units. The composition of test solution A in the paper by Moll et al. [4] is then [UO22+(aq)] = 5.6 ± 1.4 mM, [UO2(SO4)(aq)] = 26.5 ± 0.2 mM and [UO2(SO4)22−] = 17.8 ± 1.7 mM; the uncertainty ranges have been estimated using the uncertainty in the equilibrium constants. In their calculation of the composition of test solution A, Moll et al. do not seem to have taken the formation of UO2(SO4)22− into account. It seems clear that the predominant complex in the test solution is UO2(SO4)(aq), but with a significant amount of UO2(SO4)22−. It is not possible to make an accurate estimate of the speciation in test solution B because of the high ionic strength in 10 M H2SO4; however, as a result of the higher total concentration of HSO4−, it seems clear that the predominant complex is UO2(SO4)22−. It is not possible to make an accurate analysis of the speciation in test solution C using the SIT, as the test solution contains large amounts of both ClO4− and HSO4− and has a high ionic strength. To conclude: the reanalysis of the test solutions used in the EXAFS measurements of Moll et al.[4] are only to a minor degree affected by the presence of UO22+(aq), and the predominant complexes are UO2(SO4)(aq) and UO2(SO4)22−.
Specific ion interaction coefficients [1b] for the estimate of activity coefficients in the uranyl(VI)–H+–SO42− system
Ion interaction coefficient | E(i, j) (kg mol−1) | Comments |
ɛ(UO22+, HSO4−) | 0.35 ± 0.11 | The interaction coefficient is estimated from the known values ɛ(UO22+, ClO4−) = 0.46 and ɛ(UO22+, NO3−) = 0.24. |
ɛ(HSO4−, H+) | 0.15 | The interaction coefficient is estimated from the known value ɛ(HSO4−, Li+) = 0.15. |
ɛ(H+, HSO4−) | 0.15 | ɛ(HSO4−, H+) ≡ ɛ(H+, HSO4−) |
ɛ(UO2(SO4)22−, H+) ≈ ɛ(UO2(SO4)22−, Li+) | −0.02 ± 0.06 | The interaction coefficient is estimated by increasing the known value ɛ(UO2(SO4)22−, Na+) = −0.12 ± 0.06 by 0.10, the difference found between Li+ and Na+. |