Francisco Colmenero*a,
Jakub Plášilb,
Vicente Timóna and
Jiří Čejkac
aInstituto de Estructura de la Materia (IEM-CSIC), C/ Serrano, 113, 28006 Madrid, Spain. E-mail: francisco.colmenero@iem.cfmac.csic.es
bInstitute of Physics ASCR, v.v.i., Na Slovance 2, 182 21, Praha 8, Czech Republic
cMineralogicko-petrologické oddělení, Národní muzeum, Cirkusová 1740, 193 00 Praha 9, Czech Republic
First published on 27th August 2020
The determination of the full crystal structure of the uranyl sulfate mineral uranopilite, (UO2)6(SO4)O2(OH)6·14H2O, including the positions of the hydrogen atoms within the corresponding unit cell, has not been feasible to date due to the poor quality of its X-ray diffraction pattern. In this paper, the complete crystal structure of uranopilite is established for the first time by means of first principles solid-state calculations based in density functional theory employing a large plane wave basis set and pseudopotential functions. The computed unit-cell parameters and structural data for the non-hydrogen atoms are in excellent agreement with the available experimental data. The computed X-ray diffraction pattern is also in satisfactory agreement with the experimental pattern. The infrared spectrum of uranopilite is collected from a natural crystal specimen originating in Jáchymov (Czech Republic) and computed employing density functional perturbation theory. The theoretical and experimental vibrational spectra are highly consistent. Therefore, a full assignment of the bands in the experimental infrared spectrum is performed using a normal mode analysis of the first principles vibrational results. One overtone and six combination bands are recognized in the infrared spectrum. The elasticity tensor and phonon spectra of uranopilite are computed from the optimized crystal structure and used to analyze its mechanical stability, to obtain a rich set of elastic properties and to derive its fundamental thermodynamic properties as a function of temperature. Uranopilite is shown to have a large mechanical anisotropy and to exhibit the negative Poisson's ratio and negative linear compressibility phenomena. The calculated specific heat and entropy at 298.15 K are 179.6 and 209.0 J K−1 mol−1, respectively. The computed fundamental thermodynamic functions of uranopilite are employed to obtain its thermodynamic functions of formation in terms of the elements and the thermodynamic properties of a set of chemical reactions relating uranopilite with a representative group of secondary phases of spent nuclear fuel. From the reaction thermodynamic data, the relative stability of uranopilite with respect to these secondary phases is evaluated as a function of temperature and under different hydrogen peroxide concentrations. From the results, it follows that uranopilite has a very large thermodynamic stability in the presence of hydrogen peroxide. The high stability of uranopilite under this condition justify its early crystallization in the paragenetic sequence of secondary phases occurring when uranium dioxide is exposed to sulfur-rich solutions.
Uranyl sulfate minerals are invariably found in uranium deposits where uraninite commonly coexists with sulfide minerals. Uranopilite is frequently found on the surface of altered uraninite or in its vicinity.6,8 In these deposits, the oxidative dissolution processes of sulfide minerals lead to the formation of sulfur rich solutions with low pH known as acid-mine drainage.23–26 These acid solutions, generated after mining and ore processing contain high concentrations of metals and toxic elements (such as U, Cd, Cr, Pb and As) and drain from stockpiles of waste rock, mill tailings and open pits at or near the mine site, into surface streams, rivers and lakes leading to the deterioration of the water quality. The precipitation and dissolution of uranyl sulfate minerals are among the most important processes limiting the solubility of uranium in these deposits. The dissolution of these minerals leads to the formation of highly mobile aqueous uranyl sulfate complexes27–40 which may be transported long distances away from the deposit and, therefore, control the migration of uranium at the large scale.41–44 The oxidative dissolution processes of sulfide minerals and the presence of large metal water concentrations in soils with high sulfate contents can also occur in the absence of mining and may have both anthropogenic and non-anthropogenic origin.23,45–50 The presence of sulfur compounds48,49 and biogenic processes50 were fundamental for the genesis of many uranium ore deposits.51 The study of the precipitation/dissolution processes of uranyl sulfate minerals is also highly relevant for the performance assessment of spent nuclear fuel repositories.52–57
In order to make progress in all the fields related to the confinement of uranium containing wastes, a substantial amount of research should be made to perform large scale numerical modeling calculations to identify the main geochemical processes involved and to assess the effectiveness of the barrier systems used to prevent the long term migration of contaminants.58,59 Numerical simulations are the only method that allows to recognize the highly multidisciplinary nature of the problem and takes in account the complex coupled interaction between geology, hydrology, mechanics, chemistry, biogeochemistry and thermodynamics for the precise geometry of the site being studied. Appropriate numerical simulations using highly developed and augmented reactive transport codes are an essential ingredient of the long-term performance assessment of nuclear waste repositories.28,29,60–65 A serious weakness in the current knowledge base is the lack of reliable temperature dependent thermodynamic information for the most important chemical species and chemical reactions involved.31,66–68 The elaboration of a complete chemical thermodynamic database27,28,69,70 is mandatory if we desire to bring the results of such simulations in agreement with the experimental information. The stabilities and the dissolution reaction constants of a given mineral are dependent on the temperature, composition and local conditions (pH and EH), and the knowledge of its Gibbs free energy, enthalpy, and entropy thermodynamic functions of formation, including their variation with temperature, is required for their prediction. The case of the uranyl sulfate minerals is particularly important because, saving the minerals of the zippeite group,71,72 there is not thermodynamic information for the most of these environmentally important minerals even at room temperature. The need of thermodynamic data for uranyl sulfate minerals, due to its enormous relevance under acidic conditions, has been recognized in the past29,73 and has still not been addressed. The availability of accurate thermodynamic data for uranyl sulfate minerals should allow to obtain a deep knowledge about the formation processes of these minerals and of the interaction of uranium and sulfide minerals at several time scales and environmental conditions which is currently poorly understood.74 One of the main purposes of this paper is to provide precise temperature dependent thermodynamic data for one of the most important uranyl sulfate minerals, uranopilite,1–6 and, from them, to investigate its connection with the dissolution/precipitation processes occurring when uranium dioxide is exposed to oxidizing conditions in the presence of SO42− anions. Uranopilite has very high uranium content and is the earliest uranyl sulfate mineral appearing in the paragenetic sequence of uranyl-containing minerals occurring due to the corrosion of uranium ore deposits when the ore is exposed to oxidizing low pH sulfur-rich solutions. It should be noticed that a closely related mineral phase, referred to as metauranopilite, is usually described as its dehydration product. However, this mineral phase, described from a single specimen by Nováček75,76 in 1935, is dubious.8,77 Another, closely related mineral is jáchymovite,78 (UO2)8(SO4)(OH)14·13H2O. Despite of the absence of structural data, it is clear that this mineral differs significantly from uranopilite in the U:S ratio and total water content.
In this work, the infrared spectrum, the elasticity tensor and the phonon spectra of uranopilite were determined from the computed equilibrium crystal structure by using large-scale first-principles calculations (see the Appendix A of the ESI†). The experimental infrared spectrum of uranopilite was collected from a natural mineral specimen originating in Jáchymov (Czech Republic). The collected spectrum was found to be in excellent agreement with the theoretical spectrum and, therefore, all the infrared bands were assigned to specific vibrational atomic motions by performing a normal coordinate analysis of the theoretical vibrational results. The study of the spectral signatures of uranium containing materials is crucial in many research fields ranging from nuclear forensics to nuclear waste storage for their precise identification.10,79,80 Therefore, a detailed spectroscopic characterization is given. From the computed elasticity tensor, the mechanical stability of uranopilite crystal structure was analyzed and a rich set of relevant mechanical properties of this mineral, unknown as far as we know, were determined. The fundamental thermodynamic properties of uranopilite as a function of temperature were derived from the computed phonon spectra. As the mechanical properties, the thermodynamic parameters of uranopilite have not being measured experimentally and are reported here for the first time. The computed fundamental thermodynamic functions were then employed to obtain the thermodynamic properties of formation of uranopilite in terms of the elements and the thermodynamic properties of reaction of a series of reactions involving uranopilite and a subset of the most prominent secondary phases of spent nuclear fuel (SNF). From the temperature dependent thermodynamic reaction data, the relative thermodynamic stability of uranopilite with respect to these secondary phases of SNF were obtained as a function of temperature and under the presence of different concentrations of hydrogen peroxide (H2O2).
Fig. 2 The crystal structure of uranopilite: (A) view of two uranyl sulfate layers and their interlayer space from [100] direction; (B) view of a single uranyl sulfate layer from [001] direction. |
As shown in Fig. 2B, the uranyl sulfate units share the two free vertices of the sulfate tetrahedra to form neutral chains expanding along [100] crystallographic direction. These chains are between three and four uranyl-polyhedra wide and, therefore, are better described as bands. The different bands, although slightly tilted, are placed around a given plane and linked by means of hydrogen bonds to form complex layers (or pseudo-layers) parallel to the (110) plane. The space between the layers is occupied by water molecules. Two uranopilite parallel layers and their interlayer space are shown in Fig. 2A. Fig. 2B displays a single layer from [001], the direction perpendicular to the layers.
The computed unit-cell parameters of uranopilite and the corresponding unit-cell volume and density are reported in Table 1. In this table, the experimental values obtained up to date are also given for comparison. The computed values using the DFT-D2 approach are in very good agreement with the experimental values. The deviation between the calculated and experimental volumes is only about 0.4% with respect to the value reported by Burns2 and 1.6% with respect to the value reported by Plášil et al.6 The inclusion of dispersion corrections improves substantially the computed crystal structure since, as can be observed in the first row of Table 1, if the uncorrected PBE functional is used, the error in the computed volume increases up to 4.8% with respect to the volume given by Burns.2
Parameter | a (Å) | b (Å) | c (Å) | α (deg) | β (deg) | γ (deg) | Vol. (Å3) | ρ (g cm−3) |
---|---|---|---|---|---|---|---|---|
PBE | 9.0171 | 14.3018 | 14.6067 | 93.74 | 99.67 | 101.70 | 1808.69 | 3.861 |
DFT-D2 | 8.9466 | 14.1497 | 14.2549 | 94.57 | 99.16 | 101.59 | 1733.59 | 4.028 |
DFT-D2+U | 8.9467 | 14.1499 | 14.2552 | 94.57 | 99.16 | 101.59 | 1733.67 | 4.028 |
Exp.6 | 8.8556(9) | 13.9819(15) | 14.3070(30) | 96.749(12) | 98.754(12) | 99.726(9) | 1706.91 | 4.091 |
Exp.2 | 8.896(2) | 14.029(3) | 14.339(3) | 96.610(4) | 98.472(4) | 99.802(4) | 1726.03 | 4.046 |
Exp.5 | 8.901(2) | 14.042(3) | 14.521(3) | 97.41(3) | 98.97(3) | 99.69(3) | 1744.58 | 4.003 |
Exp.3 | 8.896(6) | 14.025(9) | 14.299(6) | 96.68(4) | 98.60(6) | 99.92(6) | 1719.01 | 4.062 |
Exp.4 | 8.857(6) | 13.975(8) | 14.335(4) | 96.70(4) | 98.63(4) | 99.56(6) | 1711.55 | 4.080 |
As can also be seen in Table 1 (third row), the Hubbard correction, which was used with a value of the effective Hubbard parameter81 of Ueff = 4.0 eV, does not introduce significant variations of the computed crystal structure. The insignificance of the Hubbard correction may be understood from Fig. S.2 and S.3 of the ESI† in which the computed electron density of states of uranopilite and the frontier band functions (the highest occupied and lowest unoccupied one-electron wavefunctions) are plotted, respectively. The DFT-D2 approach predicts correctly that uranopilite is an insulator with a band gap of 2.27 eV (Fig. S.2†). Furthermore, the highest occupied band function is mainly composed of contributions localized in p orbitals of oxygen atoms and d orbitals of uranium and the lowest unoccupied wavefunction is highly localized in f orbitals of uranium atoms. Because the f orbitals in uranium are empty and separated by almost 2.3 eV from the highest occupied band function, the Hubbard correction, improving the description of the strong correlation between electrons occupying f orbitals, is not needed. As an additional check, the infrared spectrum of uranopilite was determined using the Hubbard correction and the results were essentially the same as those reported in Sec. 2.4 obtained without the Hubbard correction. Due to the small impact of this correction for uranopilite, it was not employed to obtain its mechanical and thermodynamic properties.
There are six symmetrically independent uranium atoms (referred to as U1 to U6) and only one independent sulfur atom (S1) in the structure of uranopilite. The interatomic distances derived from the computed crystal structure are reported in Table S.1 of the ESI†. The experimental values of these distances reported by Burns2 are also given in Table S.1.† Both sets of distances are in excellent agreement. The experimental average U–O distances for the upper and lower apical oxygen atoms in the uranyl pentagonal bipyramids are 1.81(2), 1.82(2), 1.80(2), 1.79(2), 1.76(2) and 1.76(2) Å and the corresponding computed values 1.80, 1.80, 1.80, 1.79, 1.80 and 1.80 Å, respectively. The computed values agree also very well with the mean value of 1.79 Å obtained from many well-refined uranyl mineral crystal structures.82 Likewise, the experimental average U–O distances for the oxygen atoms in the equatorial positions in the uranyl bipyramids are 2.41(2), 2.39(2), 2.41(2), 2.39(2), 2.39(2) and 2.41(2) Å, which can be compared with the corresponding theoretical values: 2.40, 2.39, 2.40, 2.44, 2.41 and 2.39 Å, respectively. Again, these values agree well with the value of 2.37 Å obtained by Burns et al.82 by averaging the corresponding values for a numerous set of well-refined structures involving pentagonal bipyramidal coordination on uranium atom. The theoretical and experimental average S–O distances are the same, 1.49 Å. The computed full crystal structure of uranopilite is included in the ESI† in a file with CIF (Crystallographic Information File) format.
Hydrogen bond | O1⋯O2 | H⋯O2 | A(O1–H⋯O2) |
---|---|---|---|
O1(hydroxyl)–H⋯O2 | |||
OH19–H19···OW33 | 2.668 | 1.663 | 175.20 |
OH20–H20···OW29 | 2.788 | 1.816 | 170.17 |
OH21–H21···OW28 | 2.968 | 1.991 | 175.61 |
OH22–H22···OW32 | 2.812 | 1.822 | 172.91 |
OH23–H23···OW26 | 3.030 | 2.074 | 164.47 |
OH24–H24···OW30 | 2.861 | 1.884 | 171.79 |
O1(water)–H⋯O2 | |||
OW25–H25a⋯OW31 | 2.585 | 1.622 | 158.44 |
OW25–H25b⋯OW34 | 2.664 | 1.681 | 166.97 |
OW26–H26b⋯OW31 | 2.727 | 1.745 | 166.50 |
OW26–H26a⋯OW35 | 2.738 | 1.771 | 162.33 |
OW27–H27b⋯OW25 | 2.686 | 1.694 | 171.36 |
OW27–H27a⋯OW37 | 2.571 | 1.585 | 164.43 |
OW28–H28b⋯O11 | 2.679 | 1.692 | 169.18 |
OW28–H28a⋯O5 | 2.935 | 2.103 | 141.88 |
OW29–H29b⋯OW34 | 3.221 | 2.503 | 130.65 |
OW29–H29a⋯OW38 | 2.593 | 1.587 | 168.84 |
OW30–H30b⋯OW38 | 2.736 | 1.757 | 167.34 |
OW30–H30a⋯OW36 | 2.657 | 1.648 | 176.58 |
OW31–H31b⋯OW32 | 2.703 | 1.713 | 174.45 |
OW31–H31a⋯OH21 | 2.723 | 1.735 | 170.26 |
OW32–H32a⋯OW35 | 2.716 | 1.729 | 168.97 |
OW33–H33a⋯OW34 | 2.727 | 1.740 | 170.06 |
OW34–H34a⋯O8 | 2.977 | 2.109 | 147.16 |
OW34–H34b⋯OH24 | 2.757 | 1.766 | 171.55 |
OW35–H35b⋯OW33 | 2.659 | 1.673 | 167.39 |
OW35–H35a⋯O12 | 2.913 | 2.273 | 122.43 |
OW36–H36a⋯O1 | 2.771 | 1.817 | 161.70 |
OW36–H36b⋯O5 | 2.789 | 1.831 | 164.01 |
OW37–H37b⋯O3 | 2.760 | 1.776 | 174.28 |
OW37–H37a⋯OW36 | 2.730 | 1.764 | 163.17 |
OW38–H38a⋯O10 | 2.721 | 1.793 | 155.98 |
OW38–H38b⋯OW37 | 2.682 | 1.690 | 171.67 |
From the positions of the non-hydrogen atoms, Burns2 suggested a possible hydrogen bond distribution in uranopilite using crystal chemical considerations. The hydrogen bond structure afforded by the first principles calculations is slightly different to the qualitative picture inferred by Burns.2 The differences can be mainly attributed to temperature effects, because the computed structure corresponds to zero temperature and the experimental one to room temperature, and to the qualitative nature of the picture inferred by Burns.2 The determination of the energy-optimized positions of the hydrogen atoms in a given crystal structure is a significant addition because it specifies the hydrogen bond structure in a precise way. Furthermore, from the theoretical point of view, it allows to determine the vibrational spectra and material properties using the first principles methodology.
Fig. 3 X-ray powder diffraction patterns of uranopilite. The upper diffractogram was determined from the computed geometry and the lower one from the experimental one.2 |
Fig. 4 Experimental and theoretical infrared spectra of uranopilite: (A) region 600–1000 cm−1; (B) region 1000–1200 cm−1; (C) region 1450–1750 cm−1; (D) region 2700–3700 cm−1. |
The infrared spectrum of uranopilite samples from several different localities was measured in a previous work by Frost et al.84 The vibrational frequencies obtained in the present work are compared with those reported by Frost et al.84 in Table S.3 of the ESI.† While the obtained frequencies are quite consistent, the results provided in this work are apparently more exhaustive and the rigorous assignment of the infrared bands afforded from the first principles calculation is quite different from the empirical assignment performed by Frost et al.84 In fact, the assignment by Frost et al.84 is very incomplete and the authors called the attention on the great difficulty of assigning many of the infrared bands due to its overlap with other bands. Furthermore, the existence of an overtone and a large number of combination bands (described in detail here) was not noticed by Frost et al.84
The experimental and computed infrared band vibrational frequencies and the corresponding calculated intensities and assignments are given in Table S.4 of the ESI.† Additionally, Fig. S.5 (ESI†) shows images of the vibrational atomic motions in a selection of the infrared active normal modes of uranopilite. The detailed analysis of the infrared spectrum was divided in five main different spectral regions: (i) The region from 2200 to 3750 cm−1 (Fig. 4D, 5C and D); (ii) the region from 1300 to 1750 cm−1 (Fig. 4C and 5B); (iii) The region from 1000 to 1200 cm−1 (Fig. 4B and 5A); (iv) the region from 600 to 1000 cm−1 (Fig. 4A); and (v) the region from 400 to 600 cm−1. The last region, from 400 to 600 cm−1, is analyzed using the theoretical approach only. The study of the bands from this region is needed to identify the origin of some bands from the other spectral regions. The individual infrared bands pertaining to each spectral region are provided in separate sections of Table S.4.† The detailed analysis of the infrared spectrum of uranopilite and its assignment are given in Appendix B of the ESI.†
Fig. 5 Overtone and combination infrared bands in uranopilite: (A) region 1000–1020 cm−1; (B) region 1300–1350 cm−1; (C) region 2200–2500 cm−1; (D) region 3660–3740 cm−1. |
Property | Value | |
---|---|---|
B | Bulk modulus | 21.52 |
G | Shear modulus | 10.94 |
E | Young modulus | 28.06 |
N | Poisson ratio | 0.28 |
D | Ductility index | 1.97 |
H | Hardness index | 0.67 |
AU | Universal anisotropy | 6.65 |
EOS | ||
B | Bulk modulus | 20.74 ± 0.44 |
B′ | Bulk modulus first derivative | 7.30 ± 0.58 |
B′′ | Bulk modulus second derivative | −2.21 ± 0.50 |
(1) |
(2) |
The fitting parameters, B, B′ and B′′ in eqn (1), are the bulk modulus and its first and second derivatives with respect to pressure. The values found are reported in Table 3. The bulk moduli obtained from the EOS and from the computed elastic constants, B = 20.74 ± 0.44 GPa and B = 21.52 ± 1.32 GPa, respectively, are in very good agreement.
The underlying structural mechanism responsible for the NLC phenomenon in uranopilite and the precise values of the compressibility are studied in the next section by means of the analysis of the deformation of uranopilite crystal structure under pressure. It must be noticed that uranopilite shows a perfectly normal mechanical behavior under isotropic pressure. As shown in Fig. S.6 of the ESI,† the volume and unit-cell parameters of uranopilite decrease under isotropic compression up to 5 GPa. The calculated values of the volumetric compressibility, kV = −1/V(∂V/∂P)P, in the pressure range from −0.5 to 4.0 GPa are provided in Table S.6 and plotted in Fig. S.7.† The presence of negative compressibility under anisotropic pressure and its absence under hydrostatic pressure has been also observed recently for the zinc and cadmium oxalates.100
Fig. 7 Unit cell volumes and compressibilities (kV) of uranopilite calculated as a function of the external pressure applied along the minimum compressibility direction. |
In order to understand the structural variations leading to the NLC effect in uranopilite, the variation of the interatomic distances and angles and lattice parameters for the optimized structures corresponding to different pressures applied in the minimum compressibility direction were analyzed. The changes observed in the interatomic distances and angles within the uranium and sulfur coordination polyhedra and the hydrogen bond parameters were quite small. The largest structural variations occurring under the application of pressure were observed in the interlayer distance. The c lattice parameter in uranopilite is exactly equal to two times the interlayer space. As can be seen in Table S.7,† from P = 0.0029 GPa to P = 0.0073 GPa, the unit cell volume decreases mainly due to the decrease of the a and c lattice parameters. However, from P = 0.0073 GPa to P = 0.0131 GPa, while the a lattice parameter decreases, the interlayer space (and the c lattice parameter) increases leading to the increase of the total unit-cell volume. Therefore, the increase of the interlayer space is the main reason for the NLC effect in uranopilite. It is a physical effect related to the increase of the non-bonding distances along the [001] crystallographic direction under pressure.
The structure of uranopilite was also optimized for twelve different pressures directed along the minimum Poisson's ratio direction, ULν-min. The resulting unit cell volumes are reported in Table S.9† and plotted in Fig. S.8 of the ESI.† In this case, the unit cell volume increases from P = 0.014 GPa to P = 0.024 GPa. Hence, negative values of the compressibility are present in this pressure range as can be observed in Fig. S.8.† The calculated values of compressibility from P = 0.003 GPa to P = 0.029 GPa are given in Table S.10.† The minimum value of the compressibility is kV = −34.2 TPa−1 and occurs at P = 0.019 GPa. Therefore, uranopilite displays also significant negative compressibilities for external pressures applied along the direction of minimum Poisson's ratio. The variation of the geometric parameters under pressure was also analyzed and it was found that the increase of the interlayer space was, as in the case of pressures applied along the direction of minimum compressibility, the main responsible of the increase of volume.
T (K) | ΔfH | ΔfG | LogK |
---|---|---|---|
298.15 | −2163.76 | −1905.40 | 333.81 |
300 | −2163.03 | −1903.02 | 331.34 |
350 | −2146.80 | −1842.35 | 274.95 |
400 | −2135.74 | −1787.25 | 233.39 |
450 | −2128.41 | −1736.25 | 201.53 |
500 | −2123.91 | −1688.40 | 176.38 |
550 | −2121.64 | −1643.04 | 156.04 |
600 | −2121.17 | −1599.71 | 139.26 |
650 | −2122.19 | −1558.05 | 125.20 |
700 | −2124.46 | −1517.82 | 113.26 |
750 | −2127.81 | −1478.79 | 102.99 |
800 | −2132.09 | −1440.81 | 94.07 |
UO3(cr) + 1/6SO3(g) + 17/6H2O(l) → 1/6(UO2)6(SO4)O2(OH)6·14H2O(cr) | (A) |
UO3·H2O(cr) + 1/6SO3(g) + 11/6H2O(l) → 1/6(UO2)6(SO4)O2(OH)6·14H2O(cr) | (B) |
UO2CO3(cr) + 1/6SO3(g) + 17/6H2O(l) → 1/6(UO2)6(SO4)O2(OH)6·14H2O(cr) + CO2(g) | (C) |
1/2(UO2)2(SiO4)·2H2O(cr) + 1/6SO3(g) + 11/6H2O(l) → 1/6(UO2)6(SO4)O2(OH)6·14H2O(cr) + 1/2SiO2(cr) | (D) |
Fig. 9 Computed Gibbs free energies of reaction and corresponding reaction constants for the reactions (A) to (F) as a function of temperature. |
1/6(UO2)6(SO4)O2(OH)6·14H2O(cr) + 7/48O2 → 1/8[(UO2)8O2(OH)12]·12H2O(cr) + 1/6SO3(g) + 7/24H2O(l) + 7/24H2O2(l) | (E) |
1/6(UO2)6(SO4)O2(OH)6·14H2O(cr) + 7/6H2O2(l) → (UO2)O2·4H2O(cr) + 1/6SO3(g) + 1/12 O2(g) | (F) |
Under H2O2-free conditions (Fig. 10A), the most stable phase from 300 to 500 K is the uranyl silicate soddyite. At 300 K, uranopilite is the fourth more stable mineral phase after soddyite, rutherfordine and becquerelite. Hence, uranopilite should be replaced by other phases in the absence of H2O2 and in the presence of, for instance, silicate or carbonate ions. However (Fig. 10B), uranopilite is highly stabilized under the presence of H2O and H2O2 and becomes the most stable phase among those studied in this paper. Lastly, as shown in Fig. 10C, under high concentration of H2O2, the uranyl peroxide hydrate phase studtite is the most stable one in line with the results obtained in previous papers16,17,20,21,67,68 and uranopilite becomes the second most stable secondary phase of SNF among those considered in the present paper. However, this is true only for temperatures lower than 148 °C (421 K), and for larger temperatures uranopilite is more stable than studtite. Due to the importance of these results, the precise experimental measurement of the standard state enthalpy of formation of uranopilite in terms of the elements appears to be required to reduce the uncertainty in the estimate used here. Although at temperatures higher than 148 °C, it may be thought that some degree of dehydration of uranopilite may occur, it should be emphasized that the present results refer to the condition of presence of H2O2 and absence of H2O. Under this condition, the most hydrated phases appear to be more stable than the corresponding phases appearing due to its dehydration at H2O2-free conditions. For example, as showed in Fig. 10E, studtite is more stable than metastudtite from 300 to 500 K68 and the same occurs for schoepite and metaschoepite.17 Furthermore, mineral metauranopilite is highly dubious8,77 and no one report on the existence of this phase has been published since 1941.75,76 Theoretical computations aimed to find this phase have also being unsuccessful so far. The thermal decomposition curve of uranopilite, reported in Appendix C of the ESI† does not support the existence of this phase due to the monotonous dehydration thermogravimetric (TG) curve.
The access of the energy-optimized crystal structure facilitated the determination of the elasticity matrix and phonon spectra of uranopilite employing the first principles methodology. From the calculated elasticity matrix, the mechanical stability of the crystal structure of uranopilite was stablished and an abundant set of elastic properties were determined. Uranopilite was found to have a large mechanical anisotropy and to exhibit the negative Poisson's ratio and negative linear compressibility phenomena. The mechanism responsible for the negative linear compressibility effect of uranopilite under the effect of external pressures applied along the direction of minimum compressibility was investigated by analyzing the deformation of the crystal structures under pressure. The anomalous expansion of the interlayer space was found to be the main cause of the increase of the unit cell volume of uranopilite under the effect of compressive pressure. The minimum value of the compressibility was found to be kV = −77.9 TPa−1, for an applied external pressure of 0.011 GPa.
The fundamental thermodynamic functions of uranopilite, unknown so far, were obtained from the computed phonon spectra as a function of the temperature. These thermodynamic functions were used to determine the thermodynamic functions of formation of uranopilite in terms of the elements and the thermodynamic functions of reaction for a set of reactions relating uranopilite and a subgroup of the most prominent secondary phases of spent nuclear fuel. The thermodynamic reaction data was used in order to obtain the relative thermodynamic stability of uranopilite with respect to these phases as a function of temperature and under different concentrations of hydrogen peroxide. The results show that uranopilite is highly stable under the presence of hydrogen peroxide and becomes the most stable phase under the simultaneous presence of water and hydrogen peroxide. Besides, under the presence of hydrogen peroxide and absence of water, uranopilite is the second most stable phase after studtite at low temperatures and become the most stable phase at 148 °C. The high stability of uranopilite under these conditions explains its early crystallization in the paragenetic sequence of secondary phases occurring when natural uraninite is exposed to low-pH sulfur-rich solutions. From the point of view of nuclear waste disposal, sulfur-free conditions should be insured in a repository to avoid the possible formation of uranyl sulfate minerals. With the attenuation of the intensity of the strong radiation fields associated to the SNF over time, the concentration of hydrogen peroxide should decrease significantly and these minerals should dissolve into water forming very mobile uranyl sulfate complexes.
Footnote |
† Electronic supplementary information (ESI) available: (a) Calculated positions of the atoms within the unit cell of uranopilite. (b) Experimental and first principles methods employed in this work; (c) detailed description of the experimental and computed infrared spectra; (d) thermal decomposition curve (Fig. S.1); (e) computed density of states and frontier band functions; (f) interatomic distances; (g) most intense reflections in the X-ray powder pattern; (h) resolution of the infrared bands into single components; (h) depictions of the vibrational motions associated to a subset of the most important infrared active normal modes (i) comparison of the experimental infrared band wavenumbers and those from the literature (j) computed elastic constants; (k) calculated unit-cell volume and parameters as a function of the applied external pressures; (l) compressibility values as a function of the external pressure; (m) calculated fundamental thermodynamic functions; (n) thermodynamic properties of reaction. CCDC 1961689. For ESI and crystallographic data in CIF or other electronic format see DOI: 10.1039/d0ra04596a |
This journal is © The Royal Society of Chemistry 2020 |