Introduction

The share of nuclear energy has been continuously growing in order to match the high demand of energy with low CO2 emission. Consequently, there has been a growing need of exploitable natural uranium resources1. Due to its abundance, uraninite (UO2+x) is the uranium-bearing phase that is mostly used in the extraction of uranium from ores. The majority of processes used to extract uranium from uraninite involves an acid dissolution step coupled with an oxidation of the tetravalent uranium2,3. The economic impact of this step is very important. Its optimization is an absolute prerequisite for the cost-effective exploitation of the uranium deposit. With the exhausting easy exploitable ores, other uranium containing minerals are actually under consideration. Among them, brannerite (UTi2O6) is considered to be the third mineral of tetravalent uranium after uraninite and coffinite (USiO4). Consequently, it could be considered as another exploitable source of uranium4.

Brannerite was found in several places around the world such as Mount Isa (Australia)5,6, Elliot Lake (Ontario, Canada)7,8, Crocker Well (Australia)9, in some different uraniferous deposits in the Witwatersrand area (South Africa)10 and in the exploited uranium mines of Ranger and Olympic Dam (Australia)6. Natural brannerite specimens are found in an amorphous state because of self-radiation (i.e. metamictisation)11,12. The recovery of the crystal network requires thermal annealing at 1000 °C. Brannerite mineral is able to incorporate various elements, among which Pb, Ca, Th, Y, and Ce in the U-site and Si, Al, Fe in the Ti-site. It is worth noting that the presence of lead is mainly radiogenic, consequently to the decay chains of U and Th12,13. The general formula is close to [U, Ca, Y, Ce, La][Ti, Fe]2O611,14,15,16. The crystallized form of brannerite adopts the monoclinic system in the C2/m space group. The structure of brannerite is characterized by 2D infinite layers formed of TiO6 octahedra linked to each other by UO6 octahedra through the c axis.

In the literature, brannerite is reported to be difficult to dissolve. Indeed, the recovery of uranium requires the use of aggressive acidic conditions. The major part of research activities was dedicated to dissolution of natural brannerite, and few were dedicated to synthetic samples11,17,18,19,20,21,22,23,24,25,26,27,28,29. Three parameters have been shown to influence the dissolution rate of natural or synthetic brannerite: temperature, sulfuric acid concentration and redox potential of the leaching solution. From the results reported in the literature, it appeared that the complete recovery of uranium from brannerite could be reached within 1 h in 1 mol/L H2SO4 solution at high temperature (i.e 90 °C) and in the presence of 2.8 g/L (i.e., 0.05 mol/L) of Fe(III)20.

Based on the results reported, the operating temperature plays a major role in the rate of extraction of uranium and titanium. Using a natural sample from Jakutia deposit (Russia), Gogoleva25 determined two different activation energies (EA) depending on the temperature range for the dissolution of brannerite in sulfuric acid solution (0.5 mol/L) containing 0.01 mol/L of Fe(III). Between 15 and 35 °C, the activation energy calculated from the U release in solution was 50.5 kJ/mol, which indicated that the rate-limiting step was a chemical reaction occurring at the solid/solution interface. On the contrary, between 35 and 90 °C, the activation energy decreased down to 30.3 kJ/mol, suggesting the control of the dissolution kinetics by a diffusion step. When studying natural sample from Dieresis deposit (Spain), Gilligan and Nikoloski20 showed that the activation energy of the brannerite dissolution in sulfuric acid solution doped with Fe(III) also depended on the element considered (U or Ti). Indeed, at low temperatures, the activation energy related to the brannerite dissolution and calculated from U release in solution was 36 kJ/mol, which was significantly lower than that obtained for Ti (i.e., 48 kJ/mol). At higher temperatures, both activation energies were similar and lower (23 kJ/mol), which again suggested the existence of diffusion phenomena controlling the dissolution rate. The transition between these two kinetic regimes occurred at lower temperature for the higher sulfuric acid concentrations. Consequently, the dissolution of brannerite at low temperature, in sulfuric acid media and in the presence of iron could be described by the two following reactions:

$$\begin{array}{l}{\mathrm{UTi}}_2{\mathrm{O}}_6\left( {\rm{s}} \right) + 2\,{\rm{FeSO}}_4^ + \left( {{\rm{aq}}} \right)\leftrightarrows 2\,{\rm{TiO}}_2\left( {\rm{s}} \right) + {\rm{UO}}_2\left( {{\rm{SO}}_4} \right)_2^{2 - }\left( {{\rm{aq}}} \right)\\ \qquad\qquad\quad+\, 2\,{\rm{Fe}}^{2 + }\left( {{\rm{aq}}} \right)\end{array}$$
(1)

With EA(U) = 36 kJ/ mol

$${\rm{TiO}}_2\left( {\rm{s}} \right) + 2\,{\rm{H}}^ + \left( {{\rm{aq}}} \right) + {\rm{SO}}_4^{2 - }\left( {{\rm{aq}}} \right)\leftrightarrows {\rm{TiOSO}}_4\left( {{\rm{aq}}} \right) + {\rm{H}}_2{\rm{O}}$$
(2)

With EA(Ti) = 48 kJ/mol

The first reaction corresponds to the oxidative dissolution of brannerite, U(IV) being oxidized by Fe(III) and solubilized as U(VI) species at the solid/solution interface. This reaction leads to the precipitation of TiO2 (anatase) as a secondary phase. The second reaction, which is the dissolution of highly refractory TiO2 layer at the sample surface is controlled by surface reaction involving protons. Consequently, this reaction shows a greater dependence on the concentration of sulfuric acid. At high temperature, the rapid formation of the TiO2 layer through reaction (1) leads to the passivation of the brannerite surface. The diffusion of reactive species through this layer becomes the rate-limiting step.

Besides the effect of temperature, increasing the sulfuric acid concentration has a positive impact on the dissolution rate of uranium, and has even a stronger effect on the dissolution rate of titanium. It is worth noting that Gogoleva25 reported a partial order of U dissolution reaction related to the total sulfuric acid concentration of 0.69 at 70 °C in the presence of Fe(III) ions (with a concentration of 0.01 mol/L). At low sulfuric acid concentrations (i.e., from 0.1 to 0.25 mol/L) and especially at elevated temperature (96 °C), Gilligan et al.20, noticed that the concentration of titanium in the solution decreased markedly after an initial release. Titanium was found to precipitate as anatase (TiO2) through the hydrolysis of TiO+(aq), which is thermodynamically favored at low pH. This observation was consistent with the results obtained earlier by Costine et al.29, who identified a strong correlation between the amount of unleached uranium and the presence of anatase in brannerite dissolution residues.

Costine et al.29 used natural brannerite samples from the Dieresis deposit to investigate the effect of redox potential and the particle size on brannerite dissolution. The effect of redox potential was investigated by varying the Fe(II)/Fe(III) ratio at a fixed total iron concentration. At 25 and 45 °C, they showed a weak dependence of the uranium extraction on the redox potential while at 65 °C, the extent of uranium extraction reached a maximum at Eh = 460 mV vs. Ag/AgCl. Increasing or decreasing the redox potential led to lower uranium extraction rate. An optimized Fe(III) concentration of 0.01 mol/L was also determined by Gogoleva25 in 0.5 mol/L H2SO4 solution at 70 °C. Below this value, the partial order of the overall U dissolution reaction related to the total Fe(III) concentration was found to be 0.5. Increasing the concentration above this value led to the decrease of U dissolution rate.

Considered as an issue for uranium extraction from ores, the chemical durability of brannerite is also viewed as a major advantage for the immobilization of actinides in the field of radioactive waste management. Thus, brannerite was also considered as a potential waste form and entered the composition of several Synroc families30,31,32 or glass-ceramics33. Several studies28,32,34 investigated the leaching of brannerite-bearing waste forms under conditions representative for waste disposal. The release of uranium from synthetic brannerite was modeled as a function of pH and carbonate concentration in oxygenated water by analogy with UO2. The mechanism consisted in two reaction steps: oxidation of U(IV) by dissolved O2 at the surface of brannerite, then subsequent detachment of U(VI) species in the solution. Multiparametric rate laws have been established in order to quantify the promoting effect of surface coordination with protons or carbonate species. The partial order of the overall reaction related to proton activity was found between 0.4 and 0.45 for temperature ranging between 50 and 90 °C. Incongruent dissolution was observed, aqueous Ti concentration being limited by the solubility of Ti oxides or hydroxides over the whole range of conditions studied. TiO2-rich surface alteration products were evidenced at the solid/solution interface. However, under these conditions, the existence of saturation phenomena did not influence strongly the rate of release of uranium.

The majority of the reported studies were dedicated to uranium extraction from natural samples in sulfuric acid conditions and in the presence of Fe(III) used as an efficient oxidant of U(IV). Although these studies were of great importance to optimize U extraction from ores and strongly improved our understanding of the dissolution mechanism, dissolution rate laws were rarely established and when it was done, the kinetic parameters obtained appeared to be deposit-specific. In parallel, the kinetics of dissolution of synthetic brannerite considered as a potential nuclear waste form was studied under much less aggressive conditions. In this field, rate laws were derived that allowed to predict the behavior of brannerite-bearing waste forms over a wide range of repository conditions. In this paper, we decided to combine both approaches and to perform a multiparametric study of brannerite dissolution using pure synthetic sample. We investigated the impact of H2SO4 concentration and temperature with or without Fe(III) added to the solution. The aim was to determine kinetic parameters and rate laws characteristic of the brannerite crystal structure, which is hardly achievable with natural samples. The metamictisation, the presence of associated gangue minerals, of oxidized uranium and of substitutions that modify the structure of brannerite all have potential competing or synergetic impacts on the uranium extraction rate that can be further quantified against the results obtained in this work with pure synthetic brannerite.

Results and discussion

Evolution of the normalized weight loss during brannerite dissolution: general behavior

The shape of the brannerite dissolution curves was similar whatever the experimental conditions. As representative examples, the evolutions of the relative mass loss and of the normalized weight loss during the dissolution of synthetic brannerite in 2 mol/L H2SO4 solution in the presence of Fe(III) at 25 and 80 °C are presented in Fig. 1a, b, respectively. At room temperature (Fig. 1a), the presence of an initial dissolution pulse was clearly evidenced. A similar phenomenon was observed during the dissolution of UO2 in diluted nitric acid solutions35. By analogy, this initial fast release of uranium and titanium was attributed to the presence of a small fraction of more reactive powder. The higher reactivity of this fraction of brannerite material could be caused by a lower grain size, and/or by an oxidized layer. Grinding of the sample favors the surface oxidation of U(IV) materials and sieving generally does not allow to avoid the presence of a dissolution pulse. This initial release always corresponded to a relative weight loss in the range between 0.1 and 1%. As this part of the dissolution curve appeared not reproducible, it was not considered anymore to determine the normalized dissolution rate. Then, the normalized weight loss increased linearly until a relative mass loss of 50% was obtained. During this steady step, the specific surface area of the sample could be considered as constant and the normalized dissolution rate was determined. Afterwards (Fig. 1b), a decrease of the normalized dissolution rate was finally observed especially when using very aggressive conditions. In that case, high dissolution progress was obtained within short time scale and this third step was reached. As Fe(III) ions were introduced in solution in large excess, the decrease of the dissolution rate for Δm/m0 ≥ 50% could not be caused by a significant decrease of Fe(III) concentration. This phenomenon was rather attributed to the significant decrease of the reactive surface area of the sample, which is not considered in the calculation of the normalized weight loss (Eq. 4). As, the powdered sample of brannerite was not sieved, a significant fraction of the sample could have a coarse grain size with a low reactivity25. The decrease of the normalized weight loss could be also attributed to passivation caused by the surface precipitation of refractory TiO2. Nevertheless, the formation of TiO2 as a secondary phase appeared in contradiction with the congruent character of the dissolution. Additionally, it was possible to collect and to analyze by PXRD insoluble residues at the end of some of the leaching experiments. Examples of the collected PXRD patterns are presented in Fig. 2 and are compared with the pattern of the as-synthesized brannerite. The apparition of a very slight diffraction peak around 28° in 2θ corresponding to the potential presence of a very low amount of UO2. Besides, the presence of neoformed phases such as TiO2 was not evidenced. The variation of the intensity between the different patterns were attributed to the presence of preferred orientation of the remaining large brannerite grains.

Fig. 1: Released U and Ti during the dissolution of brannerite in 2 mol/L H2SO4 with 2.8 g/L of Fe(III).
figure 1

Evolution of the relative mass losses Δmi/m0 (%) and of the normalized weight losses NL(i) (g/m2) obtained during the dissolution of synthetic brannerite samples in 2 mol/L H2SO4 solution with 2.8 g/L of Fe(III) at 25 °C (a) and at 80 °C (b).

Fig. 2: Powder X-Ray Diffraction (PXRD) patterns of dissolved brannerite samples.
figure 2

PXRD patterns corresponding to the brannerite samples obtained after leaching experiments labeled 2, 5, and 11 (see Table 3).

The normalized dissolution rates and the congruence ratios determined in the second step (i.e., during the linear evolution of NL(i)) are gathered in the Table 1 for all the conditions investigated. From the results reported in this table, it was found that addition of 2.8 g/L of Fe(III) led to the strong increase of the dissolution rate of brannerite for all the experimental conditions examined. However, the influences of temperature and acid concentration were different in the presence (or not) of Fe(III).

Table 1 Normalized dissolution rates RL(U) and RL(Ti) obtained during the dissolution of synthetic brannerite samples in various conditions (with and without Fe(III)). Associated congruence ratios, R = RL (U)/RL(Ti).

In order to evidence the impact of temperature on the dissolution rate in both cases, two sets of dissolution experiments were conducted with various acid concentrations at 25, 40, 60, and 80 °C. The evolution of the normalized weight losses NL(U) and NL(Ti) obtained in 2 mol/L H2SO4 solution without and with Fe(III) are plotted in Figs. 3 and 4, respectively. The results showed that the normalized dissolution rate was strongly dependent on temperature either in the presence or absence of Fe(III). In addition, for the four temperatures studied, the normalized weight losses NL(U) and NL(Ti) were nearly the same, showing that the dissolution was congruent under the experimental conditions investigated.

Fig. 3: Dissolution of brannerite in 2 mol/L H2SO4.
figure 3

Evolution of the normalized weight losses NL(i) (g/m2) (a) and of the relative mass losses Δmi/m0 (%) (b) obtained during the dissolution of synthetic brannerite sample in 2 mol/L H2SO4 solution for various temperatures.

Fig. 4: Dissolution of brannerite in 2 mol/L H2SO4 with 2.8 g/L of Fe(III).
figure 4

Evolution of the normalized weight losses NL(i) (g/m2) (a) and of the relative mass losses Δm/m0 (%) (b) obtained during the dissolution of synthetic brannerite samples in 2 mol/L H2SO4 solution containing 2.8 g/L of Fe(III) for various temperatures.

The evolution of the normalized weight loss of brannerite calculated from the releases of uranium and titanium obtained at room temperature in various sulfuric acid concentrations and without Fe(III) is presented in Fig. 5. Under these conditions, the impact of an increase of the sulfuric acid concentration was found to be low. Indeed, when the acid concentration increased from 0.1 to 2 mol/L, the normalized dissolution rate RL(U) increased only by a factor of 1.6, i.e. from (4.72 ± 0.05) × 10−6 g/(m2 h) to (7.57 ± 0.07) × 10−6 g/(m2 h) (Table 1). Moreover, the examination of the RL(U) values reported in Table 1 showed that the impact of the acid concentration was enhanced with temperature. Increasing the sulfuric acid concentration from 0.1 to 2 mol/L led to the increase of the RL(U) values by a factor of 6.2 at 80 °C. The evolution of the normalized weight loss of brannerite calculated from the uranium and titanium releases obtained at room temperature in various sulfuric acid concentrations and in presence of Fe(III) is shown in Fig. 6. The impact of the sulfuric acid concentration was stronger when introducing Fe(III) in the solution. For instance, at 25 °C, the normalized dissolution rate RL(U) increased by a factor of 9.4, from (3.3 ± 0.1) × 10−5 g/(m2 h) to (3.09 ± 0.06) × 10−4 g/(m2 h), when the acid concentration raised from 0.1 to 2 mol/L (Table 1). Once again, the impact of the acid concentration was even more pronounced at high temperature since an increase of RL(U) by a factor of 41 was observed at 80 °C when the acid concentration increased from 0.1 to 2 mol/L (i.e., from (2.9 ± 0.1) × 10−3 to (1.20 ± 0.03) × 10−1g/(m2 h)). These results indicated that, either in absence or in the presence of Fe(III) in sulfuric acid solutions, the impact of the temperature and of the acidity were not independent. It is also noteworthy that in both cases, the dissolution remained congruent even for low H2SO4 concentrations.

Fig. 5: Dissolution of brannerite in H2SO4 solutions at 25 °C.
figure 5

Evolution of the normalized weight losses NL(i) (g/m2) (a) and of the relative mass losses Δmi/m0 (%) (b) obtained during the dissolution of synthetic brannerite samples in various H2SO4 solutions at 25 °C.

Fig. 6: Dissolution of brannerite in H2SO4 solutions with 2.8 g/L of Fe(III) at 25 °C.
figure 6

Evolution of the normalized weight losses NL(i) (g/m2) (a) and of the relative mass losses Δmi/m0 (%) (b) obtained during the dissolution of synthetic brannerite samples in various H2SO4 solutions containing 2.8 g/L of Fe(III) at 25 °C.

Establishment of multiparametric rate laws

The variations of RL(U) versus T and log(H3O+) are plotted in Fig. 7 when adding or not Fe(III) in the dissolution media. This figure highlights the strong impact of the addition of Fe(III) to the sulfuric acid solution on the kinetics of dissolution in all the experimental domain. The two sets of data were used independently to determine the three parameters of the multiparametric expression (Eq. 8). The fitted activation energy, EA, the normalized dissolution rate constant, k, and the partial order of the reaction relative to the proton activity, n, are summarized in Table 2. In order to show the agreement between experiments and modeling, the data simulated on the whole experimental domain using Eq. (8) and the fitted parameters are plotted as colored maps in Fig. 7. Additionally, Supplementary Fig. S1 of the supporting data shows the comparison of simulated and experimental data.

Fig. 7: Experimental and simulated dissolution rates of brannerite.
figure 7

Simulated maps obtained by fitting the multiparametric expression of the normalized dissolution rate of synthetic brannerite samples based on the data obtained without (red dots) or with (blue dots) Fe(III) in sulfuric acid solutions.

Table 2 Values of log k, EA and n fitted from the RL(U) values obtained in the presence or in the absence of Fe(III) with various acid concentrations and temperatures.
Fig. 8: PXRD and SEM images of as-prepared brannerite sample.
figure 8

Observed, calculated, and the difference patterns after the Rietveld refinement of the crystal structure of brannerite (left panel), SEM image of the brannerite grains used for the dissolution experiments (right panel).

From the analysis of the fitted parameters (Table 2), it appears that the introduction of Fe(III) in the system strongly increased the dissolution rate constant, k, by 5 orders of magnitude. It is now acknowledged that in the presence of oxidative species in solution regarding to U(IV), the dissolution of U(IV)-bearing phases must be considered as a corrosion reaction, in which the oxidant is consumed to convert the insoluble U(IV) species to the much more soluble U(VI) species (as UO22+ (aq)) in the solution. The rate of this reaction should depend on redox conditions, as already demonstrated for uraninite36,37,38,39,40,41,42.

The thermodynamic driving force for such a corrosion process is the difference existing between the redox potential of the solution and the potential associated to the U(IV)/U(VI) couple at the brannerite/solution interface. Increasing the redox potential of the dissolution medium leads to an increase of this driving force. For this reason, oxidants like Fe(III) are added to the leaching solution to improve the kinetics of dissolution of refractory U(IV) phases compared to solutions in equilibrium with atmosphere. However, Costine et al.29 showed a complex dependence of the dissolution rate of brannerite with redox potential of the solution.

In the absence of Fe(III) and operating in less aggressive conditions, Thomas and Zhang34 reported that uranium release was dependent on the activity of dissolved dioxygen. By analogy with UO2, they assumed that the brannerite dissolution mechanism consisted in two successive reactions. The first was the oxidation of U(IV) to U(VI) by O2 (aq) at the solid/solution interface. The second reaction involved the adsorption of protons and/or surface complexation of uranium (VI) sites at the surface, followed by the detachment of uranium (VI) based species. Under the conditions investigated by Thomas and Zhang34, the first reaction was the rate-limiting step, especially when carbonate ions were added to the solution. Carbonates are considered as ligands, which promote the detachment of U(VI) species from the brannerite surface. The comparison of complexation constants of UO22+(aq) with CO32−(aq) or SO42−(aq) gives insights for the magnitude of the ligand-promoted detachment phenomenon. Guillaumont et al.43 reported that log K°(298 K) of the reaction UO22+(aq) + 2 L2−(aq) \(\leftrightarrows\) UO2(L)22−(aq) equals 16.61 ± 0.09 and 3.02 ± 0.38 for L2− = CO32−(aq) and SO42−(aq), respectively. Assuming that the tendency to form more soluble surface complexes increases with the tendency to form aqueous complexes, sulfate ions are less efficient than carbonates in increasing the detachment step. Similar analysis was made by Gilligan et al.24 to explain the lower efficiency of the ferric chloride/hydrochloric acid system than the acidic ferric sulfate media to extract uranium from brannerite. It could be concluded that the oxidation of U(IV) by O2(aq) at the surface sites was the rate-limiting step in our study in the absence of Fe(III). Adding Fe(III) increased the redox potential of the leaching solution and might increase the rate of the first oxidation step.

The magnitude of the EA is also an important indicator of the nature of the rate-limiting reaction44. The obtained values of EA were 71 ± 4 and 91 ± 6 kJ/mol in the absence and in the presence of Fe(III), respectively for the whole temperature range investigated (i.e. 25–80 °C). Such high values were not compatible with diffusion-controlled mechanism, which are usually relatively independent of temperature (with EA < 20–25 kJ/mol). Such EA values were also significantly higher than those reported in the literature9,20,25 for natural brannerite samples in the presence of Fe(III), but in agreement with the values obtained by Zhang et al.32 for synthetic brannerite in the absence of Fe(III) (i.e., 70 kJ/mol). EA value obtained in this work without Fe(III) added to the dissolution media was found to be in the upper range, but still compatible with surface reactions controlling phenomena at the solid/liquid interface (typically between 40 and 80 kJ/mol). This strong discrepancy could be related to the metamictization of brannerite natural samples, which destroys the crystal structure. The complete loss of the crystalline structure could influence the kinetics of brannerite dissolution, or change the nature of the rate-limiting step. Due to radiation damage, U and Ti could be located in independent entities at the atomic scale. Thus, only U-rich entities might be reactive, their specific dissolution leaving behind a surface layer of Ti-rich entities acting as a diffusion barrier for the species in solution responsible for U(IV) oxidation.

The partial order of reaction related to the protons activity obtained in this work (n = 0.42 ± 0.09) in the absence of Fe(III) was in agreement with the value reported by Thomas and Zhang34 when working on synthetic brannerite in less aggressive conditions (i.e. between 0.4 and 0.45 in a temperature range of 50–90 °C). Gogoleva25 also reported a partial order of U dissolution reaction related to the total sulfuric acid concentration of 0.69 at 70 °C in ferric media (Fe(III) = 0.01 mol/L), whereas Gilligan et al24. reported 0.43 and 0.57 for U and Ti, respectively in 0.1 to 2 mol/L H2SO4 solutions at 52 °C. These values were not directly comparable to the value obtained in this work (n = 0.84 ± 0.08). Unlike the formers, the latter was valid in the whole temperature range investigated and was related to the (H3O+) activity in the leaching solution (not to the acid concentration). This n value indicates that the adsorption of protons occurs onto the surface of the brannerite, then protonated U and/or Ti surface species are detached from the brannerite lattice.

As conclusion, multiparametric dissolution experiments were conducted on pure synthetic brannerite samples. Normalized dissolution rates were evaluated under various temperatures (from 25 to 80 °C) and sulfuric acid concentrations (from 0.1 to 2 mol/L), either in the absence or in the presence of 2.8 g/L of Fe(III) in the dissolution media. Multiparametric dissolution rate laws were established with and without Fe(III) and the modeled parameters were found to be valid in the whole temperature and acidity ranges. As expected, the introduction of Fe(III) in the solution strongly increased the normalized dissolution rate of the sample and modified the activation energy and the partial order of the overall dissolution reaction related to protons activity. Contrary to what has been reported in the literature for natural brannerite samples, the dissolution was found to be congruent in the whole experimental domain examined. Both U and Ti releases in solution were controlled by chemical reactions involving the oxidant (either O2(aq) or Fe(III)) and protons at the solid/solution interface. The role of sulfate ions was not explicitly studied, but they are expected to kinetically favor the step of detachment of U(VI) species from the surface of the solid. The two-steps mechanism described from experiments performed with natural samples, which involved the dissolution in acidic media of a neoformed TiO2 layer at the surface of brannerite was not valid for synthetic samples. This strong discrepancy could be related to the metamictization of brannerite natural samples, which destroy the crystal structure. Due to radiation damage, U and Ti could be located in independent entities at the atomic scale. Thus, only U-rich entities might be reactive, their specific dissolution leaving behind a surface layer of Ti-rich entities acting as a diffusion barrier for the species in solution responsible for U(IV) oxidation. The analysis of the alteration layer of natural brannerite grains at the atomic scale could bring key information to support this assumption.

Methods

Synthesis

NH4OH and titanium alkoxyde Ti[OCH(CH3)2]4 were obtained from Sigma Aldrich and were of analytical grade. The uranium metal was obtained from CETAMA, France. However, it was used and the form uranium (IV) tetrachloride solution obtained by dissolving the metal chips in 6 M of HCl solution. The final uranium concentration is about 0.5 mol/L based on ICP-AES analyses.

Pure brannerite (UTi2O6) was obtained in the application of the method developed by Mesbah et al.45. It was inspired from the synthesis protocol proposed by Martinez et al.46 for the synthesis of uranium dioxide. The synthesis consisted in mixing U stock solution (nU = 2.33 mmol) and Ti-alkoxyde (nTi = 4.66 mmol + 3% molar excess) in 17 mL of deionized water. Then, an excess of NH4OH (46.6 mmol) leading to the formation of a brownish gel. The formed precursor was stirred for 30 min then washed twice in water and once in ethanol. The final recovered powder poured in 100 ml of ethanol and dried under vacuum at 40 °C leading to a final precursor with a specific surface area higher than 200 m2/g. The obtained sample was compacted into a green pellet, and placed in a furnace for 72 h at 1300 °C under Argon continuous stream.

Solid characterization

The obtained pellets, were mechanically milled using a ball jar at 30 Hz during 30 min. The specific surface area (SBET) of the ground brannerite pellets was determined using a MJ Micromeritics ASAP 2020 apparatus (BET method) under krypton at 77 K. The sample was previously heated at 333 K for 24 h. The measurement was performed three times and the obtained values ranged from 0.1 to 0.3 m2/g.

The synthesized powder and insoluble residues of brannerite collected at the end of some of the leaching experiments were analyzed by Powder X-Ray Diffraction (PXRD) over the use of the Bruker D8 advance diffractometer (Cu Kα1,2, λ = 1.54184 Å). Special sample holders (dome) were used in order to avoid any radioactive contamination. Data were collected in an angular range comprised between 5 and 100° (2θ) with steps of 0.019 for a counting time of 3 h per pattern. The crystal structure of the synthetic brannerite was refined by applying the Rietveld method using the Fullprof_package47. The obtained data confirmed the preparation of pure brannerite sample with a = 9.8126(2) Å, b = 3.7688(1) Å, c = 6.9243(1) Å, β = 118.94(1)°, V = 224.10(1) Å3. Observed, calculated and the difference patterns after the crystal structure refinement is viewed in Fig. 8 (left panel). Moreover, the morphology of the brannerite sample used for the dissolution experiments was observed by scanning electron microscopy (SEM, Quanta 200, ESEM-FEG, FEI Company). The collected image presented in Fig. 8 (right panel), showed the presence of relative large grains of about 5–10 µm and of smaller debris of about 1 µm. Despite the presence of these smaller particles, the specific surface area of the sample measured by krypton adsorption and the BET method was low (i.e. between 0.3 and 0.1 m2/g). These values of specific surface area corresponded to the geometric surface area of brannerite spheres of 3–9 µm in diameter, respectively, which indicated that the great majority of the brannerite grains were in this range of size.

Dissolution experiments

Samples of synthetic brannerite were dissolved to conduct a multiparametric study. 25 or 50 mg of powdered brannerite were put in contact with 25 or 50 mL H2SO4 solution, with concentrations ranging from 0.1 to 2 mol/L in polytetrafluoroethylene container. Fe2(SO4)3 was heated overnight in order to evacuate water molecules, then was added to H2SO4 solutions to obtain an Fe(III) concentration of 2.8 g/L (i.e., 0.05 mol/L). The experimental conditions considered for all the dissolution tests are summarized in Table 3. All the containers were hermetically closed and placed in different ovens at various temperatures (25, 40, 60, and 80 °C) under agitation for the duration of the experiment.

Table 3 Experimental conditions considered for the dissolution tests performed on synthetic brannerite.

During the dissolution experiments, solution aliquots were collected with a decreasing frequency. 800 μL of solution were sampled then centrifuged at 14,000 rpm for 2 min in order to avoid the presence of colloids in the solution. Then, 500 μL of the supernatant were collected and diluted with 4.5 mL of 1 wt. % HNO3 solution. 300 µL of the residual solution were put into the dissolution reactor along with 500 μL of fresh dissolution medium. By this way, it can be noticed that only 3% of the solution was replaced at each uptake. These conditions avoided possible strong disturbance of the surface over volume ratio during the experiment.

The elemental concentrations of uranium and titanium in solution were measured using inductively coupled plasma-atomic emission spectroscopy (ICP-AES). ICP-AES was first calibrated with U and Ti standard solutions prepared by dilution of a certified standard solution (CU = CTi = 1000 ppm) with 0.2 mol/L HNO3 solution. Four and five emission rays of U were selected (279.394, 367.007, 385.958, and 409.014 nm) for U and (307.864, 323.452, 334.0187, 334.941, and 336.121 nm) for Ti, respectively. The elemental concentrations and associated errors were estimated as the average value and twice the standard deviation value of the 3 replicates of all emission lines, respectively.

Once the concentrations were measured, the relative mass loss, \(\frac{{{\Delta}m}}{{m_0}}\left( t \right)\) (expressed in wt. %) and the normalized weight loss, NL(i,t) (expressed in g/m²) were calculated using the initial surface area of the brannerite sample in contact with the solution and the theoretical mass fraction of the elements using the following equations:

$$\frac{{{\Delta}m}}{{m_0}}\left( t \right) = \frac{{m_i(t)}}{{f_i \times m_0}}$$
(3)
$$N_{\rm{L}}\left( {i,t} \right) = \frac{{m_i(t)}}{{f_i \times S\left( t \right)}}$$
(4)

mi(t) (g) corresponds to the mass of element i collected in the solution at time t; m0 (g) is the initial mass of sample introduced in the system (Table 3). S (m²) is the surface area of brannerite in contact with the solution and fi (g/g) represents the mass ratio of element i in the studied solid.

The surface area of the sample in contact with solution is defined as:

$$S\left( t \right) = S_{{\rm{SA}}}\, \times \left( {m_0 - \frac{{m_i\left( {t - 1} \right)}}{{f_i}}} \right)$$
(5)

where SSA (m²/g) corresponds to the initial surface area mentioned above. The mass of element i dissolved during the experiment was calculated using the concentration of each element, Ci (g/L) according to:

$$m_i = C_i\, \times V$$
(6)

where V (L) is the volume of the acidic solution introduced in the reactor as mentioned in Table 3. The normalized dissolution rate, RL(i) is obtained as the time-derivative of the normalized mass loss:

$$R_{\rm{L}}\left( i \right) = \frac{{{\mathrm{d}}\,N_{\rm{L}}\left( i \right)}}{{{\mathrm{d}}t}}$$
(7)

For all the dissolution experiments, the normalized dissolution rate are determined for 0.1–1% ≤ \(\frac{{{\Delta}m}}{{m_0}}\) ≤ 50%. This domain corresponds to a linear variation of the normalized weight losses unlike the conditions used. Finally, the congruence ratio, R was evaluated as the ratio RL(U)/RL(Ti).

Multiparametric expression of the normalized dissolution rate

The multiparametric description of the dissolution reaction rate was proposed by Lasaga48 based on experimental results. In this study, only the influence of the activity of the proton and of the temperature on the dissolution rate of brannerite was considered. Thus, the general form of the rate law was simplified as follows:

$$R_{\rm{L}}\left( i \right) = k \times {\mathrm{exp}}\left( { - \frac{{E_{\rm{A}}}}{{RT}}} \right) \times \left( {{\rm{H}}_3{\rm{O}}^ + } \right)^n$$
(8)

where k (g/(m2 d)) is the rate constant, EA (kJ/mol) is the apparent activation energy of the overall reaction, R is the gas constant, and T is the absolute temperature. (H3O+) is the activity of proton, and n is the partial order of the reaction related to the proton.

The rate law parameters k, EA and n were fitted using the entire set of RL(U) data obtained either in the presence or in absence of Fe(III) in solution using the JMP V14 statistical software. The nonlinear regression was done using the analytic Gauss-Newton algorithm.

The pH values were calculated from the molalities of the prepared H2SO4 solutions with the geochemical speciation software PHREEQC49 associated with the Thermochimie SIT V9 thermodynamic database50,51. In this database, the log (298 K) associated to the protonation equilibrium H+ + SO42− \(\leftrightarrows\) HSO4 and the enthalpy of reaction reach 1.98 and 22.44 kJ/mol, respectively. The molar to molal conversion was performed using the densities of 0.1 mol/L, 1 and 2 mol/L H2SO4 solutions at the corresponding temperature. The density of acid solutions were calculated using the expression developed by Novotny and Sohnel52. Then, the activity corrections were performed using the Specific Interaction Theory (SIT) as implemented in PHREEQC. The proton activities of the various solutions calculated using PHREEQC software are summarized in Supplementary Table S1.