Quantum mechanical calculation of aqueuous uranium complexes: carbonate, phosphate, organic and biomolecular species.

BACKGROUND: Quantum mechanical calculations were performed on a variety of uranium species representing U(VI), U(V), U(IV), U-carbonates, U-phosphates, U-oxalates, U-catecholates, U-phosphodiesters, U-phosphorylated N-acetyl-glucosamine (NAG), and U-2-Keto-3-doxyoctanoate (KDO) with explicit solvation by H2O molecules. These models represent major U species in natural waters and complexes on bacterial surfaces. The model results are compared to observed EXAFS, IR, Raman and NMR spectra.
RESULTS: Agreement between experiment and theory is acceptable in most cases, and the reasons for discrepancies are discussed. Calculated Gibbs free energies are used to constrain which configurations are most likely to be stable under circumneutral pH conditions. Reduction of U(VI) to U(IV) is examined for the U-carbonate and U-catechol complexes.
CONCLUSION: Results on the potential energy differences between U(V)- and U(IV)-carbonate complexes suggest that the cause of slower disproportionation in this system is electrostatic repulsion between UO2 [CO3]3(5-) ions that must approach one another to form U(VI) and U(IV) rather than a change in thermodynamic stability. Calculations on U-catechol species are consistent with the observation that UO2(2+) can oxidize catechol and form quinone-like species. In addition, outer-sphere complexation is predicted to be the most stable for U-catechol interactions based on calculated energies and comparison to 13C NMR spectra. Outer-sphere complexes (i.e., ion pairs bridged by water molecules) are predicted to be comparable in Gibbs free energy to inner-sphere complexes for a model carboxylic acid. Complexation of uranyl to phosphorus-containing groups in extracellular polymeric substances is predicted to favor phosphonate groups, such as that found in phosphorylated NAG, rather than phosphodiesters, such as those in nucleic acids.

Background

The toxicity and radioactivity of U makes it a potentially hazardous element in the environment. In areas of high U concentrations, understanding U chemistry is imperative in order to predict its fate, transport, and risk. Uranium is capable of forming a wide variety of aqueous and surface complexes. Furthermore, redox reactions, mainly between U(VI) and U(IV), are common in subsurface environments (e.g., [1]).

Research has focused on the environmental chemistry of U with the goal of managing and remediating U-contaminated sites in the most effective manner ([2-4] and references therein). Recent studies have probed the molecular-level structures and processes that influence the overall behavior of U in the environment (e.g., [5]). Both analytical and theoretical studies have discussed complexation with numerous ligands [6-9] and the redox reactions between U(VI) and U(IV) (e.g., [10-15]). Computational chemistry is an important complement to experimental studies of U chemistry because this methodology can provide information that is not available via experiment, especially for transient species and those with short kinetic lifetimes. In order for the molecular modeling to be useful however, one must demonstrate that the computational methodology produces accurate results compared to known experimental data.

Before one can simulate structures, thermodynamics and kinetics with confidence, a computational methodology must be tested against observation. Environmental chemists are interested in U complexation and redox reactions, so this study focused on evaluating the ability of quantum mechanical calculations to reproduce experimental data on aqueous U complexes and redox chemistry. Specifically, models of aqueous U(VI), U(V), and U(IV) were generated and compared with experiment and previous calculations. Uranium complexes with inorganic (carbonate and phosphate), organic (oxalate and catechol), and biological (phosphodiester, phosphorylated glucosamine, and the 2-Keto-3-deoxyoctanoate) ligands were modeled and analyzed in light of previous observations. The model results are compared to interatomic distances from EXAFS, observed vibrational frequencies, and 13C and 17O NMR chemical shifts. Calculations on the observed oxidation of catechol by U(VI) are also presented.

Experimental

Computational

Hybrid density functional theory calculations were performed on all model systems using the program Gaussian 03 [16]. The basis set 6-31G(d,p) [17-20] was used for H, C, and O and the Stuttgart pseudopotential ECP60MWB and the corresponding ECP60ANO basis set [21,22] were used for U. This relativistic pseudopotential uses 60 electrons as the "core" electrons and 32 as the valence electrons. The Becke 3-parameter exchange [23,24] and Lee, Yang and Parr [25] correlation functionals were used for energy minimizations, frequency analyses and Gibbs free energy calculations. The Hartree-Fock method was used for NMR chemical shielding calculations. Excellent results were obtained by de Jong et al. [6] using a similar method. All atoms were allowed to relax during energy minimizations, and no symmetry constraints were applied.

The models were created including explicit H2O molecules around the complex to account for H-bonding by aqueous solutions. In this paper, a H-bond is considered to exist if the H---O distance is less than or equal to 2.0 Å and if the O-H---O angle is greater than 120°. These criteria are similar to those used by others (e.g., [26]) and are useful for identifying significant shifts in the calculated O-H stretching frequencies [27]. In general, initial models of solvation were created by positioning H2O molecules with either their H or O atom at approximately 1.8 Å from a O or H atom on the solute model with a O-H---O angle between 120 and 180°.

Previous work [28,29] has shown that including the H2O molecules in the primary solvation shell of UO22+ is important for obtaining accurate structures, vibrational frequencies and energetics. This study (as in [29]) investigates the effects of adding a second solvation shell to the hydrated UO22+ cation. The number of H2O molecules was chosen to be at least the minimal number necessary to form one H-bond to each of the possible H-bonding atoms in the U coordination sphere (e.g., 2 H2O molecules for each U-OH2 group). In some cases (e.g., UO2-oxalate), an increasing number of H2O molecules were included in the model to assess the effects of explicit solvation on the predicted interatomic distances, vibrational frequencies, and NMR chemical shifts. Energy minimizations were generally carried out with the default criteria in Gaussian 03. When imaginary frequencies were calculated from an energy minimized structure, a re-optimization of the structure was performed with the "Opt = Tight" option until a structure with no imaginary frequencies was found. Although we have obtained potential energy minima, there is no guarantee that each configuration is the global minimum because the potential energy surface of these models will be complicated due to many possible H-bonding configurations. Any energy minimization scheme is unlikely to find the global potential energy minimum, so molecular dynamics simulations would be useful in the future to investigate configuration space at the temperature of interest and determine average configurations for these models.

Calculated results were compared to observed EXAFS, IR, Raman and NMR spectra. Model interatomic distances were directly compared to the values extracted via analysis of EXAFS data and vibrational spectra. Theoretical vibrational frequencies were compared to observed values for uranyl without a scaling factor applied because the appropriate value is not known for this computational methodology. For the vibrations of the ligands, a scaling factor of 0.96 was applied as determined by Wong [30] for B3LYP/6-31G(d) with the assumption that the p-functions added to the H atoms do not significantly affect the vibrational frequencies of the C-C and C-O bonds. This assumption is justified based on the work of Andersson and Uvdal [31] who showed scaling factors varying only between 0.96 and 0.97 for the 6-311G basis set with addition of polarization and diffuse functions.

Chemical shieldings were calculated with the gauge-including atomic orbital (GIAO) method [32]. The 13C chemical shift values, δ13C, were calculated relative to the chemical shielding calculated for C in tetramethyl silane (TMS) using a model TMS structure energy minimized with the B3LYP/6-31G(d,p) method to be as similar to the uranium complex results as possible. The 17O reference values were the average of O shielding values in a 19 H2O cluster which worked well in the past for predicting 17O chemical shifts [33] and (CD3)2CO which has been used experimentally as well [34]. The 13C and 17O chemical shieldings were calculated using the HF/6-31G(d,p) method [17] for both the TMS and uranium complexes. Tests were run using B3LYP/6-31G(d,p) and B3LYP/VTZ [35] for the aqueous UO22+ and uranyl-oxalate complexes, but the agreement with observed δ13C values was not as good as for the HF/6-31G(d,p) method. However, the B3LYP/VTZ method resulted in significantly better agreement with the observed δ17O values. In general, the accuracy of calculated chemical shifts with different methods vary among elements (e.g., [36]), but a reason as to why this should be is a matter for further investigation and beyond the scope of this paper. There are also issues especially regarding the electric field gradient when including U in the NMR calculations due to the relativistic nature of the electrons in this atom [37,38].

Atomic charges were calculated with the natural bond orbital analysis (NBO) method of Glendening et al [39]. A modified version of natural population analysis was used for the explicitly solvated uranyl-catechol models that includes U 6d electrons in the valence space [40].

Starting structures were based on previous interpretations of experimental data where available [41-47]. The initial aqueous U models were UO22+(OH2)6•20(H2O), UO22+(OH2)8•18(H2O), UO2+OH(OH2)5•18(H2O), UO2(OH)2(OH2)4•18(H2O), UO2+(OH2)6•20(H2O), and U4+(OH2)10•20(H2O) (in both singlet and triplet states). In this notation, the first solvation sphere is denoted by "(OH2)n" and the second solvation sphere is denoted by "•n(H2O)". Thus, two initial coordination numbers (CN) were tested for U(VI): the first with CN=8 and the second with CN=10. Comparisons of the calculated structures and potential energies generated from these initial configurations were used to test which coordination state should be favored. Starting the aqueous structure in a higher coordination state is a test of whether the lower coordination state observed (i.e., CN=7) can be reproduced or whether a local minimum would be found.

Similarly, the inorganic uranium complexes UO2(OH)2(CO3)22-•28(H2O) (cis and trans), UO2(CO3)34-•28(H2O), UO2(CO3)35-•28(H2O), Ca2UO2(CO3)3•28(H2O) (two configurations), UO2(PO4)-(OH2)4•30(H2O) (mono- and bidentate), UO2(HPO4)(OH2)4•30(H2O), UO2(HPO4)22-(OH2)2•33(H2O), and UO2(H2PO4)(OH2)4+•30(H2O) were modeled. The organic uranium complexes included UO2(OH2)4(C2O4) (C2O42- = oxalate), UO2(OH2)4(C2O4)•11(H2O), UO2(OH2)4(C2O4)•30(H2O), outer-sphere UO22+(OH2)4(C6H4(OH)2)•28(H2O), singlet, triplet and quintuplet state UO22+(OH2)4(C6H4(OH)2)•28(H2O), and U(OH)2(OH2)4(C6H4O2)•28(H2O) (C6H4(OH)2 = catechol; C6H4O2 = quinone). To model extracellular bacterial ligands, complexes of uranium with a phosphodiester (UO2(OH2)4(O2P(OCH2C4H6OOH)2)+•23(H2O)=UO2OrgPO4•27(H2O), phosphorylated N-acetylglucosamine (UO2-GlcNPO4•26(H2O)), and 2-Keto-3-doxyoctanoate (UO2-KDO•26(H2O)) were studied.

Gibbs free energies of each model were calculated using the integrated equation formalism polarized continuum model (IEFPCM) of Cancès et al. [48]. The Gibbs free energies include the electrostatic energy interaction between the model and the continuum as well as the cavitation, dispersion and repulsion energy terms. The thermal and zero-point energy contributions to the Gibbs free energy as estimated in the gas-phase frequency calculations were added to the IEFPCM Gibbs free energies to correct for temperature.

Raman

A Nicolet Almega model dispersive Raman spectrometer from Thermo Scientific with a 785 nm laser was used to obtain spectra from 3446 to 111 cm-1. The final spectra for each sample were the result of 128 scan averaging at high resolution (4.8–8.9 cm-1) with a 3.1 micron spot size. Liquid Raman spectra were collected with 100% laser intensity and were baseline corrected by subtracting a water spectrum from the background.

NMR

Solid-state 13C CP/MAS NMR spectra of a freeze-dried 1:1 uranyl-catechol solution were obtained at 100.6 MHz with a 400 MHz Varian Inova spectrometer with an optimized tangential ramp of the transverse 13C field over a contact time of 2 ms. The sample was spun at 8.0 kHz in a 5 mm Varian Jakobsen-style probe assembly. The 13C chemical shifts are referenced to TMS via an external adamantane standard set to δ13C = 38.6 ppm. After re-dissolution of this sample in DMSO-d6 solution-state 13C and 1H NMR spectra were obtained on a 600 MHz Varian Inova spectrometer fitted with a 10 mm broadband probe, using pulse delays of 1 and 10 s respectively.

Results

U aqueous

The aqueous species of U(VI), U(V) and U(IV) were modeled to test the ability of the computational methodology to predict accurate interatomic distances as compared to those derived via EXAFS spectroscopy. The various bond lengths are listed in Table 1[49,50]. In general, observed and calculated interatomic distances are in good agreement. The U(VI)aq model is discussed first followed by U(V)aq then U(IV)aq.

Table 1

Coordination numbers, interatomic distances (in Å) and aqueous free energies (G in Hartrees/molecule) of inorganic aqueous U species.

Model	CN	U=O	U-O(H2)	U-OH	G
Aqueous Uranium
UO2H2O6•14(H2O)	7	1.79	2.44	-----	-1730.2813
"(6-311+G(d,p))	7	1.77	2.47	-----	-----
UO2OHH2O5•14(H2O)	7	1.80	2.44	2.19	-1729.8137
UO2OH2H2O4•14(H2O)	7	1.83	2.43	2.23	-1729.3376
U(V)O2H2O6•14(H2O)	7	1.87	2.52	-----	-1730.5124
Expt [41]	1.76 ± 0.02	2.41 ± 0.02	-----	-----
U(IV)H2O8•22(H2O) singlet	8	-----	2.42	-----	-2343.7790*
U(IV)H2O9•21(H2O) triplet	9	-----	2.47	-----	-2343.8518
Expt [41]	-----	2.41 ± 0.02	-----	-----
Carbonates
	CN	U=O	U-O(H2)	U-O(C)	U---C	U-Ca	G
UO2H2O2(CO3)2•28(H2O)
-trans	8	1.80	2.63	2.47	2.93	-----	-3022.6274
-cis	8	1.80	2.62	2.46	2.94	-----	-3022.6241
UO2H2O2(CO3)2•28(H2O)
6-31+G(d,p) -trans	8	1.80	2.66	2.47	2.94	-----	-3023.6182
6-31+G(d,p) -cis	8	1.79	2.67	2.47	2.93	-----	-3023.6268
Expt [42]	11	1.80	-----	2.43	2.88	-----	-----
				4.17 (distal)
UO2(CO3)3•28(H2O)	8	1.82	-----	2.50	2.95	-----	-3133.8283
				4.20 (distal)
UO2(CO3)3•28(H2O)	8	1.82	-----	2.52	2.97	-----	-3134.8247
6-31+G(d,p)				4.23 (distal)
Expt [43]	8	1.80	-----	2.43	2.89	-----	-----
				4.13 (distal)
Expt [49]	8	1.81	-----	2.44	2.90	-----	-----
Expt [50]	8	1.81	-----	2.44	2.92	-----	-----
				4.23 (distal)
UVO2(CO3)3•28(H2O)	8	1.89	-----	2.61	3.05	-----	-3134.0109
				4.30 (distal)
Expt [43]	8	1.90	-----	2.50	2.94	-----	-----
Expt [44]	8	1.94	-----	2.47	2.90	-----	-----
Expt [50]	8	1.91	-----	2.50	2.93	----	-----
Ca2UO2(CO3)3•28(H2O)A	8	1.84	-----	2.46	2.93	3.68	-4488.8621
Ca2UO2(CO3)3•28(H2O)B	8	1.81	-----	2.49	2.95	4.05	-4488.9081
Expt [45]	8	1.78	-----	2.43	2.86	4.07	-----
Phosphate
	CN	U=O	U-O(H2)	U-O(P)	U---P	G
UO2PO4•30(H2O)	7	1.80	2.49*	2.40	3.10	-3443.3837
UO2HPO4•30(H2O) mono	--	-----	-----	-----	-----	-3520.2502§
UO2HPO4•30(H2O) bi	7	1.79	2.48	2.445	3.09	-3520.6344
UO2(HPO4)2•33(H2O) bi	7	1.80	2.45	2.50	3.11	-4163.0025
" mono**	"	"	"	2.22	3.78	"
UO2(HPO4)2•33(H2O) bi	8	1.79	2.49	2.55	3.18	-4163.0409
" bi#	8	1.80	2.58	2.52	3.16	-4163.0190
UO2H2PO4•30(H2O) mono	7	1.78	2.51	2.28	3.64	-3520.2679
" bi	7	1.79	2.46	2.46	3.11	-3520.2712
Expt [46]	6	1.78	2.49	2.30	3.16	-----

* – One U-OH2 deprotonated to form a U-OH at 2.28 Å and a H3O+ in the solvation sphere

§ – No frequency analysis was performed for this unstable configuration.

** – The UO2(HPO4)2 complex has one HPO42- group attached in a bidentate configuration and the second in a monodentate configuration. This entry accounts for the U-O(P) and U---P distances in the monodentate HPO42- ligand.

# – Second configuration with O atoms of HPO42- bonded to UO22+ in parallel rather than perpendicular as in the lower energy configuration

U(VI) aqueous

The experimental U=O and U-OH2 distances of U(VI)aq are 1.76 ± 0.02 and 2.41 ± 0.02 Å, respectively, with a coordination number (CN) of 7 to 8 counting both the axial and equatorial O atoms [41,51-53], and the calculated values are 1.79 and 2.44 Å, respectively, with a CN=7 (i.e., 2 axial doubly bonded O atoms and 5 H2O molecules). Interpretations of the experimental data have also been confirmed by high-level quantum calculations with continuum and explicit solvation by Gutowski and Dixon [29]. These model first solvation shell values are also similar to DFT calculations [54] although the model results presented in this study are slightly more accurate. Second solvation shell U---O distances have been observed at 4.5, 7.0 and 8.7 Å [55]; the second solvation sphere average U---O distance calculated in UO22+(OH2)6•20(H2O) is 4.47 Å in excellent agreement with observation.

Gal et al. [56] observed IR peaks at 963, 253 and 160 cm-1 and Raman peaks at 874 and 198 cm-1 for this species associated with uranyl motion. Model frequencies involving U-OH2 motion were predicted at 152 to 188 cm-1 (these modes were predicted to have minimal IR intensity; Fig. 1a), UO22+ motion at 246 to 267 cm-1, U=O stretches at 833 to 848 cm-1 (only the 848 cm-1 mode had significant IR activity), and U=O stretches at 949 cm-1. A number of vibrations in the 150 to 300 cm-1 range have uranyl motion character, but the Raman activity of these modes is weak. The only strong Raman vibration is a triplet predicted at 834, 840 and 848 cm-1 (Fig. 1b). Peaks at 879 and 895 cm-1 have minor motion of the uranyl ion and are weak Raman scatterers. The present results are similar to previously calculated U=O stretching frequencies (≈820 and 890 cm-1 – [54]; 915 and 1010 cm-1 – [29]). This correspondence is considered a fair reproduction of the observed frequencies, and therefore the model is a reasonable representation of the solvation environment.

Figure 1

Calculated (a) IR and (b) Raman spectra in the 0 to 2000 cm. Model IR spectra of (c) UO2(CO3)34-•28(H2O) and (d) UO2(PO4)-(OH2)4•30(H2O) complexes exhibit reasonable correspondence to observed vibrational frequencies.

The observed δ17O in UO22+(aq) is between 1110 and 1120 ppm (depending on U concentration and relative to water) for axial O atoms and -482 ppm (relative to (CD3)2CO) for equatorial O atoms at 25°C [34]. The calculated δ17O values in this study are 1036 and 1147 ppm (average = 1092 ppm) for the axial O atoms relative to a 19 H2O cluster and -509 to -568 ppm (average = -539 ppm) for the equatorial O atoms. The absolute accuracy of the calculated values is not as high as one would like (23 and 57 ppm), but the percentage errors are 2% and 12%, respectively.

Sources of error are the choice of basis set, electron correlation method, solvation and differences in the experimental and model reference state. Using the same structure and the B3LYP/VTZ method [35], the calculated Oax and Oeqδ17O values are 1128 and -522 ppm, respectively; both of which are closer to observation. Thus, basis set and electron correlation have a more significant effect for the calculated δ17O, and this helps confirm the reasonable accuracy of the model structure.

Kubicki and Sykes [57] were able to reproduce δ17O values in aluminosilicate glasses to within a few ppm using the B3LYP/6-311+G(d,p) method, so the error here may be associated with relativistic effects due to bonding with the U atom. Solvation plays a significant role as evidenced by the large range of model δ17O values depending on details of the U-O and H-bonding. Lastly, the experimental reference states are the bulk liquid phases of water and (CD3)2CO whereas the model reference states are a 19 H2O cluster and gas-phase (CD3)2CO. The gas-phase to bulk liquid shift for 17O chemical shielding can be as large as 36 ppm [58], so this discrepancy contributes to the overal errors mentioned above as well.

Two deprotonated U(VI) monomer species were also modeled. These models are useful for estimating the behavior of aqueous U(VI) as pH increases because UO22+ undergoes hydrolysis

above pH 5 [59,60]. These species are also the likely precursors of uranium oligomers such as (UO2)2(OH)22+ that form in solution and as such can be used in future studies of uranium polymerization reactions.

As expected, the U-OH bond lengths (Table 1) are intermediate between the U=O and U-OH2 bond distances. These different bond lengths may not be distinguishable as distinct peaks in EXAFS spectra but would decrease the average U-Oeq bonds by approximately 0.05 Å. Recently, Müller et al. [61] have observed an infrared peak at 922 cm-1 in water at pH 4 and assigned this peak to monomeric UO2(OH)20. The IR-dominant peak calculated here for UO2OH2(H2O)4•14(H2O) is found at 880 cm-1 which significantly underestimates the frequency of the observed peak. However, the dominant IR peak in the UO2OH(H2O)5•14(H2O) is calculated to be at 909 cm-1 – closer to the observed peak. Since the assignment in Müller et al. [61] was based upon an interpolation of this mode between UO22+(aq) and U(OH)42-(aq), the actual assignment may be the UO2(OH)+ species.

As is commonly observed, the calculated ΔG's for hydrolysis have significant discrepancies with experiment. Using a ΔGhyd for H+ as +1098 kJ/mol [62] and the calculated G's in Table 1, the ΔGhydrolysis are equal to +130 and +282 kJ/mol as compared with the experimental values of +29 and +59 kJ/mol [6]. Such a large discrepancy is large for this level of theory and may indicate that oligomers are being formed simultaneously with hydrolysis. Because we have not modeled the oligomerization reaction, if the oligomers have a lower free energy the observed ΔGhydrolysis would be less than calculated here consistent with the discrepancy above. Oligomerization is possible depending on the concentration of uranyl in solution [61].

U(V) aqueous

A literature search did not result in any published EXAFS studies of aqueous U(V). This is likely due to the instability of this species when uncomplexed. Renshaw et al. [44] have published interatomic distance for U(V)O2CO3 complexes in solution and Chen et al. [63] have determined the crystal structure of K(UO)Si2O6 in which U is pentavalent. Mizuouka and Ikeda [64] also found stable U(V) in dimethyl sulfoxide (DMSO) solution. Using these species as a guide, one expects U=Oax distances of 1.82 to 2.07 Å, respectively. The calculated changes in these parameters for aqueous U(V) (Table 1) compared to U(VI) are consistent with this magnitude of elongation. One reason for the elongation of the U=Oax bonds in these calculations is that the H-bonding to the O atoms in these bonds has increased significantly over that found in the U(VI)aq model where only one H-bond of 1.92 Å formed. H-bonding to the Oax atoms increases the U=Oax bond length as electron density is shifted from the U=Oax bond to the H-bond. In contrast, the U(V) model had two H-bonds of 1.76 and 1.86 to one Oax atom and another of 1.58 Å to the second Oax atom.

The increase in H-bonding is a result of the increased electron density of the Oax atoms in the U(V) versus U(VI) species (the atomic charge as calculated with natural population analysis changes from -0.93 in U(VI) to -1.15 in U(V)). Although U(V)aq can still be represented as a UO2+ ion, these changes in U=Oax bond length and H-bonding foreshadow the change to a aquo-U(IV) species where only H2O molecules are bound to the U atom. The U(V)-Oeq bonds were found to be 2.47 Å in [44]. With U=Oax of 1.87 Å and U-Oeq of 2.52 Å, the model results appear to be a reasonable prediction for this species and consistent with MP2 (customized basis set) results (U=Oax = 1.81 and U-Oeq = 2.51) of Tsushima et al. [65] and the DFT (small-core ECP for U and 6-31G(d) for other atoms) results of Shamov and Schreckenbach [54] (U=Oax = 1.83 to 1.91 and U-Oeq = 2.52 to 2.59).

Vibrational spectra are available for U(V)aq, in Best et al. [66] who observed a band attributed to aqueous U(V) at 910 cm-1. No calculated vibrational modes associated with U(V)-O vibrations with significant IR intensity were found in these model calculations. Best et al. [66] noted that the presence of this band was highly sensitive to pH and that U(V) could be stabilized by dimerization with U(VI) species, but neither effect was modeled in these calculations. Hence, this assignment remains uncertain.

Our calculated IR frequencies are compared to U(V) in other phases to test whether the calculations are missing this vibration or whether the observed peak could have other origins. U=O vibrations are likely to be de-coupled from other bonds, so their frequencies should be similar (although not exactly the same) in the aqueous and crystalline phases. Chen et al. (2005) reported an IR frequency in a uranium silicate near 941 cm-1 that was assigned to a U-O stretching mode of UO2+. Mizuouka and Ikeda [64] observed an IR peak at 775 cm-1 that was assigned to asymmetric O=U=O stretching in U(V) in DMSO. A doublet with significant IR intensity exists in the model spectrum at 779 + 783 cm-1. This mode does have a component O=U=O asymmetric stretch coupled to H2O motions which was also predicted in Shamov and Schreckenbach [54]. The closest IR peak with dominantly O=U=O asymmetric stretching motion in the theoretical spectrum is predicted to occur at 826 cm-1. This vibration was predicted to occur at 786 cm-1 in a UO2+(OH2)5•12(H2O) model by Shamov and Schreckenbach [54], and this model frequency was noted to change significantly with solvation. The fact that the observed modes in these various phases have close counterparts in the U(V)aq model calculations suggest that our assumption of de-coupling for these modes is justified; however, the model calculations fail to predict any peaks for U(V)aq near 910 or 949 cm-1, so the extension of the U=O bond from the U(VI) to U(V) species may be exaggerated in the model results, thus lowering the bond force constant more than observed.

Although calculated frequencies generally are scaled for comparison with experiment to account for anharmonic and basis set effects, the correct scaling factor for uranyl vibrations and the method used here is not known, so no scaling was used. However, the close agreement of the observed and model results suggests that the model realistically represents the UO2+ ion, especially considering that the aqueous UO2+ model is being compared to a crystalline uranium silicate or U(V) complexes in DMSO.

U(IV) aqueous

The observed U-OH2 bond lengths for aqueous U(IV)are 2.41 ± 0.02 Å [41,67], and there are 9 U-O(H2) bonds averaging 2.47 Å from these computations (Table 1). (Note: This value was found for the triplet state. In the singlet state, U(IV)aq is predicted to have 8 U-O(H2) bonds averaging 2.42 Å, but the energy is approximately +190 kJ/mol higher in the singlet state.). A discrepancy on the CN exists, however, in that Moll et al. [67] concluded that U(IV) should be 10-fold coordinated.

In consideration of this interpretation, the model U(IV)aq species was generated with 10 H2O molecules initially bonded to the U(IV) cation. A minimum energy configuration for a CN of 10 could not be found. Attempts were made to find a stable U(IV) with 10 H2O molecules bonded using B3LYP/6-311+G(d,p) and with the LANL2DZ/6-31G(d,p) methods, but the O---O distances in these complexes were so short that a stable electron density could not be found. The resulting optimized structure for the B3LYP/6-311+G(d,p) U(IV)aq was similar to the B3LYP/6-31G(d,p) structure and experiment (Hennig, 2005), so we conclude that the predicted structure is not a strong function of basis set above 6-31G(d,p) (U-O(H2) = 2.44 and 2.42 Å, respectively). The LANL2DZ/6-31G(d,p) structure resulted in average U-O(H2) distances of 2.10 Å, so we conclude that the LANL2DZ is not as accurate as the ECP60MWB. In all three cases, the 9-fold coordinate U(IV)aq resulted in a stable energy minimum with no imaginary frequencies and the 10-fold coordinate U(IV)aq was unstable.

One would expect that these energy minimizations at 0 K would tend to overestimate the stability of higher coordinate species because thermal energy would favor the release of H2O from the coordination sphere which has a positive enthalpy [68]. Energy minimization for this model within a self-consistent reaction field could be performed to test whether long-range solvation could stabilize the higher coordination state.

UO2-Inorganic complexes

Uranyl complexes with a variety of inorganic anions (e.g., carbonate, phosphate, silicate). The uranyl-carbonate complexes are particularly important in the environment (e.g., [69,70]), so this study focused mainly on these U-inorganic complexes. In addition, phosphate or phosphoryl groups may be important in the chemistry of uranium with bacterial surfaces [71], so the uranyl-phosphate complex was also modeled here as a precursor to the organophosphate calculations.

UO2(CO3)22-

Uranyl dicarbonate is a common complex in typical groundwaters around pH 8 [70]. UO2(OH2)2(CO3)22-•28(H2O) was modeled in two configurations with respect to the carbonate groups: trans and cis. Both optimized structures produced interatomic U-O (2.43 Å) and U---C (2.88 Å) distances in close agreement with experiment (Table 1). A similar computational methodology was used in Majumdar et al. [72] to calculate structures of the UO2(CO3) complex, but these results predicted U-O bond lengths of ≈2.3 Å. Majumdar and Balasubramanian [73] used B3LYP, MP2 and coupled cluster doubles (CCD) methods with COSMO solvation [74] to model UO2(CO3)22- (with and without Li+ and Na+) and predicted U-O(C) bonds of 2.34 Å without Li+ or Na+. This distance increased to 2.35-2.37 Å with the charge-balancing alkalis (No U---C distance was reported in this paper). The difference of approximately 0.1 Å between these models and observation is the result of neglecting the H-bonding due to explicit solvation of the carbonate group. Lack of H-bonding causes a stronger attraction of the [CO3]2- group to the uranyl and will overestimate the complexation energy as well as underestimate the bond lengths. These differences highlight the need to include at least one solvation sphere when modeling aqueous anions to account for the strong H-bonding that occurs to these species [75].

The trans uranyl-dicarbonate species was predicted to be -9 kJ/mol lower in Gibbs free energy (Table 1). The trans configuration should be more stable because it reduces the electrostatic repulsion between the negative carbonate groups, but the free energy difference is not large compared when one considers thermal and entropic effects. Because including diffuse functions in the basis set is relatively more important for anionic species, potential energy minimizations were also run with the B3LYP/6-31+G(d,p) method on the H, C and O atoms for these two configurations. In this case, the calculated ΔG = +23 kJ/mol in favor of the cis configuration. In both cases, the calculated ΔG is within the range of error in the computational methodology employed, so we expect that both configurations could exist in solution.

UO2(CO3)34-

Fitting of EXAFS spectra has resulted in U=Oax, U-Oeq, U-Odistal, and U-C distances of approximately 1.8, 2.4, 4.26 and 2.9 Å, respectively, for UO2(CO3)34-(aq) [43,49,50,76]. In comparison, the model UO2(CO3)34-•28(H2O) (Fig. 2a) predicted distances of 1.82, 2.50, 4.20, and 2.95 Å, respectively (Table 1; Addition of a diffuse function, e.g., the 6-31+G(d,p) basis set resulted in distances of 1.82, 2.52, 4.23 and 2.97 Å, respectively – relatively insignificant changes. The average H-bond distance was the same – 1.87Å – in both energy minimizations.). Even U-Odistal distances (i.e., the O atoms on the CO32- groups not directly bonded to the U atom) are accurate to within 0.07 Å (Table 1). Gagliardi et al. [77] obtained similar values using second-order perturbation theory (MBPT2) calculations including a reaction field Hamiltonian with a spherical cavity to account for solvation effects.

Figure 2

Model structures of the aqueous species (a) UO. (a) Strong H-bonding to the carbonate groups weakens the U-carbonate bonding by approximately 0.1 Å to bring the calculated value in better agreement with observation (Table 1; [42]). H-bonds to the O atoms of the UO22+ group are relatively weak. Reduction of the U atom to U(V) causes a slight twisting of the carbonate ligands as proposed by Docrat et al. [43]. H-bonding to the uranyl O atoms becomes relatively stronger compared to the analogous U(VI) complex. (c) Addition of Ca2+ ions to charge-balance this model aqueous species results in a configuration close to the observed crystal structure of the mineral liebigite [45].

A δ13C of approximately 165.8 ppm has been observed for UO2(CO3)34- at pH 8 [78] and 168.2 ppm at pH 12 [79]. Differences in ionic strength between these two experiments could also contribute to this change in observed δ13C. Also, a value of 170.5 ppm was observed by Brucher et al. [80] at pH 9. The calculated value is 168.5 ppm, consistent with observation considering the accuracy of the calculations is a few ppm.

The observed δ17O values for axial O and complexed carbonate in UO2(CO3)34- are 1098 and 215 ppm, respectively [81]. Average model values in this study are 1102 and 314 ppm, respectively. The calculated U=Oax are in excellent agreement with observation, but the complexed [CO3]2- do not appear to be. However, there are two groups of 17O chemical shieldings for complexed carbonate in these calculations – positive and negative. The O atoms of the carbonate groups that are not bonded to U all have positive 17O chemical shifts (31, 36 and 60 ppm); but the O atoms involved in the U-O-C linkages can have either negative or positive chemical shifts (-69, 63, 10, 35, -40 and 33 ppm). If the positive and negative values are separated, the positive values give an average δ17O of 252 ppm and the negative values an average δ17O of 347 ppm. The former is in somewhat closer agreement with the observed 215 ppm value, and the latter is off the scale in the published spectrum [81]. One can conclude that the computational methodology and solvation model are reasonably accurate for predicting δ13C values in this complex, but that solvation effects complicate prediction of the δ17O values.

Anionic complexes such as this can be more difficult to model due to the strong H-bonding that occurs to the O atoms of the anion [75]. The increased solubility of these U-carbonate complexes is due to this strong H-bonding to the carbonate groups. In our model, both the Oeq and the Odistal atoms form H-bonds to the surrounding H2O molecules. Furthermore, the increased electron density of this complex over the U(VI)aq species led to the formation of shorter H-bonds (≈1.86 Å) to the U-Oax than were predicted for the uncomplexed UO22+ ion. These model values are slightly less accurate compared to experiment than those in de Jong et al. [6] based on small core Stuttgart relativistic effective core potentials and local density approximation DFT calculations with diffuse functions added to the basis sets, but the results in this study are still considered in reasonable agreement with observation.

UO2(CO3)3 Ca 2

The addition of Ca2+ to aqueous solutions of uranyl-carbonates has a dramatic effect on increasing their stability and decreasing reactivity [49,82-84]. To model this effect, two Ca2+ ions were added to the UO2(CO3)34-•28(H2O) model to form Ca2UO2(CO3)3•28(H2O) (Fig. 2c). Two structures were energy minimized that varied in their Ca2+ positions: one with both Ca2+ ions approximately in the plane of the carbonate groups (Fig. 2c) and another with both Ca2+ out of the equatorial plane. Both had similar structures with regard to the uranyl-carbonate interatomic distances (Table 1). One exception was the U---Ca distances where the complex with two Ca2+ ions in-plane predicted distances of 4.05 Å, much closer to the experimental 3.95 Å compared to the 3.68 Å of the out-of-plane complex. The calculated stability of the in-plane structure is also -120 kJ/mol lower in Gibbs free energy than the latter (Table 1). This suggests a fairly strong regular structure for this ternary complex which will influence how it can react with mineral and bacterial surfaces [71,84].

UO2 [CO3]35-

Docrat et al. [43] reported the structure of aqueous UO2 [CO3]35- (i.e., the U(V) version of the uranyl triscarbonate species; Fig. 2b) based on EXAFS, and the interatomic distances are listed in Table 1. Starting with the model optimized structure of UO2 [CO3]34-, an electron was added and energy minimization performed. The calculated interatomic distances are reasonably close to the experimental values (Table 1; [43,44,50]). The model results are in closer agreement with the values of Docrat et al. [44] and Ikeda et al. [50] who found shorter U=Oax bonds and longer U-Oeq, U---C and U---Odistal bonds than Renshaw et al. [44], but the discrepancies are still generally larger between observation and model than between the two experimental studies.

The observed δ13C NMR chemical shift for UO2 [CO3]35- is 106.7 ppm [79] compared to the 113.6 ppm average value calculated in our model which is in reasonable agreement for a paramagnetic complex considering that pseudo-contact shifts were not considered. These differences with experiment indicate that the U(V)-triscarbonato complex is not modeled as accurately as the corresponding U(VI) complex, but calculating NMR shifts can be problematic in U-containing systems where relativisitic effects may be significant [85], especially for an open shell system such as UO2 [CO3]35-. The modeling does reproduce the large observed decrease in δ13C from UO2 [CO3]34- to UO2 [CO3]35- (observed Δδ13C = -62 ppm [79] verus calculated Δδ13C = -55 ppm this study).

Similar changes in interatomic distances and H-bonding were predicted for the U(V)-carbonate complex as were calculated for the aqueous U(V) species compared to U(VI). This result is consistent with the suggestion of Tsushima et al. [65] and Ikeda et al. [50] that H-bonding to axial oxygens should increase as the U atom is reduced from U(VI) to U(V).

Aqueous U(V) rapidly disproportionates into U(IV) and U(VI), but the UO2 [CO3]35- complex can be stable for 2 h [43]. Recently, Ikeda et al. [50] have shown that the redox between U(VI)- and U(V)-triscarbonato species occurs "quasi-reversibly" and that U(V)-carbonate species can be stable for up to two weeks in a sealed cuvette. This change in electrochemical behavior with complexation has been suggested to be due to changes in the conformation of the uranyl-carbonate complexes with reduction from U(VI) to U(V) [43]. However, the similarity in measured interatomic distances via EXAFS [43,44] makes it difficult to explain the change in disproportionation kinetics. One possibility is that the carbonate group twists and the Oeq atoms are displaced from the U-Oeq plane defined in the UO2 [CO3]34- complex [43].

The model results indicate that such a twist does occur because the U-O-C-O torsions change from 163, 170, and 175° for the UO2(CO3)34- complex to 160, 170, and 173° for the UO2 [CO3]35- complex. This twist is not dramatic and the UO2(CO3)34- complex does not have all six Oeq atoms in the same plane to begin with. However, the static nature of these energy minimizations does not allow one to estimate the dynamic distortion that will occur in each type of complex at finite temperatures. Quantum mechanical molecular dynamics simulations of these complexes would be useful for examining this phenomenon.

The calculated Gibbs free energy differences between U(VI)aq/U(V)aq and UO2 [CO3]34-/UO2 [CO3]35- are significant according to these calculations. The UO22+(OH2)6•18(H2O) model is +607 kJ/mol higher in energy than UO2+(OH2)6•18(H2O), and the UO2(CO3)34-•28(H2O) model is +479 kJ/mol higher in energy than UO2(CO3)35-•28(H2O) (Table 1). (Note that these are reduction half-reactions and do not include the energy of the electron, so these results do not imply that the reduced species are more stable overall. See [65] for methods and results on reduction energies between U(VI) and U(V).) The calculated reduction energies are not expected to be highly accurate at this level of theory and size of basis set, but the large relative change should be qualitatively correct. The direction of the change is unexpected because it would suggest that the U(V)aq species is relatively more stable with respect to re-oxidation and/or disproportionation compared to the carbonate complex. The lower calculated energy difference of the UO2(CO3)34-•28(H2O) complex compared to the uncomplexed U(VI)aq model relative to the U(V) version of each species is in direct contrast to the longer stability of the U(V)-carbonate complex. The observed metastability of the UO2(CO3)35-•28(H2O) complex is likely due to kinetic factors rather than thermodynamics. Surrounding the U(V) with carbonate anions may hinder the reaction between two U(V) atoms as electrostatic repulsion raises the activation energy barrier between anions. Modeling of the disproportionation reaction of complexed and uncomplexed U(V) species could help address this hypothesis.

Uranyl phosphate

The calculated structure of the uranyl-phosphate complex in Majumdar et al. [86] is close to the observed interatomic distances [46]. In this instance, the ligand has been protonated (i.e., HPO42-) in the previous theoretical study, and the U=Oax and U-O(P) distances match experiment to within 0.01 Å. The U-O(H2) distance is 0.08 Å larger than observed which may be the result of the model U existing in 7-fold coordination compared to 6-fold coordination extracted from the experimental data. One weak bond, such as the U-OH2 bonds, may break at room temperature to form fewer, but shorter, bonds compared to the 0 K energy minimization structure obtained computationally. Another discrepancy is that the U---P distance is 0.15 Å longer in the model of Majumdar et al. [86] compared to the value obtained by Locock et al. [46].

The protonation state of the phosphate ligand will vary with pH (pKa's of 2.1, 7.2, and 12.4), so for the purposes of this study, H2PO41-, HPO42- and PO43- ligands bonded to UO22+ were studied. Each type of uranyl-phosphate complex has been detected in solution depending on conditions ([87] and references therein). In terms of bond valence, bonding of O atoms in the phosphate ligand to U may approximate protonation of the O atom. Consequently, formation of a bidentate complex between HPO42- and UO22+ results in a phosphate environment similar to the H3PO4 species. H3PO4 only exists at acidic conditions (pKa1 = 2.1), so the bidentate UO2-HPO4 complex modeled in Majumdar et al. [86] may not be the most stable uranyl-phosphate species under circumneutral conditions if O bonding to U approximates protonation. Brendler et al. [88] did interpret their potentiometric and time-resolved laser-induced fluorescence spectroscopy (TRLFS) data consistent with formation of [UO2-H2PO4]+ and UO2-(H2PO4)2 under a pH range of 2 to 5. If a HPO4-UO2 complex does form, it may be a monodentate species, so this possibility was also investigated here.

Table 1 shows that both the UO2-PO4 and UO2-HPO4 model bidentate complexes result in reasonable agreement with the experimental interatomic distances reported in Locock et al. [46]. Two attempts at optimizing a monodentate UO2-HPO4 complex were made using initial configurations that should favor the monodentate complex (i.e., nearly linear U-O-P angles and strong H-bonding to the phosphate O atoms not bonded to the U). However, both resulted in conversion to a bidentate configuration (Fig. 3a). The approximate potential energy difference between the mono- and bidentate configurations was greater than -1000 kJ/mol. Even considering the partial optimization of the monodentate structure and all the approximations and errors inherent in these model calculations, this difference is extremely large and suggests that the monodentate configuration is not stable for UO2-HPO4. The possibility exists for a monodenate UO2-H2PO4 complex as discussed below.

Figure 3

Uranyl-phosphate aqueous models for (a) UO.

The possibility of a 1:2 uranyl phosphate complex also exists. Two HPO42- ions could bond to uranyl to form a monodentate complex, so this structure was used as an initial configuration. As in the 1:1 complex, one of the HPO4 groups re-oriented itself to form a bidentate configuration during the energy minimization (Fig. 3b). The second HPO4 remained in a monodentate configuration throughout the energy minimization. The monodentate configuration had a shorter U-O(P) distance (2.22 Å) and a longer U---P distance (3.78 Å) than experimentally observed (Table 2; [46]). The monodentate configuration is maintained with a nearly linear U-O-P angle of 167° and by H-bonding of H2O molecules to the free O atoms of the HPO4 group. Based on the results discussed above and the fact that the U-O(P) and U---P distances do not agree with observation, it is likely that this monodentate configuration is metastable.

Table 2

Uranium organic complex potential energies (Hartrees/molecule), interatomic distances (in Å) and calculated 13C NMR chemical shifts (in ppm).

Models	G	CN	U=O	U-O(H2)	U-O(C)	δ13C
Oxalate
Expt (This work)	-----	---	-----	-----	-----	174.5
HOx-•8(H2O)	-989.3358	---	-----	-----	-----	176.2
Ox2-•8(H2O)	-988.7157	---	-----	-----	-----	179.3
Uranyl-oxalate
Expt (This work)	-----	---	-----	-----	-----	169.0
[UO2(OH2)4]2+-Ox	-885.0151	7	1.79	2.61	2.28	164.7
[UO2(OH2)4]2+-Ox•11(H2O)	-1725.8906	7	1.79	2.481	2.41	172.0
[UO2(OH2)4]2+-Ox•30(H2O)	-3178.3494	7	1.79	2.47	2.40	167.4
Catechol
Expt [102]	DMSO	---	-----	-----	-----	116.1, 119.8, 145.6
	Solid	---	-----	-----	-----	115.7, 121.7, 142.7
H2Cat	-382.6927	---	-----	-----	-----	115.5, 121.1, 144.3
Uranyl-catechol
Expt (This work)	DMSO					115.7, 119.3, 145.3
	Solid					115.7, 122.3, 143.4
[UO2(OH2)4]2+-H2Cat	-890.0170	8	1.76	2.62	2.60	116.1, 132.4, 137.7
[UO2(OH2)4]2+-H2Cat (OS)	-3029.9112	7	1.79	2.45	≈5.9	115.7, 126.1
						123.2, 128.6
						138.5, 147.4
[UO2(OH2)4]2+-H2Cat•28(H2O) (Singlet)
	-3029.8947	7	1.79	2.45	2.48	121.8, 128.3, 141.6
[UO(OH)(OH2)4]2+-HCat•28(H2O) (Triplet)
	-3029.9386	7	1.912	2.50	2.42	138.3, 179.1,113.2, 114.2,127.6, 128.3,
[U(OH)2(OH2)4]2+-Cat•28(H2O) (Quintet)
	-3029.9477	7	2.173	2.43	2.60	192.0, 195.3
						129.4, 151.2
						123.0, 137.1
[U(OH)2(OH2)4]2+-Cat•28(H2O) (Singlet)4
	-3030.6937	--	-----	-----	-----	160.5, 181.6,
						116.7, 144.5,
						149.6, 159.4
[U(OH)2(OH2)4]2+-Quinone•28(H2O) (OS)
	-3029.9359	7	2.17	2.46	≈6	114.9, 194.3
						137.7, 166.3
						159.8, 202.6

δ13C relative to TMS HF/6-31G(d,p)//B3LYP/6-31G(d,p) chemical shielding of 201.4 and δ17O relative to water B3LYP/6-311+G(d,p)//B3LYP/6-31G(d,p) chemical shielding of 293 ppm. δ13C shifts listed in order as the two phenolic C atoms, the two C atoms α with respect to the phenol groups, and the two C atoms β with respect to the phenol groups. (OS = outer-sphere)

1 = One U-OH2 deprotonated to form U-OH and a H3O+ in the solvation sphere

2 = One U-O at 1.85 Å and one U-OH at 1.96 Å

3 = One U-OH at 2.14 Å and one U-OH at 2.21 Å

4 = aqueous Gibbs free energy in the quintet structure without ZPE correction

The relative stability of this monodentate HPO42- ligand was checked by running an energy minimization of the same model with both HPO42- groups starting in a bidentate configuration. The model with both HPO42- groups bidentate was -101 kJ/mol lower in Gibbs free energy than the mixed monodentate/bidentate model (Table 1). Again, these results indicate that the monodentate configuration is less stable for HPO42- ligands than the bidentate configuration.

A stable monodentate configuration for the 1:1 UO2-H2PO4 model complex was found. The U---P interatomic distance in this configuration was 3.64 Å in poor agreement with Locock et al. [46] (Table 1) but in excellent agreement with a value for U---P for uranyl adsorbed to phosphoryl groups on bacterial surfaces at low pH (3.64 ± 0.01 Å; [47]). The Kelly et al. [47] value was interpreted as a monodentate phosphoryl bonded to uranyl. The Gibbs free energy of the monodentate configuration was +9 kJ/mol higher than the same model system in a bidentate configuration (UO2H2PO4 (bi); Table 1). The above results all indicate that the mechanism of ligand binding to uranyl ions is sensitive to pH and the number of ligands involved. These factors need to be considered as spectroscopic data collected under specific conditions is applied to surface complexation modeling and to field studies.

IR and Raman

Quantum mechanical calculations can be useful in helping to identify speciation and protonation changes in conjunction with IR and Raman spectroscopies as well. By predicting vibrational frequencies and modes, assignment of observed spectra to specific complexes can be done with greater confidence.

Allen et al. [78] observed a ν1(O=U=O) peak at 813 cm-1 in Raman spectra of aqueous UO2(CO3)34-; this vibrational mode is predicted to be at 871 cm-1 in these model calculations (average of two modes at 865 and 877 cm-1). Model frequencies may be approximately 4% too high due to neglecting anharmonicity [30]), so this agreement is reasonable. However, recent work has suggested that the overestimate of the vibrational frequencies for the uranyl ion is more likely due to basis set and electron correlation effects rather than anharmonicity [89]. Bargar et al. [42] measured ATR FTIR spectra of aqueous UO2(CO3)34- and found peaks near 1360 and 1510 cm-1 associated with the carbonate vibrational modes (symmetric and asymmetric ν3 modes, respectively). The model UO2(CO3)34-•28(H2O) used in this study resulted in frequencies of 1302, 1372, 1427 and 1545 cm-1 for these modes with the highest IR intensity in this range (Fig. 1c).

At first, the disagreement would seem to be significant, but when one considers that the observed bands have full-widths at half-maximum of approximately 50 cm-1, then the possibility can be considered that these observed bands are mixtures of the calculated bands. The IR intensity weighted-averages of the two sets of frequencies result in calculated values 1338 and 1503 cm-1 that are as close as can be expected to the observed frequencies for modeling aqueous phase complexes at this level of theory [90]. Without including H2O molecules of solvation around the complex, de Jong et al. [6] calculated sets of peaks at 1356 + 1383 and 1508 + 1533 cm-1. Hence, either solvation is not critical for calculating the UO2(CO3)34- frequencies or there are compensating effects in the DFT approach used by de Jong et al. (2005). Previous work has shown that strong H-bonding can dramatically affect carbonate vibrational frequencies [75], so the reason for the excellent agreement in this case with the gas-phase calculations should be investigated further by performing a systematic study of explicit and implicit solvation on this model.

Although neither IR nor Raman spectra are available for uranyl-phosphate aqueous solutions, vibrational frequencies calculated (Fig. 1d) for the UO2-PO4- complex at 818 (U=Os), 888 (U=Oas), 1021, 1052, 1097 and 1104 cm-1 (PO4 stretches) correspond well with observed Raman peaks in the mineral threadgoldite (Al [(UO2)2(PO4)2](OH)(H2O)8) at 827, 913, 1019 to 1026, 1057 and 1107 cm-1 [91]. Calculated peaks around 900 cm-1 were not observed in the Raman spectrum of this mineral nor did the calculations predict observed peaks between 952 and 976 cm-1 [91]. Considering the fact that a model aqueous complex is being compared to a mineral, the level of agreement is remarkable and suggests that the bonding environment of the uranyl-phosphate complex is being modeled accurately.

UO2-Organic

Complexation of UO22+ with oxalate was modeled because oxalate is a simple organic ligand that forms strong complexes with cations. In addition, oxalate is a common organic acid in nature. Alliot et al. [92] studied the role that uranyl-oxalate complexation plays in the adsorption of uranium to alumina and found that making the charge of the complex neutral or negative would cause uranium to desorb due to higher solubility of the aqueous complex. Previous quantum mechanical calculations have been performed on uranyl-oxalate complexes by Vallet et al. [93] using Hartree-Fock gas-phase energy minimizations and single point energies based on MP2 calculations in the conductor-like polarized continuum model (CPCM) for solvent effects. Our structural results including explication solvation by H2O molecules will be compared to these earlier calculations.

UO2-Oxalate with 4, 11, and 30 H2O

Structure

Selected structural parameters for the uranyl-oxalate complex as a function of the number of H2O molecules of solvation are listed in Table 2. In each case, the UO22+ cation was initially bonded to four H2O molecules. The coordination number does not change as a function of H2O's of solvation which is significant because one might expect the non-solvated complex to be able to form stronger bonds to the H2O molecules. However, in all three cases, one of the four H2O's was kicked out of the coordination sphere and ended up H-bonding to the complex (Fig. 4). U=O bond lengths also do not change significantly as a function of solvation. This result is an indication of the weak H-bonding that occurs to the doubly bonded O atoms in U(VI) species.

Figure 4

Model structure of aqueous uranyl-oxalate complex, UO.

The U-Oeq bonds change noticeably upon introduction of H2O molecules surrounding the complex. The U-O(H2) bonds shorten by 0.14 Å and the U-O(C) bonds lengthen by 0.12 Å from [UO2(OH2)4]2+-Ox to [UO2(OH2)4]2+-Ox•30(H2O). This effect is similar to that seen in Majumdar et al. [72] where U-O bonds to carbonate were decreased in length (and presumably increased in strength) by the neglect of solvation effects. Such changes are likely to have important effects when calculating the stability constants of uranyl complexes. The gas-phase energy minimizations of Vallet et al. [93] predicted U-O(H2) and U-O(C) bonds of 2.59 and 2.30 Å which were similar to the values obtained in this study without explicit solvation.

IR, Raman and NMR

The observed Raman spectrum of a 1:2 uranyl and oxalic acid solution at pH 3 and a synthetic Raman spectrum of the [UO2(OH2)4]2+-Ox•30(H2O) complex are shown in Figs. 5a and 5b. Observed Raman frequencies are found at approximately 850, 1050, 1120, 1144, 1270, 1300, 1350 and 1450 cm-1 ([94] and this study). The peaks at 1144, 1270, 1350 and 1450 cm-1 are all observed in uranyl-free oxalic acid solutions (Note that the presence of free oxalate indicates that the 1:1 uranyl-oxalate complex is most likely dominant with the equilibria UO2 + oxalate ⇔ UO2(Ox) ⇔ UO2(Ox)2.). The 850 cm-1 is assigned to U=O stretching (see above). Consequently, the 1050, 1120, and 1300 cm-1 modes may be due to uranyl-oxalate complexes. The most intense model Raman peaks associated with uranyl-oxalate vibrational modes in [UO2(OH2)4]2+-Ox•30(H2O) are found at 1235 (C-Os) and 1380 (C-Oas); hence, there is significant disagreement between the model predictions and observed Raman frequencies. One could conclude that the computational methodology is not accurately reproducing oxalate vibrations, except that the Ox2-•8(H2O) and HOx-•8(H2O) models predict Raman peaks at approximately 1285 and 1440 cm-1 and 1250, 1390 and 1460 cm-1, respectively which correspond to the free oxalate Raman peaks. Thus, we conclude that oxalate is not complexing with the uranyl monomer and must be complexing with a hydrolyzed uranyl oligomer.

Figure 5

(a) Observed Raman spectrum of uranium-oxalate 1:1 solution at pH 10. (b) Calculated Raman spectrum of UO2(OH2)4(C2O4)•30(H2O) (UO2Ox-W30) shows a sharp peak near 1800 cm-1 indicative of C=O bonds that would be present if oxalate forms a bidentate complex with uranyl.

A model IR spectrum is shown in Fig. 5c. The strongest IR adsorbing modes group into three regions in the 500 to 2000 cm-1 part of the spectrum shown (higher frequencies are neglected because they are dominated by motions of H2O): 700 to 900, 1200 to 1400, and 1600 to 1800 cm-1. The peak at 1797 cm-1 is the C=O stretch. The group of peaks just below this sharp peak in Fig. 5a is due to H-O-H bending motions of the solvating H2O molecules. The 1233 and 1380 cm-1 peaks are modes involving both C-O(U) and C-C stretches of the oxalate. The U=Os stretch occurs at 839 cm-1 and the U=Oas stretch at 909 cm-1. The other peaks in this region are due to motions of the H2O molecules. The observed O=U=O stretch [94] falls between the predicted vibrations at 839 and 909 cm-1. From the observed O=U=O stretch, Tsushima et al. [94] calculated a 1.74 Å U=Oax bond length based on Badger's rule (this study U=Oax = 1.79 Å; Table 2). Quantum mechanical calculations in Tsushima et al. [94] predicted a 1.69 Å U=Oax bond, but this underestimate is likely due to the fact that solvating H2O molecules were not included in this early study. The 1144 and 1270 cm-1 peaks compare favorably with the calculated peaks at 1184 and 1275 cm-1 (neglecting anharmonicity in the calculations tends to increase the predicted frequencies). The peak at 1450 cm-1 in the experimental spectrum has no equivalent in the model uranyl-oxalate complex and is due to uncomplexed oxalate. A peak at 1483 cm-1 has been observed for aqueous Na-oxalate by attenuated total reflectance FTIR [95].

IR spectra have been collected by Giesting et al. [96] on a crystal of (UO2)2C2O4(OH)2(H2O)2 but the uranyl-oxalate configuration is not similar to the aqueous model because the oxalate ligand is bonded to two uranyl groups. Thus, the vibrational spectra will be completely different from the case where two C-O(U) and two C=O bonds exist in the oxalate. A literature search revealed no other IR or Raman spectra of uranyl-oxalate complexes, so these calculations stand as predictions in the event this type of study is conducted in the future.

Values for δ13C of the uranyl-oxalate complexes are listed in Table 2 and the NMR spectra shown in Figure 6. The HF/6-31G(d,p) method was used to calculate chemical shifts relative to TMS and the experimental values were measured relative to trimethylsylyl propionate, but the difference in chemical shielding between TMS and TMP is only -0.08 ppm [97]. Hartree-Fock calculations have proven to be more accurate for reproducing δ13C values in other calculations of organic complexes [98], so this method was used instead of a DFT method such as B3LYP. (B3LYP was tested but results did not match experiment as closely for test cases.) The observed δ13C of oxalic acid in solution is approximately 174.5 ppm; model calculations on HOx-•8(H2O) and Ox2-•8(H2O) result in predicted values of 176.2 and 179.3 ppm, respectively (Table 2). Thus, the liquid-state NMR experiment appears to be observing the singly deprotonated oxalic acid species, a time-averaged value of the singly and doubly deprotonated oxalic acid, or the method is accurate to only 5 ppm. The experimental uranyl-oxalate δ13C (169 ppm) and calculated (167 ppm) values are in agreement within computational error for the [UO2(OH2)4]2+-Ox•30(H2O) model. Since the measurement was at pH 3, then the spectrum corresponds to H2Ox (pKa = 1.27) and HOx- (pKa = 4.28) in rapid exchange equilibrium (time-average). Having an observed chemical shift that is more shielded than calculated for HOx- is consistent with this.

Figure 6

NMR spectra of oxalic acid and uranium-oxalate aqueous solutions.

Solvation appears to be important for obtaining more accurate chemical shifts because the models with no or few H2O molecules are less accurate using the same basis set and method. Note that the model value is the average of two distinct values for the oxalate C atoms: 171 and 164 ppm. Reproduction of the observed vibrational frequencies and 13C chemical shifts is evidence that the oxalate binds to UO22+ in the bidentate fashion as previously described by Vallet et al. [93].

UO2- Catecholate

Uranyl has been shown to interact with catechol in aqueous solutions. Catechol is an organic ligand generated by plants [99] and an analog for common functional groups within humic substances [100]. Consequently, modeling uranyl-catechol complexation is important for environments with significant amounts of natural organic matter present. Additionally, Eng [101] has shown that U(VI) can oxidize catechol and form organic polymers. Although the U(VI) was subsequently re-oxidized in these aerobic experiments, one can presume that the U(VI) is reduced to U(V) or U(IV) in the process of oxidizing the catechol. Calculations on the structures and 13C NMR chemical shifts of UO22+-catechol complexes are presented here and results on the uranyl-catechol redox reaction are discussed.

Solvated with 4 and 28 H2O

The catechol ligand does not deprotonate in aqueous solution with UO22+ as evidenced by NMR spectroscopy (Table 2; Fig. 7). Signal from both phenol protons on the catechol are present as indicted by the approximately 1:2 ratio of intensity relative to the C-H (Fig. 7) which is suggestive of an outer-sphere uranyl-catechol complex. This interpretation is supported by observation of 13C chemical shifts for freeze-dried uranyl-catechol solutions that are nearly identical (within 0.7 ppm) to solid catechol. Outer-sphere complexation leads to the formation of significantly weaker bonds between the catechol and uranyl compared to the oxalate ligand, and the U-O distances are approximately equal between the U-O(H2) and U-O(C) bonds in the UO22+-catechol complexes (Table 2). Notably, addition of 28 H2O molecules around the bare [UO2(OH2)4]2+-catechol complexes decreases both of the equatorial bond types by approximately 0.15 Å. This is opposite to the effect predicted for carbonate and oxalate ligands which lose electron density to surrounding H2O molecules, and we have no explanation for this bond shortening at this time. The weakest of the U-O(H2) bonds is also broken with the addition of H2O's of solvation to lower the CN to 7 in the solvated case (Table 2).

Figure 7

NMR spectra of uranyl-catechol complexes. (a) 13C{1H} CP/MAS spectrum of freeze-dried pH 3.5 solution of 1:1 uranyl:catechol. (b) and (c) solution-state 13C and 1H spectra, respectively, of sample from (a) re-dissolved in DMSO-d6.

Catechol has three pairs of equivalent carbon atoms with distinct δ13C values (Table 2). For uncomplexed catechol in DMSO solvent, these are found at 116.1 (C atoms next to the COH groups), 119.8 (for the C atoms opposite the COH groups), and 145.6 ppm (for the phenolic C atoms), respectively, which are very similar to values of 115.7, 121.7, and 142.7 ppm for the neat solid [102]. HF/6-31G(d,p) NMR calculations on the B3LYP/6-31G(d,p) optimized structures result in predicted values of 115.5, 121.1, and 144.3 ppm for each of these peaks (Table 2), so this methodology is expected to reproduce δ13C values to approximately 2 ppm. Calculated values for the inner-sphere [UO2(OH2)4]2+-H2Cat•28(H2O) complex are in disagreement with the observed 13C chemical shifts by 4 to 6 ppm (Table 2). Thus, this model with direct bonding of the phenol O atoms to the uranyl cation is unlikely to exist in solution. The other inner-sphere models, [UO(OH)(OH2)4]2+-HCat•28(H2O) and [U(OH)2(OH2)4]2+-Cat•28(H2O), result in even larger discrepancies from experiment predicting six separate chemical shifts ranging from 113 to 342 ppm, so another type of uranyl-catechol complex must give rise to the observed 13C NMR spectra that show δ13C values shifted by approximately 0.5 ppm from the catechol-only solution (Fig. 7).

Raman spectra were collected on uranium-catechol aqueous solutions as well. The Raman spectra are complicated by the oxidation and polymerization of catechol during these experiments and models of the initial steps of this process will be discussed below. However, at pH 2 before oxidation and polymerization occurs, observed Raman frequencies of uranyl-catechol solutions are 580, 774, 874 (uranyl), 1040, 1050 (nitrate), 1160, 1274 (nitrate), 1355, 1468, 1497 and 1604 cm-1 (Fig. 8). Except for the nitrate peaks and peaks at 1355 to 1497 cm-1 that appear to be associated with polymerization of the catechol, the model Raman spectrum (Fig. 8) has corresponding calculated values at 595 (ring), 788 (CH wag), 1070 (CCH bend), 1187 (CCH bend), 1290 + 1307 (COH + CCH bends), 1639 cm-1 (C-C stretch).

Figure 8

Calculated (a) IR and (b) Raman spectra of UO.

In a search for the structure of the uranyl-catechol complex, the structure of a model uranyl-catechol outer-sphere complex was also calculated. The energy minimized structure is shown in Fig. 9d. Agreement between model and observed δ13C values is much improved over the inner-sphere and quinone models (Fig. 9e). Although significant errors remain for three of the δ13C values, the other three C atoms have calculated chemical shifts that are reasonably close to observation (Table 2). Furthermore, the stability of this species is calculated to be approximately -80 kJ/mol more favorable in potential energy (Table 2) than the inner-sphere [UO2(OH2)4]2+-H2Cat•28(H2O) complex shown in Fig. 9a. Apparently, the uranyl bonds to the phenolic groups are weaker than the UO22+-OH2 bonds that are broken to form them.

Figure 9

(a) [UO. Note the single bond to catechol in (a) and the H transfer to the axial uranyl O atoms in (b) concomitant with the changing electronic state and uranyl reduction (Table 3). The most stable state calculated is (e) where the uranyl has been reduced and a quinone has generated. This result is consistent with the experimental observations of [100] presuming that the U(IV) produced in this reaction was re-oxidized by O2 in these aerobic experiments.

One can conclude from the NMR and energy results that the [UO2(OH2)4]2+-H2Cat (OS) model is the closest representation of those studied here of the actual uranyl-catechol interactions in low pH experimental solutions. Modeling the structure of an outer-sphere complex is more challenging than inner-sphere complexes that are dominated by covalent bonding because the outer-sphere complex is likely to be much more flexible and experience a wider range of configurations. Thus, conformation searching for the uranyl-catechol outer-sphere complex would be useful.

U(VI) reduction and catechol oxidation

As mentioned above, oxidation of catechol has been observed in conjunction with uranyl-catechol complex formation [101]. To begin modeling this redox reaction, the [UO2(OH2)4]2+-H2Cat•28(H2O) complex that had been energy minimized in the singlet state (Fig. 9a) was also energy minimized in the triplet (Fig. 9b; [UO(OH)(OH2)4]2+-HCat•28(H2O)) and quintet states (Fig. 9c; [U(OH)2(OH2)4]2+-Cat•28(H2O)). These electronic states allow the possibility of one- and two-electron transfers from the catechol to the U(VI) atom. The Gibbs free energies of the triplet and quintet states were -72 and -96 kJ/mol lower (i.e., more favorable) than the [UO2(OH2)4]2+-H2Cat (OS) structure using our computational methodology (Table 2). The lower energy of the quintet state is consistent with the prediction above that the U(IV)aq species is in the triplet rather than singlet state.

The computed charge on the U atom changes from 3.42 to 3.29 to 3.19 using natural bond orbital analysis [103-107] as the multiplicity of the complex changed from singlet to triplet to quintuplet. The charge on the UO22+ group was predicted to decrease from +1.10 to +0.65 to +0.31 (Table 3) reflecting a reduction of the uranyl cation (Note that the NPA charges on the uranium atom only change by 0.1 electron.). The atomic spin densities on U also increased from 0 to 1.27 to 2.10 in the singlet, triplet and quintet state, respectively. The remaining spin density in each case was located on the C and O atoms of the original catechol. This result is consistent with the observation of polymerization in these solutions [101] via a radical mechanism.

Table 3

Uranium-catechol aqueous Gibbs free energies (Hartrees/molecule) and charges as calculated with the B3LYP/6-31G(d,p) basis set on H, C and O and the Stuttgart ECP60MWB on U using the natural bond orbital analysis program [103-107] with Gaussian 03 (Frisch et al., 2004).

Models	U Charge	UO22+Charge
[UO2(OH2)4]2+-H2Cat (OS)•28(H2O)
	+2.45	+1.11
Singlet
[UO2(OH2)4]2+-H2Catechol•28(H2O)	+2.46	+1.10
Triplet
[UO(OH)(OH2)4]2+-HCatechol•28(H2O)	+2.47	+0.65
Quintet
[U(OH)2(OH2)4]2+-Catechol•28(H2O)	+2.28	+0.13
[U(OH)2(OH2)4]2+-Quinone•28(H2O) (OS)	+2.52	+0.32

A modified l607 routine was used for the explicitly solvated models to include the U 6d electrons in the valence space for the natural population analysis [40]. The multiplicities (singlet, triplet and quintet) correspond to U formal oxidation states of +6, +5 and +4. The "U Charge" and "UO22+Charge" headings indicate the actual charge calculated using NBO. Note that the protonation state of the catechol (H2Catechol, HCatechol and Catechol) indicates a H+transfer to the uranyl group simultaneous with the electron transfer.

These calculations predict a thermodynamic driving force for the reduction reaction; but, because the U(VI)-catechol complex is stable enough to give rise to Raman and NMR spectra in low pH solutions, the activation energy barrier must be significant at room temperature as well consistent with the observation that oxidation/polymerization requires days to occur at room temperature [101]. Modeling the reaction pathway is beyond the scope of this paper, but the reaction mechanism is likely to involve a proton-coupled electron transfer because the stable species is [U(OH)2(OH2)4]2+-catecholate where the doubly-bonded O atoms in the uranyl cation have become protonated via transfer of H atoms from the phenol groups in the catechol (Fig. 9b). This pathway makes sense in light of the experimental evidence that a quinone-like species is produced when catechol is oxidized in this reaction [101]. Using the outer-sphere uranyl-catechol model (Fig. 9d) as the reactant and a U(IV)(OH)2(OH2)5-quinone•28(H2O) (Fig. 9e) as the product, the overall reaction ΔG is estimated to be approximately -65 kJ/mol (Table 2); however, our model results suggest that the triplet and quintet states in Table 2 are at a lower Gibbs free energy than the outer-sphere U(IV)(OH)2(OH2)5-quinone•28(H2O) complex.

UO2-Biological

Adsorption of uranium to bacterial surfaces is an important aqueous geochemical process (see [108]). This adsorption is generally thought to occur to functional groups located in the extracellular polymeric substances (EPS) located around many bacteria. Water may be present throughout EPS [109], so this section of the paper examines aqueous uranium complexation with EPS functional groups.

Previous work has concluded that the phosphoryl and carboxlate groups are the most important for binding uranyl [47,110,111]. Consequently, we have modeled two types of phosphoryl groups, one that mimics phosphoryl backbones in nucleic acids (OrgPO4; see also [112]) and another that represents phosphoryl group in glucosamine (GlcNPO4) which is thought to be an important component in EPS for metal adsorption [113]. The carboxylate group was modeled in the 2-Keto-3-deoxyoctanoate (KDO) component of EPS which has also been suggested as playing a role in uranyl adsorption to bacteria [114,115].

UO2-Phosphodiester

Previous workers have suggested that monodentate complexes of uranyl and bacterial EPS phosphoryl groups are the predominant mechanism of bonding (e.g., [47,71,111]). These calculations sought to compare the relative stabilities of monodentate and bidentate configurations. Model results in this study strongly indicate that the bidentate complex should be favored to the nucleic acid phosphoryl groups. Three initial configurations were subjected to energy minimizations using the same methodology described above for the inorganic and organic U-complexes. Two initial monodentate structures were used (e.g., Fig. 10a); in both cases, the energy minimization led to a bidentate structure (e.g., Fig. 10b) that was at least -520 kJ/mol lower in Gibbs free energy. Although relaxation of atoms in the model other than those involved in forming the second bond between the phosphoryl group and the uranyl account for some of this lowering in energy, it is clear that the bidentate configuration is lower in Gibbs free energy by much more that the expected errors in the calculations. (Note: The energy minimizations beginning with three different starting configurations resulted in final structures with energies within ± 25 kJ/mol of the average. This order of magnitude is used as an estimate of the uncertainty in the Gibbs free energies resulting from these energy minimization procedures.)

Figure 10

(a) Monodentate UO[47]showing monodentate phosphoryl bonding to U on bacterial surfaces at low pH.

Although the calculated energetics clearly favor the bidentate configuration, the calculated U---P distance in these complexes is approximately 3.3 Å (Table 4) whereas the observed value is 3.64 to 3.68 Å [111,116]. This discrepancy could be due to the fact that the Kelly et al. [47,111] experiments were conducted at a pH of 1.7 to 4.8, so the phosphoryl groups were probably protonated which prevented bidentate bonding. The model structures of UO2-OrgPO4 monodentate (Table 4) are also consistent with a shorter U-Oeq bond of 2.30 to 2.37 Å measured by Koban et al. [117] using EXAFS in solutions of pH 3.5 to 5.5. Francis et al. [116] performed EXAFS on samples reacted at pH 5 and obtained similar results as in Kelly et al. [111], but the expected pKa for phosphoryl groups on bacterial surfaces is approximately 7.2 [118], so experimental conditions may not have reached a state where the phosphoryl was predominantly deprotonated.

Table 4

Model aqueous Gibbs free energies (without ZPE corrections) and interatomic distances calculated for uranium model complexes with biological ligands. U-X stands for the shortest U to P or C distance in the model.

Models	Energy	CN	U=O	U-O(P)	U-O(H2)	U-X
Expt [47]	-----	8	1.77	2.33	2.45	3.64
UO2-OrgPO4•27(H2O) bi	-3601.0619	7	1.79	2.42*	2.57	3.18
UO2-OrgPO4•27(H2O)a
mono, initial	≈-3600.4331*	7	1.79	2.30	2.53	3.40
UO2-OrgPO4a•27(H2O)
bi, final	-3601.0751	8	1.77	2.57	2.61	3.24
UO2-OrgPO4b•27(H2O)
mono, initial	≈-3600.8669*	7	1.79	2.30	2.53	3.43
UO2-OrgPO4b•27(H2O)
bi, final	-3601.0662	7	1.77	2.47	2.66	3.27
UO2-OrgHPO4•27(H2O)
mono	-3601.5378	7	1.77	2.46	2.48	3.83
UO2-GlcNPO4•26(H2O) mono	-3652.6059	7	1.81	2.23	2.51	3.61
UO2-GlcNPO4•26(H2O) bi**	-3652.6170	7	1.80	2.43	2.52	3.18
Expt [47]	-----	8	1.77	2.33	2.45	2.89
UO2-KDO•26(H2O) mono	-3029.6108	6	1.76	2.39	2.34	3.44
UO2-KDO•26(H2O) bi	-3029.6144	7	1.75	2.42	2.51	2.89
UO2-KDO•26(H2O) OS	-3029.6099	7	1.75	-----	2.47	5.12

* – Gibbs free energies are estimated because no stable energy minimum was determined.

** – One U-OH2 has deprotonated to form a U-OH and a H3O+ in the solvation sphere.

In order to investigate the possibility that a protonated phosphoryl (i.e., mimicking low pH conditions) favors the monodentate uranyl complex, a uranyl-organophosphate complex (UO2-OrgHPO4; Fig. 10c) was investigated. A stable monodentate configuration was obtained when the phosphoryl group was protonated which resulted in a U---P distance of 3.8 Å (Table 4). The observed value is between the two calculated bidentate and monodentate values, so either the monodentate complex calculation is overestimating the U---P distance by 0.2 Å or there exists a mixture of mono- and bidentate uranyl-phosphoryl complexes on bacterial surfaces that result in an average value of 3.6 Å. Furthermore, in case of the protonated phosphoryl group, the calculated U-O(P) bond was 2.46 Å, significantly longer than the 2.30 to 2.37 Å [117] or 2.29 Å [116] observed. These discrepancies may be due to the role that multiple U-O-P linkages play in U binding to polyphosphate groups on bacterial surfaces [116,119]. These complexes were not modeled in this study and should be included in future work.

UO2-GlcNPO4

Another possible explanation of the discrepancies in interatomic distances between EXAFS and these calculations may be that uranyl binds to a phosphoryl group unlike the phosphodiester model. To examine this possibility, models of the glucosamine phosphoryl groups, GlcNPO4 (Fig. 11) were also studied which represents a phosphonate compound. Phosphonates have been shown to complex U(VI) effectively [113,120]. The UO2-GlcNPO4 model calculations also predicted that the bidentate configuration had a lower Gibbs free energy than the monodentate configuration but only by approximately -30 kJ/mol as compared with a difference of -520 kJ/mol for the phosphodiester (Table 4). In addition, a monodentate configuration was found in this case that represented a local potential energy minimum, so the calculated difference is less of an estimate in this case compared to the UO2-OrgPO4 models (Table 4).

Figure 11

(a) The monodentate UO[47]. These two facts suggest that protonation of the phosphoryl group at low pH stabilizes the monodentate configuration and that the bidentate configuration should be more stable at circumneutral pH.

Furthermore, the U---P distance calculated for monodentate UO2-GlcNPO4 is in excellent agreement with observation [47]. The U-O(P) bonds for the mondodentate and bidentate configurations are 2.23 and 2.43 Å, respectively. Thus, the U-O bond shortening observed by Koban et al. [117] is predicted, and the average of the two configurations results in a U-O(P) bond length close to that observed (i.e., 2.30 to 2.37 Å). Consequently, the conclusion that uranyl will bond to phosphoryl groups in a mondentate manner at low pH and convert to a bidentate complex at circumneutral pH is supported by these results because the model structures are similar to observation at lower pH, but theoretical potential energies of the bidentate configuration are significantly lower. These results also suggest that uranyl will favor glucosamine phosphoryl groups over phosphodiester groups. Micro-environmental conditions (such as lower dielectric constant or lower pH) within the EPS region could be responsible for stabilizing the monodentate over the bidentate configuration because these factors were not considered in the model calculations of this study.

UO2-KDO

As in the case of uranyl-phosphoryl binding, uranyl-carboxylate complexes may either form mono- or bidentate structures (Fig. 12). Previously, the bidentate configuration of uranyl with acetate has been modeled [6] and compared to experimental interatomic distances and frequencies. The relative potential energies of the mono- and bidentate UO2-KDO complexes and the model U---C distances are consistent with the interpretation that the bidentate configuration is favored (Table 4), but the Gibbs free energy difference between the two configurations is smaller compared to the phosphate complexes (i.e., < -10 kJ/mol). This relatively small difference suggests that the monodentate configuration could be favored under some circumstances, especially where U-hydroxide species are the stable aqueous species because OH- is a strong ligand that could compete with carboxylate groups.

Figure 12

UO[110], but the (b) bidentate configuration is calculated to be a lower energy state using the methodology discussed in this paper. (c) An outer-sphere configuration has a higher potential energy than either of the two inner-sphere model complexes.

The calculated UO2-KDO interatomic distances have significant discrepancies with the values reported in de Jong et al. [6] for the 1:1 UO2-acetate complex (U=O = 1.75 and 1.75, C-O = 1.29 and 1.27, U-O = 2.21 and 2.51, and U-C = 2.63 and 2.89 Å, for the de Jong et al. [6] study and the UO2-KDO values calculated in this study). The main differences are for the U-O bonds to the carboxylate group and the U---C distance. However, these differences are mainly due to the fact that the uranyl coordination sphere was not completed in the 1:1 complex modeled in de Jong et al. [6]. Their reported values for the U-O and U---C distances increase toward the values reported here in the 1:3 model uranyl-acetate complex and in the EXAFS data on solids [121].

Although the monodentate complex is predicted to be of slightly higher Gibbs free energy, it may be difficult to rule out the existence of this complex based on EXAFS results. The model interatomic distances in the monodentate complex are similar to those predicted for the bidentate complex except for the U---C distance (i.e., U=O = 1.76, C-O = 1.26, U-O = 2.34, and U---C = 3.45 Å). The latter is longer by approximately 0.5 Å, but this type of interatomic distance can be difficult to measure with EXAFS in materials such as bacterial EPS.

A third possible configuration is the outer-sphere pair formed by the uranyl cation and KDO (Fig. 12c). The outer-sphere model results in a higher potential energy than both the mono- and bidentate configurations by +12 kJ/mol. This value is well within computational accuracy and suggests that an equilbrium exists between inner- and outer-sphere species. Previously, it has been common to assume that inner-sphere species are in a lower energy state than outer-sphere, but recent work by Catalano et al. [122] on arsenate adsorption to hematite has demonstrated that this is not necessarily the case. Consideration of these outer-sphere species must be included in bacterial complexation modeling.

Vibrational spectroscopy could be more diagnostic in this case because the splitting of the carboxylate vibrational modes changes depending on whether the bonding is mono- or bidentate [123]. For the bidentate complex (Fig. 13a), the calculated C-Os mode is approximately 1500 cm-1 (there are four modes involving symmetric C-O stretching coupled with CH3 motion in the KDO; Fig. 13), and the C-Oas mode is at 1612 cm-1. The C-Os and C-Oas frequencies observed via ATR FTIR for bidentate uranyl acetate in solution are 1467 and 1527 cm-1, respectively. If the theoretical frequencies are scaled by 0.96 (the scale factor for organic molecules using B3LYP/6-31G(d) method; [124]), this results in errors of 27 and 20 cm-1, respectively, for these two modes. Considering the difficulties in model aqueous uranyl-organic complexes, and the fact that the model results are for UO2-KDO rather than the experimental uranyl acetate, this size of error is satisfactory.

Figure 13

Calculated IR spectra of UO.

The symmetric and asymmetric U=O vibrations are predicted to be approximately 900 (two frequencies at 895 and 906 cm-1) and 983 cm-1 (Fig. 13a). Both modes are relatively strong in this model, but the U=Oas vibration has about twice as much IR intensity. The calculated frequencies are close to those predicted by de Jong et al. (2005) for the 1:2 uranyl-acetate complex, but they represent a much smaller splitting between the U=Os and U=Oas modes than the 1:1 complex in de Jong et al. [6]. However, the theoretical vibrational frequencies are consistent with the interpretation of de Jong et al. [6] that adding more H2O molecules of solvation would lower the calculated splitting.

Measured U=O vibrational frequencies for the 1:1 complex in solution are 861 and 954 cm-1. [125]. Hence, the values calculated here are 3 to 5% higher than observed as expected, and the experimental and model splitting between the U=Os and U=Oas modes are 93 and 83 cm-1, respectively. One can conclude that the model is representing the vibrations of the uranyl-carboxylate complex fairly accurately, especially considering that the KDO molecule has been used here instead of acetate.

For the monodentate complex, the U=Os and U=Oas modes were calculated to occur at 894 and 982 cm-1 (Fig. 13b); and hence, these vibrations may not be helpful in distinguishing between mono- and bidentate complexation because they do not change frequency between the two model complexes. Consequently, tracking changes in the U=O modes may not be useful in distinguishing complex structures. On the other hand, the C-Os and C-Oas modes were computed at 1425 (three modes between 1420 and 1428 cm-1) and 1676 cm-1 (two modes at 1646 and 1686 cm-1). Estimating a scale factor or 0.96 based on the comparison to experiment above and assuming these calculations do as accurate a job on the monodentate complex vibrations, the C-Os and C-Oas stretches should occur near 1415 and 1610 cm-1. The calculated splitting between C-Os and C-Oas in the monodentate case is almost 200 cm-1 – more than 100 cm-1 greater than the observed splitting of these bands for uranyl acetate [125]. Hence, vibrational spectroscopy can be used as a relatively quick and inexpensive method for determining uranyl carboxylate bonding configurations. The success of time-resolved laser-induced fluorescence spectroscopy (TRLFS; e.g., [113] and references therein) also suggests that future quantum mechanical calculations focusing on predicting fluorescence spectra would be worthwhile.

Conclusion

The aqueous speciation calculations demonstrate the accuracy of the computational methods employed for predicting these structures because the model results are consistent with available experimental data. The U(V)aq species maintains the UO2 moiety, but the increased H-bonding to the axial O atoms is a precursor to the protonation of these atoms in the reduced U(IV)aq species. U(IV)aq is more stable in the triplet state over the singlet state by approximately 190 kJ/mol in these model calculations.

Uranyl-carbonate model complexes reproduce experimental EXAFS and NMR results provided explicit solvation (i.e., including H2O molecules) is included. The stability and structure of the ternary Ca-UO2-CO3 aqueous complex is also predicted. Kinetic hindrance of U(V) disproportionation by the presence of carbonate was not consistent with the hypothesis that distortions of the U-O-C-O torsion angle stabilizes U(V); instead, the slower disproportionation rate is likely due to electrostatic repulsion between the highly negatively charged U(V)-triscarbonate complexes (-5 e-).

Model results are shown to be consistent with spectroscopic results on uranyl-organic complexes as well provided the first solvation shell around these complexes is included in the model. In particular, the NMR spectra collected in this study are consistent with an outer-sphere uranyl-catechol complex. The oxidation of catechol by UO22+(aq) was shown to occur through a H-radical mechanism as two phenolic H atoms are transferred in sequence to the axial O atoms of the UO22+. This results in a U(IV)aq and quinone. The intermediate quinone radical species can explain the observation of catechol oxidation and polymerization in the presence U(VI) in aqueous solutions [101].

For uranyl-cell surface complexation, uranyl is predicted to favor binding at phosphonate groups rather than phosphodiester groups. Although the inner-sphere bidentate configuration is predicted to have the lowest Gibbs free energy in these models, the differences between these configurations and outer-sphere associations is relatively small suggesting that a significant portion of the observed complexation could involve outer-sphere binding.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

JDK carried out the quantum mechanical calculations and wrote these portions of the paper. PJ collected the Raman spectra. GPH wrote the Raman methods and results sections. BLP collected NMR spectra and wrote these sections of the paper.