# Alpha dose rate and alpha decay dose impacts on the residual alteration rate regime of HLW nuclear glasses

Magaly Tribet, Caroline Marques, Sarah Mougnaud, Veronique Broudic, Christophe Jegou, Sylvain Peuget

► To cite this version:

Magaly Tribet, Caroline Marques, Sarah Mougnaud, Veronique Broudic, Christophe Jegou, et al. Alpha dose rate and alpha decay dose impacts on the residual alteration rate regime of HLW nuclear glasses. npj Materials Degradation, 2021, 5, https://doi.org/10.1038/s41529-021-00183-4. 10.1038/s41529-021-00183-4. cea-03298512

HAL Id: cea-03298512

https://cea.hal.science/cea-03298512

Submitted on 23 Jul 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

# Alpha dose rate and decay dose impacts on the long-term alteration of HLW nuclear glasses

Magaly Tribet<fnref n="1" />,<fnref n="1" />✉, Caroline Marques<fnref n="1" />, Sarah Mougnaud<fnref n="1" />, Véronique Broudic<fnref n="1" />, Christophe Jegou<fnref n="1" /> and Sylvain Peuge<fnref n="1" />

In the prospect of deep geological disposal, the long-term behavior of high-level nuclear glasses has to be investigated regarding alpha radiation induced by long-life minor actinides. The present study focuses on the effects of alpha radiation on the long-term chemical reactivity of R7T7-type glasses, by separately considering the alpha dose rate and the alpha decay dose. Old SON68 glasses doped with ²³⁸/²³⁹PuO₂ or ²⁴⁴CmO₂ were studied to simulate high alpha dose rates corresponding to an early water ingress and a high level of alpha decay doses corresponding to long-term disposal conditions. A part of the ²³⁸/²³⁹Pu-doped glass block was annealed to fully recover the irradiation-induced damage accumulated since the glass was fabricated and to dissociate the effect of the alpha dose rate from that of the alpha decay dose. The glasses were then leached under static conditions at 90 °C for several years. The results showed that the residual alteration rate is not affected by the alpha dose rate over a wide range of dose rate values expected under disposal conditions: this glass remained relatively insensitive to the alpha radiolysis phenomena at the glass-water interface. However, over the duration of the experiments, the residual alteration rate of the damaged ²³⁸/²³⁹Pu-doped glass was enhanced compared to that of the annealed glass. This result is in agreement with those obtained on the ²⁴⁴Cm-doped glass and with reported values in the literature on simplified externally irradiated glasses, indicating that the ballistic effects of the recoil nuclei are responsible for this increase in the residual alteration rate.

npj Materials Degradation (2021)5:36; https://doi.org/10.1038/s41529-021-00183-4

## INTRODUCTION

Several countries have chosen to recycle U and Pu in UOX and MOX fuels, and to vitrify the ultimate high-level radioactive waste (HLW) composed of fission products and minor actinides in silicate glasses. In France, the borosilicate glass called R7T7 is currently used to confine the HLW remaining after spent nuclear fuel reprocessing and destined for geological disposal. The safety of such repositories relies on a multi-barrier approach, based on the properties of the glass, metallic containers, and host rock. The physical and chemical stabilities of the glass over a geological timescale must thus be verified by studying the glass response to radiation emitted by the nuclear waste and its behavior when submitted to water dissolution<fnref n="1" />.

The alteration of SON68-type glass (which is the non-radioactive surrogate of R7T7 glass) has been extensively studied since the 1980s. The fundamental leaching mechanisms of this type of glass, such as water penetration and interdiffusion²⁻⁵, hydrolysis of the bonds of the glass network⁶,⁷, condensation of the hydrolyzed species⁷⁻⁹, and precipitation of secondary phases in certain cases⁹,¹⁰ are relatively well described. In the long term under disposal conditions, the alteration solution is expected to be poorly renewed, which will lead to an increase of the poorly soluble glass former elements (Si, Al, Zr...) in solution and the formation of a protective alteration layer at the glass surface. The glass will thus dissolve very slowly, by what is called the residual alteration rate regime¹¹. However, the fundamental mechanisms by which the alteration layer forms and the rate-limiting reactions governing this residual alteration rate regime are still debated and two main models are described in the literature⁷,¹²⁻¹⁸. First, simultaneous hydrolysis and precipitation mechanisms at the glass-water interface are proposed in the interface-coupled dissolution-precipitation (ICDP) model⁷,¹³⁻¹⁵. Such mechanisms lead to inward progress of the leaching front and formation of the

so-called surface alteration layer (SAL). In this model, the residual glass alteration is controlled by this SAL. Second, interdiffusion and ion exchange with water (H⁺/H₃O⁺), followed by progressive in situ hydrolysis and subsequent recondensation⁸,¹²,¹⁶⁻¹⁸ are the proposed mechanisms to explain the formation of the alteration layer. In this last case, the alteration layer region, also called gel, is constantly reorganized over time⁸,¹⁹ and acts as a diffusive barrier. Actually, this is the transport-limiting factor of this barrier which explains the residual alteration rate regime.

Furthermore, the specificity of HLW nuclear glasses is their radioactivity lasting over time and which has to be taken into account to accurately describe the glass aging under disposal conditions. Most studies on the radiation impact on long-term leaching behavior focused on the impact of beta and gamma radiation<fnref n="1" />,²⁰⁻²². However, some recent results evidenced the predominant role of alpha decay of minor actinides on the long-term radiation component, and assessed its impact on the glass structure and properties<fnref n="1" />,²³⁻²⁵. Such results thus highlighted the need for improved evaluation of the impact of this kind of radiation on glass leaching.

Due to the various complexities associated with the radioactive specimen handling and characterization, external ion irradiation is commonly used to simulate the effects of self-irradiation damage²⁵,²⁶. However, in such experiments, ion irradiation and leaching cannot occur simultaneously: the glass is first irradiated and then leached, whereby it is observed that the alteration behavior of borosilicate glasses can be modified by prior external irradiation by heavy ions²⁴,²⁶⁻²⁸ specially during the residual alteration rate regime²⁴,²⁶,²⁷. Radiation damage induced by ion irradiation was found to induce the formation of an alteration layer approximately four times thicker than for non-irradiated glass leached under the same conditions²⁷. The increase in the alteration layer thickness with irradiation dose was attributed to

<footnote n="1">¹CEA, DES, ISEC, DE2D, University of Montpellier, Marcoule, Montpellier, France. ✉email: magaly.tribet@cea.fr</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Si releases and pH evolution versus time for doped glasses leaching tests. Total silicon concentration versus time for annealed and damaged ²³⁸/²³⁹Pu-doped glass leaching tests (pink triangles and green diamonds, respectively) and ²⁴⁴Cm-doped glass leaching test (blue circles). pH ranges at 90 °C, calculated from measured pH at 25 °C, are detailed for each test in order of magnitude. The gray and hatched areas correspond to saturation values of SON68 glass with respect to orthosilicic acid (H₄SiO₄), based, respectively, on Gin et al.⁴¹ or on Jegou et al.⁵⁰ for a pH₉₀°C ≤ 8.97. Red arrows indicate sampling with filtration and ultrafiltration to check for colloids. Errors bars are of 3%.</figcaption></figure>

the effects of nuclear collisions that change the structural and
mechanical properties of the material, enhancing the leach
rate²⁴,²⁷. The rate of altered layer growth diminished with time
for both irradiated and non-irradiated glasses, with square root of
time dependency; this indicates that a diffusive mechanism occurs
irrespective of the radiation condition²⁷.

Recent solid-state characterization of irradiated and then
altered samples²⁶ shows that radiation damage can accelerate
the maturation of the alteration layer, leading to a porous gel-like
microstructure. This observation and the kinetics enhancement
previously mentioned²⁷ are concomitant, but no link between
these two statements is established. Moreover, a higher reactivity
in ¹⁸O (from enriched water) and ¹⁶O (from OH or water contained
in the alteration layer) exchanges at 25 °C is also observed in pre-
irradiated glasses within the reactive interface, compared to the
boron and alkali-depleted silica-rich gel layer (where no higher
reactivity is reported). Furthermore, experiments performed on
heavily altered glasses²⁹ have shown that ballistic collisions from
ion implantation could heal the porous gel part of the altered
glasses to create a pore-free homogeneous microstructure. These
results demonstrate that the nuclear decay dose affects the
alteration layer, and hence, the residual alteration rate, but the
precise origin of these effects has not yet been clearly established.

For radioactive specimens, very few studies that utilize a
simplified medium have been reported in the literature²⁰,³⁰⁻³⁵,
and they do not all focus on the residual alteration rate regime nor
lead to the same conclusions. In diluted and renewed solutions,
neither the alpha dose rate nor alpha decay have been seen to
impact actinide-doped SON68 glasses³⁰,³²⁻³⁴, while ‘dramatic’
differences have been recently reported on a simplified Pu-doped
glass, subject to a high alpha dose rate (~8700 Gy·h⁻¹) compared
to the reference experiment on the non-radioactive sample³¹.
Only one study on the residual alteration rate regime, which
focused on alpha dose rate (150 Gy·h⁻¹), has been reported²⁰, and
this study concluded that no radiation effect was present.

The current study aimed to assess the impact of alpha radiation
on the long-term aqueous behavior of alpha-doped SON68
radioactive glasses by considering the effects of alpha dose rate

and alpha decay dose separately. An old, damaged ²³⁸/²³⁹Pu-
doped SON68 glass was used because of its significant alpha
decay dose considering the structural modifications (14 MGy<sub>nucl</sub>)²³
and its alpha dose rate considering its doping in ²³⁸/²³⁹Pu (a value
of 3500 Gy·h⁻¹ is expected in R7T7-type glass few decades after
elaboration). A leaching experiment was also performed under the
same conditions on a ²⁴⁴Cm-doped SON68 glass (alpha decay
dose of 35 MGy<sub>nucl</sub> and alpha dose rate of 23,500 Gy·h⁻¹) to
complete the alpha radiation description. The dissociation of the
radiation parameters (dose rate versus decay dose) on the
leaching behavior was achieved by annealing a part of the ²³⁸/
²³⁹Pu-doped glass to fully recover the irradiation-induced damage;
the annealed glass had the same alpha dose rate (3500 Gy·h⁻¹) as
the non-annealed one, but a very low alpha decay dose
(<0.1 MGy<sub>nucl</sub>). Actually, the damage on this annealed glass
remained low enough for years, which allowed a targeted study
of pure radiolytic effects on the long-term glass leaching behavior.

RESULTS AND DISCUSSION

Details of leaching experiments, including reference experiments,
are given in the Leaching device and conditions section.

pH and total silicon concentration stabilization in the
leachates

The pH ranges measured in the doped glass leaching tests and
then calculated at 90 °C are presented in Fig. 1. They are compared
to the pH₉₀°C values measured in the reference experiments,
which were in the range of 9.25 ± 0.30 ('ref. 1' experiment) and
8.95 ± 0.30 ('ref. 2' experiment). For these systems under irradia-
tion, the measured pH is around 0.5–1 pH unit lower than that
those measured in the reference experiment leachates; however,
once the silicon releases are stabilized, the pH values also remain
stable over time. This acidification of the irradiated glass leachates
has already been observed in pure water in previous studies³⁶⁻³⁸
and under environmental conditions³⁹,⁴⁰. This could be due to
radiolysis effects, mainly air radiolysis, which is known to induce
nitric acid production in solution³⁶. This could occur when the
glass powders were transferred in water from the preparation cell
to the leaching cell before the closure of the titanium reactor and
the beginning of the experiment. However, in this pH range
(8.0–9.5), no effect of the pH value on the residual alteration rate
regime was expected, according to the literature⁴¹.

The concentration of silicon first increased in the leachates and then stabilized after less than 100 and 400 days for ²³⁸/²³⁹Pu-doped glasses and ²⁴⁴Cm-doped glass, respectively. The time necessary to reach this stabilization varied for the ²⁴⁴Cm-doped glass, which may be attributable to its radiological characteristics (Table 1), composition (very slightly different from the other glasses studied here, see Supplementary Table 1 in Supplementary Data), and the particle size distribution of the prepared powder, which could be slightly coarser. The last two factors could cause a slight delay in reaching the residual alteration regime. Thus, the residual alteration rate regime was considered to be reached after 100 and 400 days for the ²³⁸/²³⁹Pu-doped and ²⁴⁴Cm-doped glasses, respectively.

Once stabilized, the Si concentrations in solution were in
agreement with the saturation values of the SON68 glass with
respect to orthosilicic acid in solution (H₄SiO₄aa) proposed by Gin
et al.⁴¹ and deduced experimentally under the same leaching
conditions (gray area obtained for pH₉₀°C ranging from 7 to 9.5, in
Fig. 1). Furthermore, no colloids were observed by comparing the
filtered and ultrafiltered concentrations in solution (corresponding
to various Si concentrations, indicated by red arrows in some
samplings in Fig. 1). These data indicate that the altered layer
properties, in terms of solubility, are not significantly affected by
alpha self-irradiation. The present results are in agreement with
<table><caption>**Table 1.** Alpha dose rate and cumulative dose values of the doped glasses considered in this study.</caption><thead><tr><td>Alpha-doped glass</td><td>Alpha dose rate (Gy h<sup>-1</sup>)</td><td>Alpha cumulative dose at the beginning of leaching experiment</td></tr></thead><tbody><tr><td>0.4 wt% <sup>244</sup>Cm-doped glass</td><td>23,500</td><td>~ 3 × 10<sup>18</sup> α decay.g<sup>-1</sup>–35 MGy<sub>nucl</sub></td></tr><tr><td>Damaged <sup>238/239</sup>Pu-doped glass</td><td>3500</td><td>~ 1.2 × 10<sup>18</sup> α decay.g<sup>-1</sup>–14 MGy<sub>nucl</sub></td></tr><tr><td>Annealed <sup>238/239</sup>Pu-doped glass</td><td></td><td>&lt; 4 × 10<sup>16</sup> α decay.g<sup>-1</sup> (negligible)</td></tr><tr><td>Previously studied <sup>239</sup>Pu-doped glass<sup>20</sup></td><td>150</td><td>~ 3.7 × 10<sup>16</sup> α decay.g<sup>-1</sup> (negligible)</td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Fig. 2** Boron releases versus time for doped glasses leaching test and references experiments. Boron normalized mass loss (NL B) versus time for the <sup>238/239</sup>Pu- and <sup>239</sup>Pu-doped glass leaching tests (pink triangles and inversed purple triangles, respectively) and reference experiments ('ref. 1' with open circles and 'ref. 2' with open stars). The residual alteration rate value (Rr) associated with each experiment is given. Error bars are of ±10% and their calculation is detailed in the "Solution analyses" section.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 3** Boron normalized mass loss (NL B) versus time for annealed (pink triangles) and damaged (green diamonds) <sup>238/239</sup>Pu-doped glass leaching tests. Solid lines represent residual alteration rate. Error bars are of ±10% and their calculation is detailed in the "Solution analyses" section.</figcaption></figure>

those previously published on the <sup>239</sup>Pu-doped glass leachate and are corroborated by SEM/TEM observations of the microstructure of the alteration layer formed on the <sup>239</sup>Pu-doped glass, which was similar to that formed on non-radioactive samples<sup>20</sup>.

## Release of tracer elements

Glass alteration can be traced using soluble elements, especially boron, sodium, and lithium, released from the glass into the solution. Supplementary Figs 2–5 in Supplementary Data present the variations in normalized mass losses calculated from these releases for the three radioactive glasses and the reference experiment 'ref. 1' respectively. For each leaching test on actinide-doped glasses, the data show that, similar to non-radioactive experiments (Supplementary Fig. 5), congruence is conserved between the three tracer elements (B, Li, and Na), which indicates that the congruent behavior of these three elements is not affected by alpha self-irradiation. Moreover, the absence of colloids, including B, Li, or Na, was confirmed, and this was in agreement with similar observations in Si as previously described in the last paragraph. Thus, for simplicity, only boron releases are shown in Figs 2, 3, and 4 hereafter.

## Impact of the alpha dose rate

The boron releases as a function of time are presented in Fig. 2 for the actinide-doped glasses with very low alpha decay doses, that were subjected to only significant alpha dose rate parameters. The results for the annealed <sup>238/239</sup>Pu-doped glass were compared with those of the two reference experiments 'ref. 1' and 'ref. 2' and the previously studied <sup>239</sup>Pu-doped glass<sup>20</sup>.

The evolution of the boron releases was very similar for all the experiments: a rate drop was observed after 30–50 days and stationary releases over time attained after ~100 days. The boron normalized mass loss values were very close for the experiments on the annealed <sup>238/239</sup>Pu-doped glass and non-radioactive glass 'ref. 2' while, by comparison, the boron normalized mass loss values were initially higher for the experiments on the <sup>239</sup>Pu-doped glass and reference experiment 'ref. 1'. For glasses of the same composition (the two Pu-doped glasses and the 'ref. 2' glass), this difference could be linked to the absence of the precise knowledge of the specific surface area of the <sup>239</sup>Pu-doped glass powder. Indeed, the specific surface area was measured for powders from the non-radioactive glass and the annealed <sup>238/239</sup>Pu-doped glass (Table 4). The boron normalized mass loss values were, thus, accurately calculated. Otherwise, the <sup>239</sup>Pu-doped glass powder was obtained by re-sieving, between 63 and 125 µm in size, some powder batches prepared in the '1980s' in the Vulcain laboratories at CEA Marcoule. The particle size of this powder was therefore not measured but was assimilated to that of a non-radioactive powder prepared under the same conditions, as detailed in the previous study<sup>20</sup>. The resulting boron normalized mass loss values were thus less accurate than those of the other experiments. For the reference glass 'ref. 1' boron normalized mass loss values were higher than for the 'ref. 2' experiment, possibly due to the slightly different glass chemical composition (see Supplementary Table 1 in Supplementary Data), which could induce initial releases of tracer elements slightly increased.

It is worth noting that the maximum boron normalized mass loss values (around 1 g.m<sup>-2</sup>) correspond to an alteration layer thickness of 370 nm. The good correspondence between the thickness expected from the boron releases in solution and the one observed by TEM was previously confirmed on the <sup>239</sup>Pu-doped glass<sup>20</sup>. These TEM characterizations also confirmed that
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Time and square root of time dependence of boron releases for the doped glass leaching tests. Boron normalized mass loss (NL B) versus time (a) and square root of time (b), with pink triangles, green diamonds, and blue circles representing annealed <sup>238/239</sup>Pu-doped glass, damaged <sup>238/239</sup>Pu-doped glass, and <sup>244</sup>Cm-doped glass, respectively. In (a), the residual alteration rate values (Rr) associated with each experiment are given (with the relative Pearson's R-value) and expressed in g m<sup>-2</sup> d<sup>-1</sup>. In (b), NL B versus √t can be linearly fitted, either passing through the origin (dashed line) or with a positive y-intercept (solid lines). Error bars are of ±10% and their calculation is detailed in the "Solution analyses" section.</figcaption></figure>

the alteration layer was very adherent to the underlying unaltered glass, as commonly observed on non-radioactive glasses leached under similar conditions (high S/V ratio, no water renewal, for several years).

The residual alteration rates were calculated after stabilization of the tracer element releases in solution, that is, after ~100 days for these experiments. They are of the same order of magnitude as the considered experiment, ranging from (1.2 ± 0.3) × 10⁻⁴ g·m⁻²·d⁻¹ to (2.8 ± 0.4) × 10⁻⁴ g·m⁻²·d⁻¹. Moreover, the mean residual alteration rate values determined from the reference experiments on the non-radioactive glasses were close or even higher compared to those determined from the two Pu-doped glass leaching experiments. It should be noted that the residual alteration rate values with their associated uncertainties are very close in the powders for which the specific surface area has been measured (i.e., all but the <sup>239</sup>Pu-doped glass). The slightly lower residual alteration rate value calculated from the boron releases of the <sup>239</sup>Pu-doped glass is due to the uncertainty in the glass powder specific surface value and does not reflect any beneficial effect of the low alpha dose rate of this glass. An impact of radiolytic species created by alpha radiation of the leachate close to the glass surface was suspected in the literature to explain higher leaching rates<sup>34,37</sup>, when the leaching tests were performed in static conditions. By focusing on the results obtained here from glass powders for which the specific surface area has been measured (i.e., all, but the <sup>239</sup>Pu-doped glass), as the residual alteration rate of the most radioactive Pu-doped glass (3500 Gy·h⁻¹) was not higher than the residual alteration rates of the two reference experiments, it can be concluded that radiolytic phenomena at the glass-solution interface do not seem to modify the overall mechanisms and kinetics of glass alteration. Thus, even under these conditions where radiolytic species can accumulate with time, no impact of alpha radiolysis is expected. Therefore, this result is still in agreement with previous results on the residual

alteration rate regime<sup>20</sup>, which did not show any alpha dose rate value effect on the kinetics of <sup>239</sup>Pu-doped glass. We can therefore deduce that an alpha dose rate of up to 3500 Gy·h⁻¹ has no impact on the alteration kinetics of a R7T7-type glass in the residual rate regime. The present results allow us to further elaborate on the conclusion drawn in a previous paper<sup>20</sup> where the alpha dose rate had no impact on glass leaching behavior. The present study indicates that this is also the case for a higher dose rate value, and therefore, this behavior can be determined to cover a wide range of the alpha dose rate values expected under disposal conditions, even with early water ingress. The present results are also in agreement with those obtained for SON68 glasses altered under gamma irradiation<sup>21</sup>, which means that whatever the irradiation field and related energy deposition values at the glass-water interface<sup>42</sup>, the SON68-type glass remained relatively insensitive to the mechanism by which radiolysis of both the materials (glass, alteration layer) and the water was induced. The hypothesized protective role of pore water in limiting and/or healing the radiation defects created by radiolysis of the solid, as discussed in the literature concerning gamma external radiation of nanoporous glasses<sup>43,44</sup> and electron external irradiation of simplified glasses and corresponding alteration layers<sup>22</sup>, is invoked here to explain the absence of a dose rate impact on the alteration behavior of the glass in the residual alteration rate regime. In such conditions, the alteration layer remained very thin (<400 nm), adherent to the glass and thus, was protecting the glass from further alteration as observed on the non-radioactive glasses altered in the same conditions. Furthermore, these results did not support the 'dramatic' impact of the alpha dose rate reported in one study<sup>31</sup>, which could instead be related to the leaching conditions chosen in that study (aerated conditions, very low S/V, water renewal, and Teflon™ vessel) combined with a different glass composition (higher content of Na) rather than to the alpha dose rate value. Indeed,
the alteration layer developed on that $^{238}$Pu-doped glass was thick (10 µm), mechanically unstable, and probably composed of aluminum hydroxide, therefore very different from that expected in this study or the related one$^{20}$. No data concerning the pH was specified in that publication$^{31}$ and the glass leaching was monitored in solution by the Pu releases, which could not be considered as a tracer element without any experimental comparison with B and Na.

## Impact of the alpha decay dose

Since there was no discernable impact of the dose rate for the annealed $^{238/239}$Pu-doped glass, the results obtained on the annealed glass were, therefore, taken as a reference for comparison with the doped glasses with a significantly higher alpha decay dose. Indeed, this annealed glass has the same composition as the higher dose damaged glass, and all experiments were performed under exactly the same conditions (shielded cell, titanium leach reactors, radiolysis of glass-water interface).

The evolution of the boron normalized mass losses as a function of time is presented in Fig. 3 for both the damaged and annealed $^{238/239}$Pu-doped glasses. The two followed the same trend; the boron release initially greatly increased with time and then tended to stabilize after around 100 days, which indicates the beginning of the residual rate regime with the formation of a protective alteration layer$^1$. Once the residual alteration rate regime was reached, both the normalized mass loss and residual alteration rate values were significantly higher for the leaching experiment on the damaged glass than those of the annealed glass. This could be due to the alpha decay dose, as this is the only difference between these two experiments (same glass block before annealing and same leaching conditions).

However, in order to complete the reasoning, the results for the annealed and damaged $^{238/239}$Pu-doped glasses were compared to those of the $^{244}$Cm-doped glass (Fig. 4). The boron release values, as well as the residual alteration rates (slopes), remained higher for the two doped glasses having a significant alpha decay dose (Fig. 4a) than those of the annealed $^{238/239}$Pu-doped glass, considered as the reference here. The cumulative dose, prior to the leaching experiment, induced some modifications in the alteration behavior of the glass, which was manifested by an increase in the residual alteration rate values.

Furthermore, a slight decrease in the boron release rate could be observed for the $^{244}$Cm-doped glass (Fig. 4). The evolution of the normalized mass loss values was plotted as a function of the square root of time (Fig. 4b). For the three doped glasses considered here, the evolution is linear versus $\sqrt{t}$ either from the beginning of the leaching experiment ($^{244}$Cm-doped glass) or after a period of time ($^{238/239}$Pu-doped glasses). This time period lasted about 90 days (around 2750 s$^{1/2}$), which corresponds to the period of establishment of the residual alteration rate regime. The non-linearity with $\sqrt{t}$ within the first 90 days could thus be due to the predominance of another mechanism inducing a greater glass alteration for this $^{238/239}$Pu-doped glass (annealed and damaged) before the residual regime was established.

Then, as described in the literature, the evolution of the altered glass thickness is consistent with a residual alteration rate regime governed by a global diffusive mechanism$^{27,45-47}$. Thus, the apparent diffusion coefficient, $D_{app}$, was calculated from the altered thickness by applying Eq. (1), which was derived from Fick's second law by considering a semi-infinite medium. This apparent diffusion coefficient (in $m^2.s^{-1}$) is an overall value, which integrates different mechanisms such as the mobility and reactivity of the species. The $D_{app}$ values, derived from the plots displayed in Fig. 4b, are given in Table 2.

A $D_{app}$ value of $(1.5 \pm 0.2) \times 10^{-22} m^2.s^{-1}$, calculated from the leaching experiment on the annealed $^{238/239}$Pu-doped glass

<table><caption>**Table 2.** Apparent diffusion coefficient ($D_{app}$) calculated from the linear regression obtained by fitting NL B versus $\sqrt{t}$ in Fig. 4b.</caption><thead><tr><th></th><th>D<sub>app</sub> (m<sup>2</sup>.s<sup>-1</sup>)</th></tr></thead><tbody><tr><td><sup>244</sup>Cm-doped glass</td><td>(2.23 ± 0.05) × 10<sup>-21</sup></td></tr><tr><td>Damaged <sup>238/239</sup>Pu-doped glass</td><td>(2.70 ± 0.30) × 10<sup>-21</sup></td></tr><tr><td>Annealed <sup>238/239</sup>Pu-doped glass (acts as reference)</td><td>(1.5 ± 0.2) × 10<sup>-22</sup></td></tr></tbody></table>

considered as a reference, is of the same order of magnitude as that reported in the literature for non-radioactive SON68 powders $(1.4 \times 10^{-22} m^2.s^{-1})^{48}$ and 5–10 times lower than that obtained on the ISG glass$^{27,48}$. It is worth noting that this value of $(1.5 \pm 0.2) \times 10^{-22} m^2.s^{-1}$ is on the order of magnitude of the values expected for diffusion in solids but still ~3 orders of magnitude lower than the water diffusion values previously determined in borosilicate glasses and extrapolated at pH = 9$^5$, and ~50 times lower than the interdiffusion coefficient value calculated in the literature$^{49}$ and used to fit the experimental data on pre-altered SON68 glasses in saturated Si conditions$^{50}$ (90 °C, pH$_{90°C}$ = 9, S/V = 10 cm$^{-1}$). These observations indicate that various mechanisms are involved in this apparent diffusion coefficient.

$$ \text{Altered thickness} = \frac{NL_B}{\rho} = 2\sqrt{\frac{D_{app} \times t}{\pi}} \quad (1) $$

However, $D_{app}$ is 15 to 18 times higher in the damaged radioactive glasses compared to the reference value (Table 2). This increased value calculated for damaged actinides doped glasses is similar to that (×20) calculated from a leaching experiment performed on gold ion-damaged ISG glasses$^{27}$. Moreover, the nuclear doses induced by the alpha decays in the damaged $^{238/239}$Pu and $^{244}$Cm-doped glasses were 14 and 35 MGy$_{nucl}$, respectively (Table 1), at the beginning of the leaching experiments, and these values are similar to that needed to reach the plateau (i.e., around 14 MGy$_{nucl}$) in the alteration enhancement described for gold ion-damaged ISG glasses$^{24,27}$. These similarities suggest that the ballistic effects of the recoil nuclei are responsible for this increase in residual alteration rate observed for the damaged radioactive glasses. However, this increase was slightly lower in the case of the doped glasses than for externally irradiated simplified glasses. This could be explained either by differences in the glass composition and the higher sensitivity of simplified glass compositions toward irradiation$^{25}$ or by the partial annealing effect of the alpha particles in the case of alpha decay$^{51-54}$ as described in the literature for various borosilicate glass compositions. Moreover, the slightly lower boron release values seen in the case of the $^{244}$Cm-doped glass compared to the damaged $^{238/239}$Pu-doped glass were not inconsistent with a higher cumulative alpha decay dose value because the saturation value of 14 MGy$_{nucl}$$^{27}$ was reached for both doped glasses. This observation could, thus, be explained by either the slightly different glass composition or the fact that the specific surface area is not accurately known for the $^{244}$Cm-doped glass and could be slightly overestimated.

## Origin of the higher alteration kinetics of the damaged glasses

It should be noted that the alpha decay dose and the heavy ion irradiations$^{24,27}$ lead to similar changes in long-term glass alteration behavior. As proposed in the literature$^{24,27}$, this can be related to the changes in glass structure and properties. This link between glass structure and leaching behavior was also reported recently for glasses with simplified compositions whose structure was modified either by fast quenching$^{55,56}$ or by external irradiation by high-energy heavy ions$^{28}$. Some authors have
<figure><img src="image_7.png" /><figcaption>**Fig. 5** Schematic of glass alteration layer and three hypotheses detailed in text. The first hypothesis concerns the glass because it involves an initial water penetration. Then, the interdiffusion process induces the transformation of the glass in what is called the 'reactive interface'.</figcaption></figure>

highlighted the role of the B(III) content in glass to explain the decrease in glass chemical durability⁵⁵,⁵⁶. However, no significant impact of the nuclear dose was identified during the initial alteration rate regime of both actinide-doped borosilicate glasses³⁰,³²⁻³⁴ and glasses externally irradiated with heavy ions³², as measured by dynamic leaching tests. This initial alteration rate regime, occurring in very diluted solutions, is well known to be governed by the hydrolysis of the silicate network. These results, compared to the present one, suggest that, depending on the main mechanism(s) involved in a given medium (very diluted or more concentrated solutions), the chemical durability of a borosilicate glass is more or less impacted by its structure.

Based on these results, we developed three hypotheses to explain the impact of the nuclear contribution of the alpha decay dose on the long-term chemical durability of borosilicate glasses occurring during the residual alteration regime, with the formation of the protective alteration layer. Each hypothesis was related to a specific zone (glass, interfacial zone, and gel) and the corresponding possible mechanism(s) that could control the alteration during the residual alteration rate regime⁷,⁸,¹²⁻¹⁸. A schematic of the glass alteration profile is proposed in Fig. 5 to better illustrate the hypotheses described below.

1. Under irradiation, the glass is restructured, which led (among other effects) to a modification of the free volumes²³ and could induce an increase in the mobility of water in this irradiated glass close to the reactive interface¹⁸,⁵⁷ and consequently, an increase in the alteration kinetics of the interdiffusion mechanism. Depending on whether the most important mechanism affecting the glass structure is solely water penetration or the overall interdiffusion mechanism, including alkali migration and subsequent glass reorganization and boron release³, this first hypothesis could be closely linked to the second one.

2. The internal energy of the damaged glass is known to increase under irradiation; the alpha decay dose in an R7T7-type glass led to a stored energy of ~110 J.g⁻¹²³,⁵⁸ due to the formation of defects in its structure and the resulting increased disorder. A damaged glass could thus be seen as a glass having higher energy and entropy than pristine glass (i.e., non-irradiated glass). A possible higher reactivity of the glass at the reactive interface in contact with water could thus be hypothesized due to the higher energetic state of the damaged glass⁵⁹. Indeed, it could be schematically seen

<figure><img src="image_8.png" /><figcaption>**Fig. 6** Schematic of hypothesis 2, which describes energy of system (glass) as a function of reaction coordinate. This implies that the transition states of reactions occurring from the pristine and damaged glasses are similar.</figcaption></figure>

(Fig. 6) as a decrease in the activation energy necessary for the chemical reactions to hydrolyze the glass bonds⁶⁰ in this confined medium, assuming that the transition state remained similar in both cases (for pristine and damaged glasses). As no impact of the alpha decay dose was observed in the literature in a dilute medium governed by a total hydrolysis of the silicate network³⁰,³²⁻³⁴, we can consider that the confined nature of the medium where hydrolysis reactions could occur here is a crucial point. This hydrolysis could affect either weakly linked elements such as alkalis or former elements such as boron or silicon, which is preferentially hypothesized. For example, in the case of the damaged glass, despite a constant chemical composition, the B(III)/B(IV) ratio is higher than that of pristine glass²³,²⁴. Moreover, the preferential reactivity of boron (III) was recently shown and explained by the observation that B (III) acts as a Lewis acid, whereas water acts as a Lewis base⁶¹. This hypothesis of increased reactivity does not presume the partial/total ability of the involved hydrolyzed
reactions in the case of the silicate network hydrolysis, which could be partial and progressive (for example, one to three of the four Si–O bonds impacted in SiO$_4$) or total (all four Si–O bonds broken).

3. The gel could preserve the structure of the parent glass from which it was formed if its formation involved progressive hydrolysis-condensation reactions of the former elements⁸. In this case, one could think that the gel formed from a glass damaged by recoil nuclei could have different structures and thus different protective properties than that formed from a non-irradiated glass.

Some recent experimental and simulation studies can provide further information to discuss these three hypotheses. First, Jan et al.⁶² explored the direct impact of a modified glass structure (boron coordination, Qn species distribution, strength of Si–O–Si linkages) on leachability using a Monte-Carlo method coupled with molecular dynamics. They concluded that these modifications of the glass short-range order have a limited impact and indicated that, potentially, the medium-range order may more significantly affect the leachability of the glass⁶². They also used molecular dynamics to explore the mobility of water by building glasses containing 10 mol% H₂O and reproducing the recoil nuclei impact by ultrafast quenching⁶³. They qualitatively showed that the water mobility (i.e., diffusion) increased with the glass disorder.

Additional characterizations were recently carried out on irradiated and altered ISG glasses²⁷. First, the evolution of the gel porosity was observed by TEM in cryogenic mode as a function of the history of the parent glass²⁶. The authors showed no observable porosity in the gel obtained from a non-irradiated glass leached for 2 months and from an irradiated sample leached for a few days, whereas a porosity with average pore sizes of (1.5 ± 0.2) nm was observed on the pre-irradiated and 2 months leached ISG samples. The appearance of porosity during the alteration was recently reported in the literature⁸,¹⁹, and this would indicate a slow reorganization of the gel over time (several years), with a tendency to develop a more porous but also more protective structure with respect to glass dissolution by water. It seems possible, based on the current state of knowledge, that the microstructure of the irradiated glass gel would evolve faster than that of the non-irradiated glass, without making any link between this early "maturation" and the increased alteration kinetics of pre-irradiated glasses. In addition, the healing of pores was evidenced in simplified glass compositions by irradiating with low-energy heavy ions and monitoring by TEM in real time²⁹, demonstrating the ability of the gel to reorganize itself. This healing ability of the gel material could be hypothesized in the case of radioactive glasses that are subject to their own irradiation fields. However, as these radioactive glasses behave as externally irradiated glasses in terms of chemical dissolution²⁷, the low impact of gel porosity on glass chemical behavior could be hypothesized. Furthermore, ¹⁸O enriched water penetration in the alteration layer and the associated exchanges involving the pore water and the silanol entities were explored at room temperature²⁶. The ¹⁸O/¹⁶O exchanges were similar within the uncertainty limits in the major part of the gel, regardless of the parent glass history (irradiated or not), suggesting that the two gels had similar transport and exchange properties with respect to water regardless of their microstructures observed by TEM. However, in the reactive interface, the relative ¹⁸O/¹⁶O exchanges were significantly higher in the pre-irradiated glass compared to the non-irradiated one, suggesting that the previous irradiation modified the properties and reactivity of the glass and/or the reactive interface.

All these recent studies favor hypotheses 1 and 2 but not hypothesis 3. They also strongly suggest the need to refine our knowledge and description of the interfacial zone in order to identify the mechanism(s) responsible for the increase in glass alteration due to the nuclear impact of the alpha decay dose.

Finally, concerning the impact of the present results on long-term behavior prediction of nuclear glasses under disposal conditions, the level of alpha decay accumulated in a nuclear glass will exceed the dose required to significantly damage it (several 10¹⁸a.g⁻¹, i.e., 12 MGy$_{nucl}$) before water arrival, expected after several thousand years¹. Therefore, the glass damage has to be taken into account on the leaching studies in order to predict as accurately as possible the glass long-term chemical durability²⁴. As described in the present work, it is expected that the residual alteration rate value at a given time would be increased. But, the alteration rate evolution with square root of time indicates that the alteration rate still diminishes with time. Consequently, the impact of the alpha decay dose on the leaching rate will diminish with time under disposal conditions and the alteration rate of the damaged glass will become as low as the one observed on non-damaged glasses.

To conclude, actinide-doped SON68 glass powders were leached at 90 °C and at an S/V ratio ranging from 23 to 76 cm⁻¹ for up to 5 years to study the residual alteration rate regime. A part of an old ²³⁸/²³⁹Pu-doped glass block was annealed to fully recover the irradiation-induced damage for more than 30 years. Thus, the effect of the alpha dose rate parameter on leaching behavior could be dissociated from the effect of the alpha decay dose parameter on this old ²³⁸/²³⁹Pu-doped glass (14 MGy$_{nucl}$). A complementary study concerning a 0.4 wt% ²⁴⁴Cm-doped SON68 glass displaying a significant alpha decay dose (35 MGy$_{nucl}$) and two reference experiments on non-radioactive glasses were also performed.

The annealed ²³⁸/²³⁹Pu-doped glass showed similar behavior in water to the corresponding non-radioactive glasses and also similarly to a previously studied ²³⁹Pu-doped glass²⁰. Thus, the residual alteration rate value is not affected by the alpha dose rate, and this conclusion is demonstrated for the wide range of dose rate values expected under disposal conditions (≤3500 Gy·h⁻¹), even in the eventuality of an early water ingress. This result is also in agreement with the one obtained for SON68 glasses altered under gamma irradiation (up to 10,000 Gy·h⁻¹), which means that, irrespective of the irradiation fields and related energy depositions at the glass-water interface, SON68-type glass remains relatively insensitive to a medium inducing radiolysis of the materials (glass, alteration layer) and water. The protective role of pore water can be hypothesized to explain this result.

The residual alteration rate of the damaged ²³⁸/²³⁹Pu-doped glass was clearly enhanced compared to that of the annealed glass (i.e., having a very low level of damage). This was also observed for the damaged ²⁴⁴Cm-doped glass and is in agreement with that of ISG glasses externally irradiated with gold ions²⁴,²⁷, indicating that the nuclear effects of the recoil nuclei are responsible for this increase in the residual alteration rate due to the alpha decay accumulation. The present results and recent studies from the literature²⁸,⁵⁵,⁵⁶ demonstrate a link between the glass structure and its leaching behavior. Although no mechanisms are clearly proposed to explain such a link, the role of the reactive interface is strongly suspected.

Further studies and characterizations should focus on this layer to provide a better understanding of the precise mechanisms underlying it. Indeed, the solid characterization of the actinide-doped altered glasses by SEM and TEM, which compares the microstructures of the alteration layers formed as a function of the parent glass history (with or without an alpha decay dose), will provide valuable information to complete the description of the alpha decay dose impact and the suspected involved mechanism(s). The consideration of damaged glasses provides an opportunity to focus on different glass structures with the same chemical composition and can greatly improve our understanding of long-term leaching mechanisms. The healing
ability of the gel formed under self-irradiation will also be investigated. On externally irradiated simplified glasses, fully leached fine glass powders can be characterized by CPMAS NMR in order to complete a possible correlation between elements in the glass ($^{11}$B, $^{29}$Si) and water species ($^{1}$H)$^{64}$ in the alteration layer and specifically at the reactive interface.

# METHODS

## Glass materials

Two Pu-doped glasses were fabricated in a glove box in the Vulcain facility at CEA Marcoule in 1985 by means of the melting pot process$^{35}$. A nitric acid solution containing the radionuclides was first added to a simulated inactive fission product solution and then mixed with the glass frit prepared in the form of sludge. No reducing components were added. The mixture was then heated at 1150 °C for 3 h. The glass was then annealed for 1 h at 525 °C and cooled at a rate lower than 50 °C.h$^{-1}$. Under these fabrication conditions and with a loading factor of 0.85 wt%, Pu was expected to be homogeneously dissolved in the glass matrix$^{65}$.

The chemical compositions of the Pu-doped glasses were determined after hot acidic dissolution of glass specimens and ICP-AES analyses of the non-radioactive elements and were consistent with the calculated compositions. They are detailed in Supplementary Data (Supplementary Table 1). The Pu concentrations in the glasses were determined from the thermal power measurement of glass samples by isothermal calorimetry, and the Pu isotopic composition was determined. The two Pu-doped glasses both included ~0.85 wt% of $PuO_2$ but differed in their isotopic distribution; they were thus called $^{239}$Pu and $^{238/239}$Pu-doped glasses because of their isotopic composition, as detailed in Table 3. The glass that was more enriched in $^{239}$Pu has already been extensively studied$^{20,42}$ and is used for comparison with the results presented here. The second glass was mainly enriched in $^{238}$Pu compared to $^{239}$Pu to provide a greater alpha dose rate (Tables 1 and 3, in bold).

0.4 wt% $^{244}$Cm$O_2$-doped glass (called $^{244}$Cm-doped glass) was fabricated in the high-level waste laboratory (DHA) of the Atalante facility at CEA Marcoule in 2004. The fabrication and initial characterization of this glass have already been described in detail$^{66}$. The glass sample was produced by melting oxides, carbonates, and nitrates. The first step, prior to initiating each glass melt, consisted of dissolving curium oxide in a 3 mol.L$^{-1}$ nitric acid solution and measuring the Cm concentration by isothermal calorimetry in order to adjust the solution volume to obtain the desired concentration in the glass. The isotopic composition of the curium used for glass fabrication is shown in Table 3. The content of the main isotope, $^{244}$Cm, is written in bold. The curium-doped nitric acid solution was then added to the oxides in a platinum crucible. This method was selected rather than a dry route to promote curium solubilization in the glass. The mixture was melted for 3 h at 1200 °C without stirring, then quenched to a temperature slightly above the glass transition point (600 °C), from which it was cooled more slowly to release the residual stresses from the glass. The

theoretical chemical composition of this $^{244}$Cm-doped glass is detailed in Supplementary Data (Supplementary Table 1).

The aim of this study was to assess the impact of alpha self-irradiation on long-term aqueous behavior by dissociating the impact of the alpha dose rate and the impact of the alpha decay dose.

The radioactive glasses studied here, doped with short-lived actinides ($^{244}$Cm, $^{238}$Pu), were subjected to self-radiation, inducing varying alpha dose rate values. Consequently, after several years, they reached alpha decay doses equivalent to those sustained by the R7T7-type glass several thousand years after fabrication. The alpha dose rate values and the alpha cumulative dose of these glasses are presented in Table 1. The characteristics of these two glasses were compared to those of the previously studied $^{239}$Pu-doped glass$^{20,42}$. The alpha cumulative dose values of the annealed $^{238/239}$Pu-doped glass and of the $^{239}$Pu-doped glass previously studied are considered as negligible because they induced <5% modifications to the solid structure and properties, compared to the modifications expected at the saturation dose value$^{23}$.

The alpha dose rate values of the two Pu-doped glasses (i.e., $^{238/239}$Pu and $^{239}$Pu-doped glasses) covered a wide range of the alpha dose rate values expected for a typical R7T7 glass over thousands of years$^{21}$. Moreover, it was previously demonstrated by calculation$^{42}$ that the mean dose rates calculated from the glass composition were a good approximation of the local dose rate values both in the glass side and in the water side in contact with the glass powder.

Concerning the alpha decay dose values, at the beginning of the leaching tests, the $^{238/239}$Pu-doped and $^{244}$Cm-doped glasses were already highly damaged regarding the doses needed to induce the evolution of the glass structure and properties, whose evolution becomes significant after around 12 MGYnucl and saturates at around 40 MGYnucl$^{23}$.

In order to dissociate the alpha dose rate parameter from the alpha decay dose parameter, the $^{238/239}$Pu-doped glass block was cut into two equal parts, and two monoliths were also separated from this block and polished to obtain a surface roughness of <1 µm on one side (see Fig. 7). One part of the cut block and one polished monolith were annealed for 3 h at 545 °C, that is, 30 °C higher than the Tg temperature of this glass. Then, the cooling was slow ($15$ °C.h$^{-1}$ between 545 °C and 350 °C) in order to allow the glass to relax and recover its initial structure immediately after fabrication (see Supplementary Fig. 1 in Supplementary Data). This annealing efficiency was checked by measuring the hardness of both the annealed and damaged polished monoliths. An increase of 29% in the hardness value was measured after annealing, which was in agreement with the expected hardness difference (30%) between the pristine and damaged glasses at a saturation value given in the literature$^{23}$.

To obtain a high glass surface to water volume ratio and then quickly reach the residual alteration regime during the leaching experiment, glass samples were ground and sieved to recover the 63–125 µm size fraction. Then, they were washed in water to remove the finest particles. The specific surface area of the annealed and damaged $^{238/239}$Pu-doped glass powder samples was measured by krypton adsorption using the BET method$^{67}$, and the values are given in Table 4. The $^{244}$Cm-doped glass powder was too highly radioactive to be manipulated in glove boxes and to perform such a specific surface area measurement. Thus, the specific

<table><caption>**Table 3.** Isotopic compositions of the doped glasses at the beginning of the leaching tests (normalized to the Pu and Cm contents).</caption><thead><tr><th>wt%</th><th><sup>239</sup>Pu-doped glass</th><th><sup>238/239</sup>Pu-doped glass</th><th><sup>244</sup>Cm-doped glass</th></tr></thead><tbody><tr><td><sup>238</sup>Pu</td><td>0.1</td><td><strong>21.8</strong></td><td></td></tr><tr><td><sup>239</sup>Pu</td><td><strong>75.3</strong></td><td><strong>47.8</strong></td><td>0.2</td></tr><tr><td><sup>240</sup>Pu</td><td>21.7</td><td>18.7</td><td>36.7</td></tr><tr><td><sup>241</sup>Pu</td><td>0.7</td><td>1.0</td><td></td></tr><tr><td><sup>242</sup>Pu</td><td>0.6</td><td>5.0</td><td></td></tr><tr><td><sup>241</sup>Am</td><td>1.5</td><td>5.6</td><td></td></tr><tr><td><sup>243</sup>Cm</td><td></td><td></td><td>0.6</td></tr><tr><td><sup>244</sup>Cm</td><td></td><td></td><td><strong>45.3</strong></td></tr><tr><td><sup>245</sup>Cm</td><td></td><td></td><td>8.1</td></tr><tr><td><sup>246</sup>Cm</td><td></td><td></td><td>8.3</td></tr><tr><td><sup>247</sup>Cm</td><td></td><td></td><td>0.4</td></tr><tr><td><sup>248</sup>Cm</td><td></td><td></td><td>0.4</td></tr></tbody></table>

<figure><img src="image_9.png" /><figcaption>**Fig. 7** Schematic view of the annealing of a part of the $^{238/239}$Pu-doped glass. The monoliths were polished to micron range roughness on one face.</figcaption></figure>
<table><caption>**Table 4.** Glass powder samples and leaching conditions descriptions.</caption><thead><tr><th></th><th>0.4 wt% <sup>244</sup>CmO<sub>2</sub>-<br>doped glass</th><th>Damaged <sup>238/239</sup>Pu-<br>doped glass</th><th>Annealed <sup>238/239</sup>Pu-<br>doped glass</th><th><sup>239</sup>Pu-doped<br>glass<sup>20</sup></th><th>'ref. ¹' experiment</th><th>'ref. ²' experiment<sup>41</sup></th></tr></thead><tbody><tr><td>S<sub>spe</sub> (cm<sup>2</sup> g<sup>-1</sup>)</td><td>645 ± 50 (theoretical)</td><td>720 ± 68 (measured)</td><td>2340 ± 230<br>(measured)</td><td>645 ± 50<br>(theoretical)</td><td>885 ± 62<br>(measured)</td><td>860 ± 40<br>(measured)</td></tr><tr><td>Initial S/V<br>(cm<sup>-1</sup>)</td><td>~26</td><td>23 ± 3</td><td>76 ± 7</td><td>~20</td><td>38 ± 3</td><td>50 ± 3</td></tr><tr><td>m<sub>glass powder</sub> (g)</td><td>11.8</td><td>9.6</td><td>9.6</td><td>9.4</td><td>13.2</td><td>4.651</td></tr><tr><td>V (L)</td><td>0.297</td><td>0.298</td><td>0.307</td><td>0.305</td><td>0.310</td><td>0.080</td></tr><tr><td>Leaching<br>duration</td><td>>5 years</td><td>~3 years</td><td>>3.5 years</td><td>3.2 years</td><td>2.7 years</td><td>2.2 years</td></tr></tbody></table>

surface area of this <sup>244</sup>Cm-doped glass powder sample was estimated at about 645 ± 50 cm<sup>2</sup>·g<sup>-1</sup> based on the tests performed with a ground and sieved inactive SON68 glass powder (detailed in Rolland et al.<sup>20</sup>). To ensure this specific surface area value, the grain size was qualitatively checked by timing the settling rate of this <sup>244</sup>Cm-doped glass powder in water. Using Stokes' law, it was determined that more than 80% of the glass grains were larger than 60 ± 7 µm. In Table 4, uncertainties in the initial S/V ratio are specified only when the S<sub>spe</sub> has been experimentally measured.

It is worth noting that the higher than expected specific surface value, S<sub>spe</sub>, measured for the annealed glass powder sample was probably due to less efficient washing of the powder in water. The smallest particle size was estimated to be <1.4 µm by using Stokes' law in water and measuring the maximum settling rate (>7 h in 5 cm water).

## Leaching device and conditions

Radioactive glass leaching tests were performed in a shielded cell in the Atalante facility at CEA Marcoule in titanium reactors. Before entering the shielded cell, each titanium reactor was annealed for 3 h at 450 °C to form a passivating TiO<sub>2</sub> film to minimize the possible interactions between the species generated by water radiolysis and the experimental setup. In the literature, it was shown that finely divided titanium oxides were capable of decomposing H<sub>2</sub>O<sub>2</sub> through redox reactions<sup>68</sup>. In the present case, the titanium surfaces in the reactors were polished and thus, we assume that these surfaces covered by TiO<sub>2</sub> were not sufficient to significantly affect the H<sub>2</sub>O<sub>2</sub> concentration produced by alpha radiolysis under anoxic conditions in the present systems. Indeed, this was experimentally demonstrated in a previous study with small titanium containers under gamma radiation<sup>21</sup>, where H<sub>2</sub>O<sub>2</sub> experimental concentrations were reproduced by calculations with CHEMSIMUL<sup>69</sup> for a homogeneous closed system of pure water.

For each experiment, about 10 g of glass powder was first mixed with 30 mL of ultrapure water to allow it to be easily transferred into the leaching reactor, and the leachate volume was adjusted to ~300 mL by adding ultrapure water. The leaching tests were thus carried out under static conditions at an initial glass-surface-area-to-solution volume (S/V) ratio of ~23 to 76 cm<sup>-1</sup>, as described in Table 4.

After reactor closure, the solution was degassed by argon bubbling for 3 h (argon purity = 99.9999%). Then, an argon overpressure of 4 bars was applied to prevent any ingress of air into the reactor and to allow leachate sampling at regular intervals. Finally, the temperature was adjusted to 90 °C. The total duration of each radioactive glass leaching test was at least 3 years.

Two reference experiments, called 'ref. 1' and 'ref. 2', were performed under similar conditions (90 °C, S/V values) and were used to compare the present results obtained for Pu and Cm-doped glasses to those obtained for non-radioactive glasses. The 'ref. 1' experiment was performed in a titanium reactor similar to the one used for radioactive glasses, but the SON68 glass had a slightly different composition (see Supplementary Table 1 in Supplementary Data), whereas the 'ref. 2' experiment was performed on the same glass composition compared to the Pu-doped glasses, but in a 120 mL PFA Savillex™ reactor. PFA reactors are commonly used for leaching experiments in a non-radioactive environment. This second reference experiment, already described in detail<sup>41</sup>, also provided a reference in terms of leaching devices.

## Solution analyses

The leachate pH was measured in an aliquot at room temperature using either a combination gel-filled pH electrode (Checker) in the hot cell or an ISFET probe (Metrohm) if the measurements were performed in glove boxes. The measurement uncertainty was ±0.1 pH unit. The pH at 90 °C was calculated to be about 0.4 pH units lower than the pH measured at 25 °C, according to both CHESS calculations<sup>70</sup> and experimental measurements on non-radioactive solutions<sup>71</sup>. The uncertainty of the calculated pH<sub>90°C</sub> was ±0.3 pH unit.

For each experiment, most of the leachate samples were analyzed directly, except for the annealed <sup>238/239</sup>Pu-doped glass where the leachate was systematically filtered (0.45 µm) to prevent the presence of small glass grains. However, for all experiments, some samples were filtered (0.45 µm) and ultrafiltered (10,000 Dalton) to check for colloids, as previously described by Rolland et al.<sup>21</sup>. Both fractions were then acidified with nitric acid to reach a pH ranging from 0.5 to 1. Cation analyses were performed using a Horiba Jobin Yvon JY46 ICP-AES system. The analytical uncertainty was ±3% for concentrations above 10 mg·L<sup>-1</sup>, ±5% between 1 and 10 mg·L<sup>-1</sup>, and ±20% for concentrations below 1 mg·L<sup>-1</sup>. The releases in solution were expressed in terms of normalized mass losses (NL in g·m<sup>-2</sup>) for element *i* and were calculated as follows for a leachate sample taken at time *t*:

$$NL_i^t = \frac{m_i^t}{x_i \times S} \qquad (2)$$

where

$$m_i^t = C_i^t \times \left( V - \sum_{j=1}^{t-1} v^j \right) + \sum_{j=1}^{t-1} C_j^t \times v^j \qquad (3)$$

$m_i^t$ represents the amount of element *i* released in solution at time *t* (in g), $C_i^t$ is the measured concentration of element *i* at time *t* (in g·L<sup>-1</sup>), $C_j^t$ is the concentration measured in the previous leachate sample ($t-1$) for which the sampled volume was $v^{t-1}$, $x_i$ is the mass fraction of element *i* in the glass, *S* is the glass surface area (corresponding to S<sub>spe</sub> × m<sub>glass powder</sub> in m<sup>2</sup>), and *V* is the initial volume of the solution (L). For tracer elements such as B, Na, and Li, this parameter (NL) represents the amount of altered glass. The relative uncertainty calculated on NL is ~10%, with major contributions coming from the specific area measurements/estimations (see Table 4), cation analyses (see above), and the glass formula (3% for B, Na, Li, and Si). Error bars are specified in each figure.

When the glass grain size was much greater than the glass altered thickness, this altered thickness can be directly calculated from the normalized mass losses of boron by dividing it by the density of the glass ρ (g·cm<sup>-3</sup>). This density is taken as 2.75 for all the glasses studied here.

$$\text{Altered thickness} = \frac{(NL_B)_t}{\rho} (\mu\text{m}) \qquad (4)$$

The residual alteration rate for element *i* (Rr<sub>i</sub>) corresponds to the slope of the NL<sub>i</sub> evolution versus time (in days) in a given time range, and can be described by Eq. (5).

$$Rr_i = \frac{d(NL_i)}{dt} (\text{g.m}^{-2}.\text{d}^{-1}) \qquad (5)$$

It was calculated using the linear fit function in the Origin software using 'direct weighting' mode. This mode takes into account the entire uncertainty dataset with the same 'weight'.
DATA AVAILABILITY

The data that support the findings of this study are available from the corresponding author upon reasonable request.

Received: 17 November 2020; Accepted: 7 June 2021;
Published online: 07 July 2021

REFERENCES

1. Gin, S., Jollivet, P., Tribet, M., Peulet, S. & Schuller, S. Radionuclides containment in nuclear glasses: an overview. *Radiochim. Acta* **105**, 927–959 (2017).

2. Doremus, R. H. Interdiffusion of hydrogen and alkali ions in a glass surface. *J. Non Cryst. Solids* **19**, 137–144 (1975).

3. Bouyer, F., Geneste, G. & Gin, S. Hydrogen-sodium interdiffusion in borosilicate glasses investigated from first principles. *J. Non Cryst. Solids* **352**, 3147–3152 (2006).

4. Rébiscoul, D., Bruguier, F., Magnin, V. & Gin, S. Impact of soda-lime borosilicate glass composition on water penetration and water structure at the first time of alteration. *J. Non Cryst. Solids* **358**, 2951–2960 (2012).

5. Rébiscoul, D. et al. Water penetration mechanisms in nuclear glasses by X-ray and neutron reflectometry. *J. Non Cryst. Solids* **353**, 2221–2230 (2007).

6. Advocat, T., Crovisier, J. L., Vernaz, E., Ehret, G. & Charpentier, H. Hydrolysis of R7T7 nuclear waste glass in dilute media: mechanisms and rates as a function of pH. *MRS Proc.* **212**, 57–64 (1991).

7. Hellmann, R. et al. Nanometre-scale evidence for interfacial dissolution-reprecipitation control of silicate glass corrosion. *Nat. Mater.* **14**, 307–311 (2015).

8. Gin, S. et al. Dynamics of self-reorganization explains passivation of silicate glasses. *Nat. Commun.* **9**, 2169 (2018).

9. Valle, N. et al. Elemental and isotopic (Si-29 and O-18) tracing of glass alteration mechanisms. *Geochim. Cosmochim. Acta* **74**, 3412–3431 (2010).

10. Mercado-Depierre, S., Fournier, M., Gin, S. & Angeli, F. Influence of zeolite precipitation on borosilicate glass alteration under hyperalkaline conditions. *J. Nucl. Mater.* **491**, 67–82 (2017).

11. Gin, S. et al. An international initiative on long-term behavior of high-level nuclear waste glass. *Mater. Today* **16**, 243–248 (2013).

12. Gin, S. et al. Origin and consequences of silicate glass passivation by surface layers. *Nat. Commun.* **6**, 8 (2015).

13. Geisler, T. et al. Aqueous corrosion of borosilicate glass under acidic conditions: a new corrosion mechanism. *J. Non Cryst. Solids* **356**, 1458–1465 (2010).

14. Geisler, T. et al. The mechanism of borosilicate glass corrosion revisited. *Geochim. Cosmochim. Acta* **158**, 112–129 (2015).

15. Lenting, C. et al. Towards a unifying mechanistic model for silicate glass corrosion. *npj Mater. Degrad.* **2**, 28 (2018).

16. Cailleteau, C. et al. Insight into silicate-glass corrosion mechanisms. *Nat. Mater.* **7**, 978–983 (2008).

17. Gin, S. et al. Atom-Probe Tomography, TEM and ToF-SIMS study of borosilicate glass alteration rim: a multiscale approach to investigating rate-limiting mechanisms. *Geochim. Cosmochim. Acta* **202**, 57–76 (2017).

18. Bunker, B. C. Molecular mechanisms for corrosion of silica and silicate glasses. *J. Non Cryst. Solids* **179**, 300–308 (1994).

19. Ngo, D. et al. Spectroscopic ellipsometry study of thickness and porosity of the alteration layer formed on international simple glass surface in aqueous corrosion conditions. *npj Mater. Degrad.* **2**, 20 (2018).

20. Rolland, S. et al. Tc-99- and Pu-239-doped glass leaching experiments: residual alteration rate and radionuclide behavior. *Int. J. Appl. Glass Sci.* **4**, 295–306 (2013).

21. Rolland, S. et al. Influence of gamma irradiation effects on the residual alteration rate of the French SON68 nuclear glass. *J. Nucl. Mater.* **433**, 382–389 (2013).

22. Mougnaud, S. et al. Effect of low dose electron beam irradiation on the alteration layer formed during nuclear glass leaching. *J. Nucl. Mater.* **482**, 53–62 (2016).

23. Peulet, S., Delaye, J. M. & Jegou, C. Specific outcomes of the research on the radiation stability of the French nuclear glass towards alpha decay accumulation. *J. Nucl. Mater.* **444**, 76–91 (2014).

24. Peulet, S., Tribet, M., Mougnaud, S., Miro, S. & Jegou, C. Radiation effects in ISG glass: from structural changes to long term aqueous behavior. *npj Mater. Degrad.* **2**, 23 (2018).

25. Mir, A. H. & Peulet, S. Using external ion irradiations for simulating self-irradiation damage in nuclear waste glasses: State of the art, recommendations and, prospects. *J. Nucl. Mater.* **539**, https://doi.org/10.1016/j.jnucmat.2020.152246 (2020).

26. Tribet, M. et al. New insights about the importance of the alteration layer/glass interface. *J. Phys. Chem. C.* **124**, 10032–10044 (2020).

27. Mougnaud, S. et al. Heavy ion radiation ageing impact on long-term glass alteration behavior. *J. Nucl. Mater.* **510**, 168–177 (2018).

28. Lonartz, M. I. et al. The effect of heavy ion irradiation on the forward dissolution rate of borosilicate glasses studied in situ and real time by fluid-cell Raman spectroscopy. *Materials* **12**, 1480 (2019).

29. Mir, A. H. et al. Effect of decades of corrosion on the microstructure of altered glasses and their radiation stability. *npj Mater. Degrad.* **4**, 11 (2020).

30. Tribet, M. et al. Irradiation impact on the leaching behavior of HLW glasses. *Procedia Mater. Sci.* **7**, 209–215 (2014).

31. Zubekhina, B. Y. et al. Chemical alteration of 238Pu-loaded borosilicate glass under saturated leaching conditions. *Radiochim. Acta* **108**, https://doi.org/10.1515/ract-2018-3097 (2019).

32. Peulet, S. et al. Effect of alpha radiation on the leaching behaviour of nuclear glass. *J. Nucl. Mater.* **362**, 474–479 (2007).

33. Fillet, S. *Mécanismes de corrosion et comportement des actinides dans le verre nucléaire "R7T7".* Thèse de l'Université Montpellier II, France (1987).

34. Wellman, D. M., Icenhower, J. P. & Weber, W. J. Elemental dissolution study of Pu-bearing borosilicate glasses. *J. Nucl. Mater.* **340**, 149–162 (2005).

35. Fillet, S., Nogues, J. L., Vernaz, E. & Jacquet-Francillon, N. Leaching of actinides from the French LWR reference glass. *MRS Proc.* **50**, 211–218 (1985).

36. Advocat, T., Jollivet, P., Crovisier, J. L. & del Nero, M. Long-term alteration mechanisms in water for SON68 radioactive borosilicate glass. *J. Nucl. Mater.* **298**, 55–62 (2001).

37. Weber, W. J., Wald, J. W. & McVay, G. L. Effects of alpha-radiation on the leaching of a nuclear waste glass. *J. Am. Ceram. Soc.* **68**, C253–C255 (1985).

38. Schuman, R. P. Leach testing of Idaho-National-Engineering-Laboratory waste forms in a gamma-field. *Nucl. Technol.* **65**, 422–431 (1984).

39. Lemmens, K. & Van Iseghem, P. The effect of gamma radiation on the dissolution of high-level waste glass in Boom Clay. *MRS Proc.* **663**, 227–235 (2001).

40. Lemmens, K., Van Iseghem, P. & Wang, J. The leaching of Pu, Am, Np and Tc from high-level waste glasses in clay media. *MRS Proc.* **294**, 147–154 (1993).

41. Gin, S. & Mestre, J. P. SON 68 nuclear glass alteration kinetics between pH 7 and pH 11.5. *J. Nucl. Mater.* **295**, 83–96 (2001).

42. Mougnaud, S., Tribet, M., Rolland, S., Renault, J. P. & Jegou, C. Determination of alpha dose rate profile at the HLW nuclear glass/water interface. *J. Nucl. Mater.* **462**, 258–267 (2015).

43. Le Caer, S., Renault, J. P. & Mialocq, J. C. Hydrogen peroxide formation in the radiolysis of hydrated nanoporous glasses: a low and high dose rate study. *Chem. Phys. Lett.* **450**, 91–95 (2007).

44. Klett, C. et al. H-2 production through oxide irradiation: comparison of gamma rays and vacuum ultraviolet excitation. *Int. J. Hydrog. Energy* **38**, 3889–3897 (2013).

45. Chave, T., Frugier, P., Ayral, A. & Gin, S. Solid state diffusion during nuclear glass residual alteration in solution. *J. Nucl. Mater.* **362**, 466–473 (2007).

46. Frugier, P. et al. SON68 nuclear glass dissolution kinetics: current state of knowledge and basis of the new GRAAL model. *J. Nucl. Mater.* **380**, 8–21 (2008).

47. Aréna, H., Rébiscoul, D., Garcés, E. & Godon, N. Comparative effect of alkaline elements and calcium on alteration of International Simple Glass. *npj Mater. Degrad.* **3**, 10 (2019).

48. Gin, S., Beaudoux, X., Angeli, F., Jegou, C. & Godon, N. Effect of composition on the short-term and long-term dissolution rates of ten borosilicate glasses of increasing complexity from 3 to 30 oxides. *J. Non Cryst. Solids* **358**, 2559–2570 (2012).

49. Grambow, B. & Müller, R. First-order dissolution rate law and the role of surface layers in glass performance assessment. *J. Nucl. Mater.* **298**, 112–124 (2001).

50. Jegou, C. *Mise en évidence expérimentale des mécanismes limitant l'altération du verre R7T7 en milieu aqueux. Critique et proposition d'évolution du formalisme cinétique.* Thèse de l'Université Montpellier II, France (1998).

51. Charpentier, T. et al. Self-healing capacity of nuclear glass observed by NMR spectroscopy. *Sci. Rep.* **6**, 25499 (2016).

52. Mir, A. H. *Radiation Effects on Oxide Glasses: Importance of Energy Deposition and Relaxation Processes*. Thèse de l'Université de Caen, France (2015).

53. Mir, A. H. et al. Defect recovery and damage reduction in borosilicate glasses under double ion beam irradiation. *Epl* **112**, 36002 (2015).

54. Karakurt, G. et al. Understanding of the mechanical and structural changes induced by alpha particles and heavy ions in the French simulated nuclear waste glass. *J. Nucl. Mater.* **475**, 243–254 (2016).

55. Angeli, F. et al. Effect of thermally induced structural disorder on the chemical durability of International Simple Glass. *npj Mater. Degrad.* **2**, 31 (2018).

56. Stone-Weiss, N. et al. Understanding the structural drivers governing glass-water interactions in borosilicate based model bioactive glasses. *Acta Biomater.* **65**, 436–449 (2018).

57. Wang, Q. Y. et al. Effects of 80 KeV proton radiation on the optical properties and microstructure of type-GG17 glass as rubidium lamp envelope. *Nucl. Instr. Meth. Phys. Res. B* **267**, 2489–2494 (2009).

58. Maugeri, E. A. et al. Calorimetric study of glass structure modification induced by α decay. *J. Am. Ceram. Soc.* **95**, 2869–2875 (2012).
59. Ngo, D. et al. Hydrogen bonding interactions of H₂O and SiOH on a bor-oaluminosilicate glass corroded in aqueous solution. *npj Mater. Degrad.* **4**, 1 (2020).

60. Zapol, P., He, H., Kwon, K. D. & Criscenti, L. J. First-principles study of hydrolysis reaction barriers in a sodium borosilicate glass. *Int. J. Appl. Glass Sci.* **4**, 395-407 (2013).

61. Kapoor, S. et al. Structural and chemical approach toward understanding the aqueous corrosion of sodium aluminoborate glasses. *J. Phys. Chem. B* **122**, 10913-10927 (2018).

62. Jan, A., Delaye, J.-M., Gin, S. & Kerisit, S. Monte Carlo simulation of the corrosion of irradiated simplified nuclear waste glasses. *J. Non Cryst. Solids* **519**, 119449 (2019).

63. Jan, A. *Impact des effets balistiques sur les propriétés de la couche de gel des verres nucléaires simplifiés: une approche de simulation Monte Carlo*. Thèse de l'Université de Montpellier, France (2020).

64. Brunet, F., Charpentier, T., Le Caer, S. & Renault, J. P. Solid-state NMR characterization of a controlled-pore glass and of the effects of electron irradiation. *Solid State Nucl. Magn. Reson.* **33**, 1-11 (2008).

65. Lopez, C. *Solubilité des actinides et de leurs simulants dans les verres nucléaires. Limite d'incorporation et compréhension des mécanismes*. Thèse de l'Université Paris XI, France (2003).

66. Peuget, S. et al. Irradiation stability of R7T7-type borosilicate glass. *J. Nucl. Mater.* **354**, 1-13 (2006).

67. Brunauer, S., Emmett, P. H. & Teller, E. Adsorption of gases in multimolecular layers. *J. Am. Chem. Soc.* **60**, 309-319 (1938).

68. Suh, M., Bagus, P. S., Pak, S., Rosynek, M. P. & Lunsford, J. H. Reactions of hydroxyl radicals on titania, silica, alumina, and gold surfaces. *J. Phys. Chem. B* **104**, 2736-2742 (2000).

69. Kirkegaard, P. & Ejegbakke, E. CHEMSIMUL: a Simulator for Chemical Kinetics (2000).

70. Van der Lee, J. & de Windt, L. CHESS Tutorial and Cookbook. Updated for version 3.0, Users manual N°LHM/RD/02/13 (2002).

71. Advocat, T. *Les mécanismes de corrosion en phase aqueuse du verre nucléaire R7T7. Approche expérimentale. Essai de modélisation thermodynamique et cinétique*. Thèse de l'université Louis Pasteur, Strasbourg, France (1991).

## ACKNOWLEDGEMENTS

The authors express their gratitude to CEA, Orano, and EDF (France) for funding the research within the framework of the VESTALE project. They also wish to thank all the LMPA staff for their contributions to the experiments on radioactive glasses; J. Maurin, R. Burgaud, and J. Muller for the specific surface area measurements of radioactive glass powders; and S. Gin for sharing the results from the reference experiment 'ref. 2'.

## AUTHOR CONTRIBUTIONS

M.T. supervised the entire study and wrote the paper. C.M. and V.B. were responsible for the experimental analysis of the leachates. S.M. was responsible for the leaching experiment on the ²⁴⁴Cm-doped glass; M.T., S.M., C.J., and S.P. were involved in data interpretation. All authors helped to edit the paper.

## COMPETING INTERESTS

The authors declare no competing interests.

## ADDITIONAL INFORMATION

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41529-021-00183-4.

**Correspondence** and requests for materials should be addressed to M.T.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021