Investigations towards incorporation of Eu3+ and Cm3+ during ZrO2 crystallization in aqueous solution

Scientific Reports volume 13, Article number: 12276 (2023) Cite this article

135 Accesses

4 Altmetric

Metrics details

Nuclear energy provides a widely applied carbon-reduced energy source. Following operation, the spent nuclear fuel (SNF), containing a mixture of radiotoxic elements such as transuranics, needs to be safely disposed of. Safe storage of SNF in a deep geological repository (DGR) relies on multiple engineered and natural retention barriers to prevent environmental contamination. In this context, zirconia (ZrO2) formed on the SNF rod cladding, could be employed as an engineered barrier for immobilization of radionuclides via structural incorporation. This study investigates the incorporation of Eu3+ and Cm3+, representatives for trivalent transuranics, into zirconia by co-precipitation and crystallization in aqueous solution at 80 °C. Complementary structural and microstructural characterization has been carried out by powder X-ray diffraction (PXRD), spectrum imaging analysis based on energy-dispersive X-ray spectroscopy in scanning transmission electron microscopy mode (STEM-EDXS), and luminescence spectroscopy. The results reveal the association of the dopants with the zirconia particles and elucidate the presence of distinct bulk and superficially incorporated species. Hydrothermal aging for up to 460 days in alkaline media points to great stability of these incorporated species after initial crystallization, with no indication of phase segregation or release of Eu3+ and Cm3+ over time. These results suggest that zirconia would be a suitable technical retention barrier for mobilized trivalent actinides in a DGR.

Most of the high-level radioactive waste produced by fission of 235UO2 will eventually be stored in a deep geological repository (DGR). The waste consists primarily of the spent nuclear fuel (SNF), i.e., UO2, containing various fission products and transuranics, such as Pu, Np, Am, and Cm. Several of these elements contribute to the long-term radiotoxicity of the SNF, which ultimately will return to the initial level of the natural uranium ore after 105–106 years1. For the nuclear waste disposal safety assessment, it is important to understand the chemical behavior of the long-lived, radiotoxic transuranic elements after potential water intrusion into the repository, and the subsequent weathering or corrosion of barrier materials and the SNF matrix itself. Thus, immobilizing reactions, such as the incorporation of transuranics in secondary phases will play a role in determining the fate of these radionuclides in the geosphere2,3. One of the first materials that radioactive elements from the waste matrix can interact with is Zircaloy, which is the fuel rod cladding material used in pressurized water reactors. Zirconium-based materials are promising in the context of radiological contamination prevention due to their long-term stability as well as high dopant capacity before undergoing phase separation4,5,6. Already during reactor operation, a passivating zirconia (ZrO2) corrosion layer is formed on the cladding surface7. In contact with water in a DGR environment, slow dissolution and recrystallization of the ZrO2 layer may lead to incorporation of potentially mobile transuranics from the SNF matrix. Zirconia has three naturally occurring polymorphs. At ambient temperature and pressure conditions, the monoclinic (m) phase is the thermodynamically most stable structure. The tetragonal (t) and cubic (c) polymorphs can be stabilized at high temperatures or by incorporation of various aliovalent (e.g. Mg2+, Ca2+, Fe3+, Y3+, La3+, Nb5+)3,8,9,10,11,12,13,14 and isovalent (e.g. Ce4+, Ti4+)15,16 cations. Incorporation of subvalent cations is accompanied by oxygen vacancy formation in the host structure to retain charge neutrality. These vacancies will thereby reduce the coordination number of the cations in the structure, typically of the Zr4+ host. The incorporation of subvalent dopants in ZrO2 has been extensively studied using high-temperature (1000–1500 °C) synthesis methods, including co-precipitation or solid-state reaction routes yielding a crystalline, doped ZrO2 material. Much less is known about the uptake mechanism of subvalent cations, specifically actinides, into the ZrO2 crystal structure during crystallization in aqueous solution in the context of nuclear waste management.

The formation rate of a crystalline zirconia material under aqueous conditions has been shown to depend on the temperature (90–120 °C) and the pH (pH 3–13) of the solution. In very alkaline media a mixture of the tetragonal and monoclinic zirconia structures and a full crystallization of the initially amorphous zirconia precursor has been shown to occur after less than 100 h17. Zirconia demonstrates a higher tendency to transform from the tetragonal, often referred to as meta-stable, into the monoclinic phase upon hydrothermal aging, which is influenced mostly by the crystallite size18. Shukla and Seal19 summarized the effects of various parameters such as surface energy, interfacial energy, etc. on determining the stability of the undoped tetragonal zirconia particles, ultimately resulting in a critical crystallite size upon which the free energy favors the thermodynamically more stable monoclinic phase. Some of the investigated parameters included the surface energy, the interfacial energy, strain energies and multiple other effects that play into the stability range of this phase. The critical crystallite size for the tetragonal structure to be thermodynamically stable was estimated to be in a range between 10 and 33 nm.

Upon incorporation of dopants into the zirconia crystal structure, the high-temperature t + c crystal structures are stabilized during hydrothermal synthesis. The kinetics of the crystallization from the amorphous hydrous precursor, depending on different dopants and concentrations, show slightly faster crystallization rates compared to undoped zirconia20. This is likely due to the coexistence of stabilized tetragonal and/or cubic structures in the zirconia material. Low-temperature hydrothermal approaches after successful co-precipitation yield fully incorporated species with no sign of phase segregation21.

This study aims to understand the incorporation of trivalent actinides (Cm3+) and lanthanides (Eu3+) into ZrO2 in aqueous solution starting from a hydrous zirconia precursor. A pH of 12 and a temperature of 80 °C are chosen for the aqueous syntheses in the present study as both alkaline conditions and elevated temperatures have been shown to enhance the crystallization rate and conversion of hydrous zirconia to ZrO217,22. Moreover, high pH values have been shown to arise in the near-field of a DGR following e.g. canister corrosion and weathering of cementitious backfill and sealing materials23,24,25, while elevated temperatures are expected in the presence of the heat-generating SNF matrix, depending on its age and composition26,27. The experiments are designed to study the influence of the dopants, their concentration, and the hydrothermal aging time on the incorporation mechanism. The two studied dopants have been chosen due to their + 3 oxidation state being the most stable one in the aqueous phase and their excellent luminescence properties making it possible to study their coordination environment via luminescence spectroscopy. Powder X-ray diffraction (PXRD) and spectrum imaging analysis based on energy-dispersive X-ray spectroscopy in scanning transmission electron microscopy mode (STEM-EDXS) provide further insights into the phase composition and the spatial distribution of the elements in the ZrO2 structure.

The Eu3+- and Cm3+-doped hydrous zirconia precursors were analyzed by PXRD without further treatment following precipitation. Complete uptake of the dopants was confirmed by analysis of the supernatant solution by ICP-MS for Eu and LSC for Cm, with the concentrations being below the detection limit for each method. Thereby, the targeted dopant concentrations in the zirconia solid phases of 30 ppm (Cm3+) and 500 ppm or 1 mol% (Eu3+) can be confirmed. All precursors were X-ray amorphous. With increasing aqueous synthesis time at the chosen synthesis pH and temperature of 12 and 80 °C, respectively, a phase transformation of the amorphous phase to crystalline zirconia takes place. The diffractograms for selected Eu3+- and Cm3+-doped samples are compiled in Figure S1. For the Eu3+-containing samples, the crystalline fraction of the material could be quantified with known amounts of silicon as internal standard, while the ZrO2 phase composition was obtained via Rietveld refinement. Due to the very similar reflection positions of the tetragonal (t) and cubic (c) phases, and the relatively broad Bragg reflections, a distinction between these phases could not be made. The Rietveld refinement was thus performed considering only the monoclinic and tetragonal ZrO2 modifications. Results of the refinement have been compiled in Table S1 in the SI. The overall quantity of these higher symmetry polymorphs (t + c) is given in the following text. The quantity of crystalline ZrO2 in the Eu3+-doped samples and the phase composition of the Eu3+- and Cm3+-doped samples are shown in Fig. 1. The transformation of the amorphous precursor to the crystalline material increases notably after 10 and 16 days for samples containing 500 ppm (Fig. 1a) and 1 mol% Eu3+ (Fig. 1b), respectively. In both cases, the samples undergo almost full crystallization within a few days, reaching > 80% crystallinity after 21 days (500 ppm Eu) and 32 days (1 mol% Eu). For both dopant concentrations, the samples eventually fully crystallize after a hydrothermal synthesis time of > 1 year. Immediately after the initial crystallization, all Eu3+- and Cm3+-doped samples show a mixture of both the monoclinic and the t + c modifications of zirconia (Fig. 1a–c, red filled circles). With increasing synthesis time, the monoclinic crystal structure (Fig. 1a–c, black filled squares) becomes dominant, accounting for almost 80% of the phase composition, while the t + c phases constitute the remaining 20% of the overall composition for the doped samples.

Phase composition of the ZrO2 samples containing 500 ppm Eu3+ (a), 1 mol% Eu3+ (b), and 30 ppm Cm3+ (c), with the monoclinic phase depicted in squares (black), the tetragonal and cubic phase in circles (red), and the overall crystallinity of the Eu3+-doped samples shown as triangles (blue) determined with Si as internal standard. Synthesis conditions: T = 80 °C, 0.5 M NaCl, pH 12.

The crystallite size of all samples was determined using the Debye–Scherrer Eq. (1) and is given in Table 1.

A slight increase of the crystallite size can be seen with increasing hydrothermal treatment time for all dopants and concentrations. The different ZrO2 polymorphs and the chosen dopant concentrations, however, do not have a significant effect on the crystallite size. Moreover, the crystallite size seems to be independent on the dopant element considering the similar crystallite sizes for the Eu- and Cm-doped ZrO2.

The morphology and element distribution in the ZrO2 samples treated for 32 days and 399 days with 1 mol% Eu were investigated with STEM-based analyses. Looking at the HAADF-STEM images in Fig. 2A,C,E,G, ZrO2 seems to form agglomerates which are mainly composed of smaller, rather rounded particles having a size in the nanometer range. With increasing reaction time, the agglomerates tend to grow larger. EDXS-based element mapping (Fig. 2B,D,F,H) of the agglomerates shows a homogeneous distribution of Zr (blue) and O (red) (see also Figs. S2, S3). Regarding Eu (green), there seems to be a preferential accumulation at the periphery of the agglomerates, while a lower amount of europium (compared to the Zr and O signals) is observed in the bulk structure of the material. This is the case for both samples, taken during the crystallization process and after completion of the phase transformation, indicating that accumulation of the Eu dopant at the periphery of the particles takes place during the growth and crystallization and remains there after crystallization has been attained.

STEM-based analysis for the Eu:ZrO2 samples containing 1 mol% dopant after 32 days (A–D) and 399 days (E–H) of thermal treatment time. (A,C,E,G) are HAADF-STEM images, while (B, D, F, H) show the corresponding superimposed EDXS-based element distribution maps for Zr (blue), Eu (green), and O (red).

The successful co-precipitation of the dopant cations was confirmed in experiments accounting for surface adsorption or the formation of binary phase mixtures as alternative mechanisms for the removal of Eu or Cm from solution during precipitation. Luminescence spectroscopy enables insights into the structural uptake of a dopant of interest into a solid phase, independent of its crystallinity, and was successfully applied in previous studies regarding nuclear waste forms28,29. These experiments are described in detail in the supporting information (section “Luminescence spectroscopy”, Figs. S4, S5) and confirm that the dopants become incorporated into the hydrous zirconia phase during precipitation.

To get a first impression of the speciation or speciation changes of Eu3+ and Cm3+ within the zirconia structure at different aqueous synthesis times, luminescence emission spectra were collected at room temperature following indirect laser excitation in the UV-range (λex = 394 nm for Eu3+ and λex = 396.6 nm for Cm3+). The luminescence emission spectra of the ZrO2 samples doped with 30 ppm Cm3+ (Fig. 3a) and 1 mol% Eu (Fig. 3b) at selected aqueous synthesis times up to 399 days are shown in Fig. 3.

Luminescence emission spectra of the Cm:ZrO2 samples containing 30 ppm Cm3+ (a), and the Eu:ZrO2 samples containing 1 mol% Eu3+ (b). The Cm3+ spectra are normalized to the area under the curve, while the Eu3+ spectra are normalized to the overall intensity of the 7F1 band (magnetic dipole transition, 583–603 nm). Unique features for different crystal structures are indicated with arrows (m monoclinic, t tetragonal, c cubic).

The emission spectra of all synthesized Eu:ZrO2 (500 ppm and 1 mol%) samples can be found in the supporting information (Figs. S6, S7). In the beginning of the synthesis where the solid phases are dominated by amorphous material, all Eu-doped samples show similar emission spectra independent of their dopant concentration. The signals of the 7F0-2 \(\leftarrow\) 5D0 transitions are very broad, consistent with Eu incorporation in an amorphous solid phase. Following spectral decomposition of the recorded composite spectra, using the different luminescence lifetimes of the incorporated Eu species, pure component spectra for two non-equivalent Eu species could be obtained Figure S8. The luminescence lifetimes of 364 ± 22 μs and 1010 ± 125 μs were obtained from the decay curves (Figure S9). These lifetimes can be correlated with the number of H2O/OH− entities directly coordinating to Eu, using Eq. (2). The shorter lifetime corresponds to 2.3 water molecules in the Eu coordination sphere, while the longer lifetime indicates an almost complete loss of hydration water, with 0.4 remaining hydration water molecules in the first coordination sphere. With Eu incorporated in a hydrous solid structure, the rather short lifetimes can be better described by O–H vibrational quenching from the hydrous zirconia precursor rather than from coordinating water molecules from the aqueous solvent. As the hydrous zirconia structure can be described by the presence of several non-equivalent cation sites, which differ in their coordination to water of hydration in the amorphous structure, the presence of several non-equivalent Eu species incorporated at these different host cation sites can be understood30. Consistent with the Eu-doped samples, two non-equivalent species with different luminescence lifetimes were found for the Cm-doped samples (Figure S5). These lifetimes correspond to 5.1 and 1.6 hydration water entities coordinating to the Cm3+ cation, which speaks for an enhanced quenching of the water of hydration in the amorphous solid phase in comparison to Eu3+.

With increasing synthesis time and increasing crystallinity of the solid phase, the bands in the Eu3+ emission spectra become narrower, and the Stark splitting of the 7F1 and 7F2 bands becomes visible (Fig. 3b). The 7F2 band (605–640 nm) in the recorded emission spectra show a prominent peak at approximately \(\lambda\) = 613 nm. This peak is indicative of Eu incorporation into monoclinic ZrO231,32,33. With the onset of the formation of the crystalline zirconia structures, a small peak on the blue side of the 7F2 \(\leftarrow\) 5D0 transition, at around \(\lambda\) = 607 nm becomes visible. This peak has been shown in several studies to arise from Eu3+ incorporation in the tetragonal or cubic ZrO2 modifications31,33,34,35. Eu3+ in the monoclinic environment appears to reach its maximum for the 32 d sample. At longer aqueous synthesis times, a slight decrease of the 613 nm peak intensity can be observed, while the shoulder at 607 nm increases and the bands in the spectra become broader.

Similarly, the Cm spectra (Fig. 3a) show clear narrowing in comparison to the spectrum recorded for Cm3+ in the hydrous zirconia precipitate, with increasing aqueous synthesis time, pointing toward an increasingly ordered local environment around the Cm3+ cation with increasing crystallinity of the solid phase. The emission peak narrowing is accompanied by an extremely large red-shift of the spectra, to approximately 650 nm. To the best of our knowledge, such a large shift has never been observed before for Cm3+ incorporation in a solid matrix. The large shift implies a strong perturbation of the energy levels for Cm incorporated in the crystalline ZrO2 environment, which could be caused by an exceptionally short Cm–O bond, and/or a slight covalent character of the bond36. Understanding the underlying reason for this extreme perturbation, however, is not the scope of the present study and will be explored in detail in further studies. Although the incorporation of Cm3+ in the different structural polymorphs (monoclinic vs. tetragonal or cubic) cannot be deduced from the recorded emission spectra, the strong red-shift of the spectra excludes surface sorption as an uptake mechanism under the applied synthesis conditions37. Moreover, the emission peak width of the 32 days sample is narrower than the peak width after 117 days aqueous synthesis time, which is consistent with the Eu3+ data discussed above, and points toward the presence of several Cm3+ environments in these solid phases.

For a better understanding of the local environment of the trivalent cations in the ZrO2 crystal structure, site-selective laser excitation experiments were conducted at 10 K for selected crystalline samples. In these experiments, tunable laser excitation is used to selectively probe non-equivalent species in the solid matrix. For Eu3+, the electrons are directly and selectively promoted from the 7F0 ground state to the emitting 5D0 excited state. This technique was performed on the 500-ppm-doped samples after 32 d and 460 d and the 1 mol%-doped samples after 32 d and 399 d aqueous synthesis times. For Cm3+, laser excitation from the 8S′7/2 ground state to the 6D′7/2 excited state, typically denoted with the letter A, was used. The excited A-state of the curium ion undergoes crystal field splitting, giving rise to a maximum of four crystal field states A1−A4. Selective laser excitation to all these crystal field states was performed for samples synthesized for 32 days and 117 days.

The excitation spectra of the crystalline samples for both Eu concentrations and synthesis times are shown in Fig. 4.

Excitation spectra of Eu:ZrO2 with 500 ppm Eu3+ (solid lines and filled symbols) and 1 mol% Eu3+ (dashed lines, open symbols) after 32 d (blue) and 460 or 399 d (green) aqueous synthesis times.

All samples show a broad peak around ~ 579.5 nm and a narrow peak at ~ 582.9 nm with a shoulder at lower wavelengths ~ 582 nm. The excitation spectra after 32 days as well as 460 days are almost identical for the 500 ppm concentration series indicating that once ZrO2 has almost completed crystallization, the Eu environment in the solid structure does not change. The 1 mol% concentration sample series shows the same three peak positions, however, with slightly lower contribution of the more red-shifted signals at 582 nm and beyond. Emission spectra at the corresponding excitation peak maxima for both europium concentrations are shown in Fig. 5.

Emission spectra collected at different excitation wavelengths for Eu:ZrO2 samples containing 500 ppm Eu3+ (a, solid lines) and 1 mol% Eu3+ (b, dashed lines) after 460 and 399 days.

The emission spectrum at 579.5 nm is very broad considering that the host structure is fully crystalline, and the luminescence spectroscopic investigations are conducted in a He-cooled cryostat to reduce thermal broadening of the luminescence signal. When recording a spectrum after a delay time of 8000 μs, a clear emission band belonging to the F0–F2 transition, at 607 nm, indicative of Eu3+ incorporation in a tetragonal/cubic environment, can be seen (Figure S10). This implies, that the broad excitation spectrum is characterized by a short-lived and long-lived species in the tetragonal/cubic crystal structure, which after the delay shows the emission of the longer-lived species only. In both cases, the emission spectrum shows maximum Stark splitting of the emission bands, indicative for Eu incorporation on a lattice site with low symmetry38. For the tetragonal and cubic crystal structures, reduced splitting is expected due to symmetry induced degeneracy. This shows that the Eu3+ environment in the crystal structure does not correspond to the tetragonal or cubic bulk structure due to lattice distortions caused by the difference in ionic radii of Eu and Zr, which are 1.066 and 0.84 Å for a coordination number of eight, respectively39. In addition, one oxygen vacancy is generated in the ZrO2 structure for every two incorporated europium ions, to preserve charge neutrality in the crystal structure. The presence of these vacancies are likely to decrease the local order, and consequently the local symmetry of the solid matrix.

The two peaks at higher wavelengths (582.1 and 582.9 nm) are associated with Eu incorporated into the monoclinic crystal structure, based on the emission spectra recorded at these excitation wavelengths, showing the characteristic 7F2 peak at 613 nm. Since the monoclinic crystal structure only has a single cation site that can be occupied by a dopant ion, the two detected monoclinic environments must have some fundamental differences. Both excitation peaks at 582.1 and 582.9 nm differ only minimally in their emission spectra, mainly in the relative intensities of the 7F1 and 7F2 bands. Thereby, the two Eu3+ environments are very similar, and may differ only with respect to their coordination to the surrounding oxygen anions and created oxygen vacancies, as previously found for Eu3+ incorporated in cubic ZrO234. More specifically, a coordination number of seven in the pristine monoclinic ZrO2 environment would arise when the oxygen vacancy is not directly coordinated to the Eu3+ cation, while a reduction of the coordination number to six would occur when a vacancy is directly coordinated to the incorporated Eu3+ cation. The vacancy is known to be slightly larger than the oxygen anion, thereby increasing the unit cell volume and, subsequently increasing the M–O distance (M = Zr4+ or Eu3+)40. Thereby, we assign the Eu3+ environment at 582.1 nm to Eu3+ incorporation in monoclinic ZrO2 with direct coordination to an oxygen vacancy, while the species at 582.9 nm is Eu3+ in a pristine m-ZrO2 environment.

Based on the area under the curve in the excitation spectra, the dopant shows preferential incorporation into the tetragonal/cubic structure although the monoclinic crystal structure is dominant in the samples based on XRD results after the beginning of crystallization (Fig. 1). More specifically, approximately 70 ± 5% of the Eu3+ signal can be attributed to the c + t phase, and the remaining 30 ± 5% to Eu3+ in m-ZrO2. This observation cannot be made by using UV excitation (Fig. 3b). At room temperature (UV excited spectra), all present Eu3+ species are simultaneously excited with the laser light. Radiationless relaxation from a higher-lying excited state to the emitting excited state occurs, and energy transfer processes, for example from the cubic/tetragonal phase to the monoclinic phase can occur, which is likely the reason for the overrepresentation of the monoclinic environment using this method. Alternatively, this could be explained by a lower transition probability for Eu in the cubic/tetragonal structure and a higher, and therefore overrepresented, probability of the transition for the monoclinic phase.

The excitation spectra of the Cm:ZrO2 samples for 32 days and 117 days hydrothermal treatment times are shown in Fig. 6.

Excitation spectra of the 30-ppm Cm:ZrO2 samples after 32 days (blue) and 117 days (orange) aqueous synthesis times.

The excitation spectra are characterized by three peaks. Peaks A3 and A2 are crystal field states of the 6D′7/2 term with higher energy, usually referred to as hot-band transitions, while the peak labeled A1 is the crystal field level with the lowest energy, from which luminescence emission back to the ground state occurs. Contrary to the sample crystallinity, the excitation peak of the 32 days sample is slightly narrower than for the sample subjected to a longer aqueous synthesis time. The luminescence emission spectra collected via excitation directly to the emitting crystal field state (A1) show extreme line-narrowing effects (see Figure S11), which can be explained by a continuum of very similar environments arising from the lack of long-range order in the solid structure. Direct excitation to the hot-band transition (A2), avoiding the line narrowing effect, reveals multiple Cm3+ environments for the 32 days and 117 days samples, Fig. 7.

Emission spectra of ZrO2 samples doped with 30 ppm Cm3+ after 32 days (a) and 117 days (b) collected at different excitation wavelengths between 623 and 630 nm (excitation to A2 level).

The overall luminescence intensity of the 32 days sample is poor in comparison to the 117 days sample. However, in both cases two broad peaks with emission peak positions at approximately 639 nm and 642.5 nm can be identified. These emission peak positions are consistent with Cm3+ incorporation in a t + c ZrO2 environment41. The narrower peak at 649 nm (32 days sample) or 648 nm and 649 nm (117 days sample) are thereby assigned to Cm3+ in a monoclinic ZrO2 structure, in accordance with the Eu3+ luminescence data. The fact that the peak positions of Cm3+ incorporation in m-ZrO2 are almost identical, points toward very similar coordination spheres which in turn is consistent with Cm3+ coordination to seven oxygen anions or six anions and one vacancy, respectively.

The amorphous zirconia structure containing Eu starts crystallizing after approximately 10 days at a constant synthesis pH and temperature of 12 and 80 °C, respectively (Fig. 1). Other hydrothermal zirconia studies showed that at elevated temperatures of \(\ge\) 95 °C, this initial transformation occurs much faster, sometimes already after a few hours20,22,42. This indicates a temperature dependence of the crystallization rate, which to some extent can be attributed to the slightly lower solubility of the hydrous zirconia precursor at lower temperatures43. Additionally, the dissolution/precipitation reaction is strongly pH-dependent. The reported solubility of the hydrous zirconia precursor is several orders of magnitude lower in the near-neutral pH range than for very acidic or alkaline systems43,44. According to the studies of Clearfield and Denkewicz et al., the formation of the different zirconia modifications depends on the starting pH43,45. The thermodynamically stable monoclinic zirconia modification is formed over the entire pH range, while the tetragonal and cubic phases are only observed in alkaline systems. Over time, the tetragonal structure undergoes a transformation into the thermodynamically stable monoclinic modification43. Other studies show that the dopant concentration influences the kinetics of the tetragonal to monoclinic transformation during hydrothermal aging46. Lower dopant concentrations tend to increase the transformation rate from the meta-stable tetragonal to the monoclinic phase. The dopant size also plays a crucial role in the kinetics of the transformation. Dopants with a larger size mismatch to the zirconium host were shown to have a much slower conversion rate than dopants having a similar ionic radius. In all cases, the tetragonal to monoclinic transformation reached a maximum conversion at which no further material underwent phase transformation. This value was reported to be between 80 and 85% for the material that underwent phase transformation from tetragonal to monoclinic. Higher dopant concentrations reduce that value due to lattice stabilization caused by the dopant46,47.

In the current study, no phase transformation of a meta-stable tetragonal phase was observed over the course of > 1 year. In contrast, the amount of the t + c phases increases slightly but constantly over time for the samples with 500 ppm Eu3+ or 30 ppm Cm3+ (Fig. 1). This is likely related to the small size of the crystallites in all samples (Table 1). It has been shown that the tetragonal to monoclinic phase transformation is suppressed below a certain particle size threshold, thereby stabilizing the tetragonal phase. An in-depth study by Garvie18 investigated the stability of tetragonal crystallites formed after precipitation from aqueous solution. Shukla and Seal19 summarized Garvie’s study amongst others to show the impact of various effects (i.e. surface energy, interfacial energy, strain energy, anions and more) on the phase stability. They concluded a critical crystallite size of 10–33 nm before the tetragonal phase would become meta-stable and transform into the monoclinic phase14,18,19. The maximum tetragonal crystallite size in this study is 18 nm (Eu:ZrO2, 1 mol%, 112 days), i.e. well within the reported range resulting in a stable tetragonal or t + c phase in the samples. Prolonged synthesis times of up to 460 days do not affect the crystallite size, which would explain the absence of a tetragonal to monoclinic phase transformation. The stability of the tetragonal phase can be attributed to the role of interfacial energy caused by the agglomeration and modification of the surface energy of the zirconia particles48, with the tetragonal phase remaining stable under these hydrothermal conditions. Besides mere size effects, subvalent dopants are known to stabilize the tetragonal and cubic modifications. This stabilization mechanism is attributed to the reduction of the coordination number, especially for the small Zr4+ host cation, via the introduction of oxygen vacancies following trivalent cation incorporation, to preserve charge neutrality in the crystal structure. Especially in the sample with 1 mol% Eu3+, the occurrence of the t + c modifications in the ZrO2 samples can be reasonably assumed to arise from both crystallite size effects and the above-discussed stabilization by the subvalent dopant49.

Existing studies have reported phase segregation of dopants from stabilized zirconia structures during hydrothermal conditions. For example, Zhang et al.46 reported immediate phase segregation of oversized (Nd, La) or undersized dopants (Al) from the ZrO2 matrix for dopant concentrations of 0.6–0.8 mol%48. According to the combined microscopic and spectroscopic investigations, these findings were not confirmed in the present study, using oversized trivalent dopants in concentrations ranging from 30 ppm to 1 mol%. Phase segregation would lead to either an accumulation of the trivalent cations on the ZrO2 surface as adsorption complexes, or the precipitation of the dopant hydroxides or hydroxycarbonates from oversaturated solutions. Both adsorption and precipitation processes could be excluded based on the luminescence spectroscopic investigations described in the supporting information (section “Luminescence spectroscopy”, Figs. S4, S5). For the highest Eu3+ dopant concentration of 1 mol%, our STEM-EDXS data shows accumulation of the trivalent cation close to or at the surface of the ZrO2 particle agglomerates (Fig. 2), which we will refer to as superficial dopant species in the following text. A relatively lower Eu3+ concentration can be detected within the agglomorates, i.e. in the bulk (Figs. S2, S3). The amount of superficial vs. bulk species does not increase over aqueous synthesis time, but the distribution of the dopant within the zirconia matrix seems to remain the same once crystallization of ZrO2 has started. By examining the luminescence data for the Eu3+- and Cm3+-doped ZrO2 samples with increasing aqueous synthesis time, a complex but unchanged speciation of the f-element cations in the solid structure can be deduced. From the Eu3+ luminescence data, a clear preference for incorporation within the t + c modifications over the monoclinic one can be seen, based on the area under the measured excitation spectra shown in Fig. 4, Eu3+ incorporation in the monoclinic phase amounts to approximately 30 ± 5% of the Eu3+ speciation in the crystalline ZrO2 solid material, as deduced from the combined integrals of the peaks at 582.1 and 582.9 nm relative to the overall integral of the measured spectra. The amount of Eu3+ associated with the higher symmetry polymorphs, characterized by the broad signal at approximately 579.5 nm in Fig. 4, is 70 ± 5%. This amount is further divided between two Eu3+ environments which differ in their luminescence decay times. The presence of a short-lived Eu3+ species, which is not adsorbed on the ZrO2 surface, could arise from e.g. a Eu3+ ion partially coordinating to water in the surrounding solution or from Eu3+ coordination to surface-terminated hydroxyl groups at the very superficial sites in the ZrO2 particles. The Eu3+ species with a significantly longer luminescence lifetime is indicative of complete isolation of Eu3+ from the surrounding aqueous solution, which can be attributed to Eu3+ incorporation into the c + t-ZrO2 bulk material. Based on the combined microscopic and spectroscopic results, we assign the short-lived Eu3+ species to the visible Eu3+ accumulation at the periphery of the agglomerates. Whether the accumulation can be attributed to an overall smaller size of the particles, thereby increasing the surface-to-bulk ratio, can unfortunately not be deduced from our TEM-based analysis due to the agglomerate formation preventing the unambiguous identification of individual ZrO2 particles.

Finally, Eu3+ and Cm3+ incorporation within the monoclinic crystal structure is characterized by the presence of at least two non-equivalent species with clearly narrower luminescence spectra than observed for incorporation in t + c ZrO2. This is rather surprising, as the monoclinic crystal structure only has one crystallographically unique Zr4+ cation site and, consequently, only one Eu3+ and Cm3+ environment would be expected. The Eu3+ emission spectra collected via selective excitation at the peak maxima of the two non-equivalent environments in the monoclinic modification, however, are very similar (Fig. 5b), pointing toward very subtle differences between these species. One potential explanation for the presence of two non-equivalent species with very similar emission spectra is the presence of oxygen vacancies in nearest-neighbor or next-nearest-neighbor positions with respect to the incorporated dopant. The position of the generated oxygen vacancies with respect to Eu3+ is rather well-described for Eu3+ incorporation in the cubic ZrO2 polymorph, but has hitherto not been described for the monoclinic structure34,50,51. Although a definite assignment of the different trivalent environments in the monoclinic ZrO2 structure would require additional studies, the formed species were observed to be stable and do not change much during extended periods of aqueous synthesis. We can thereby conclude, that the incorporation of trivalent f-elements in the zirconia polymorphs formed after crystallization show a high stability without any signs of dopant segregation from the crystal structure.

This study investigated the role of ZrO2, a corrosion product forming on the nuclear fuel rod cladding surface, to act as engineered barrier material via structural incorporation of potentially mobilized radionuclides in a DGR. Experimental conditions such as elevated temperature and solution chemistry were chosen to mimic conditions that may arise in a generic DGR as a result of radioactive decay and the corrosion of the canister and cementitious backfill and sealing materials. Our experiments revealed a successful co-precipitation of Eu3+ as well as Cm3+ with zirconia, forming an amorphous hydrous precursor containing multiple incorporation sites of the dopant elements at the initial stage. Upon progressing synthesis duration, the material underwent crystallization yielding a mixture of mainly the thermodynamically stable monoclinic structure and a smaller fraction of the tetragonal and cubic polymorphs. The occurrence of the tetragonal and cubic phases arises from the structural stabilization due to the introduction of oxygen vacancies reducing the coordination number of the Zr lattice as well as the small crystallite size (< 20 nm). The various phases were stable and did not undergo a phase transformation, most likely due to agglomeration of crystallites, thus modifying the surface energy. Immediately after crystallization, the trivalent dopants could be located in the different crystal structure polymorphs. Although the phase distribution showed that most of the zirconia crystallized in the monoclinic phase, for Eu, a preferential incorporation into the tetragonal and cubic structure was observed. In the case of the Cm incorporation, this preference for the high-temperature polymorphs was not evident. Over the entire synthesis duration of up to 460 days, the dopant remained structurally incorporated in the crystal structure and did not undergo phase segregation. TEM-based analysis revealed a clear accumulation of the dopant at the particle peripheries resulting in an additional dopant environment. This environment in characterized by a clearly shorter luminescence lifetime in comparison to incorporated species, either as a result of direct Eu3+ coordination to water in the surrounding solution or from Eu3+ coordination to surface-terminated hydroxyl groups at the very surface of the ZrO2 particles. Independent of the crystal structure, the dopant resides in a superficial and bulk environment. Complementary XRD, STEM and luminescence spectroscopy measurements revealed trivalent cations to be structurally embedded into ZrO2 during slightly elevated temperature hydrothermal experiments. These results indicate that ZrO2, present as cladding corrosion material of spent fuel rods in a DGR, would additionally be suitable to act as engineered barrier material for immobilization of mobilized trivalent cations by structural uptake of these mobilized species, in case of water intrusion into a DGR. Future studies will focus on the role of pH towards the immobilization of trivalent actinides within the different zirconia polymorphs. In addition, spectroscopic investigations of Eu3+- and Cm3+-doped monoclinic ZrO2 obtained via high-temperature calcination routes are planned to elucidate the origin of multiple dopant environments in the m-ZrO2 crystal structure and the unusually large spectral shifts acquired for the Cm3+-doped material.

M3+-doped hydrous zirconia samples were prepared by dissolving ZrOCl2·8 H2O (Sigma-Aldrich, > 99.99%) in 0.01 M HCl using a 1:2 ratio mass (g)/volume (mL), resulting in a Zr4+ concentration of 1.55 M. The trivalent cation was added to the Zr solution from a stock solution of 0.1 M EuCl3·6 H2O in 1 mM HCl or 0.1 mM 248Cm3+ in 0.01 M HClO4 to obtain the desired dopant concentrations of 500 ppm, 1 mol% (Eu), and 30 ppm (Cm) in the hydrous zirconium precursor. The concentration of the dopants refers to the number of M3+ per total M3+ + Zr4+ cations in doped ZrO2, assuming quantitative incorporation of the cations from solution into the formed co-precipitate. Subsequently, this mixed solution was added dropwise under constant stirring to a 0.5 M NaCl solution in a ratio of approximately 1:20 for Eu3+ doped and 1:28 for Cm3+ doped samples (M3+/M4+ solution: NaCl). The pH was readjusted to its initial value of 12 with 2 M NaOH. The suspension was left for 24 h at room temperature before the solid phase was recovered via centrifugation of the suspension for 30 min at 4000 rpm. An aliquot of the supernatant was taken for Inductively-Coupled Plasma Mass Spectrometry (ICP-MS for Eu) or Liquid Scintillation Counting (LSC for Cm) analysis to ensure complete uptake of the dopant. In all cases, the dopant concentration in the supernatant was below the detection limit, thereby confirming the quantitative uptake of the dopants by the formed solid phase. The rest of the liquid phase was discarded, and the wet solid was washed twice with MilliQ water (18.2 MOhm). The washing step included addition of 40 mL of MilliQ water into the Greiner tubes, full resuspension of the solid phase, and immediate centrifugation of the suspension, subsequently discarding the washing solution. After the second washing step, the wet solid was dried for 24 h at 80 °C. The dry solid was ground in an agate mortar into a fine powder.

To probe whether the M3+ dopant was (1) adsorbed onto the Zr-containing precipitate or (2) formed a separate phase, such as a hydroxide or hydroxy carbonate rather than a co-precipitate, additional samples were prepared. Adsorption of Eu3+ on the amorphous precipitate was investigated at pH 12. For this purpose, non-doped hydrous zirconia was precipitated in the absence of Eu3+, as described above. 60 mg of the precipitate were thereafter suspended in 40 mL of 0.5 M NaCl at pH 5. 25 µL of a 0.01 M EuCl3 \(\cdot\) 6 H2O stock solution were added. The pH 12 was adjusted gradually over the period of two days by manual addition of small aliquots of NaOH under continuous stirring to avoid large pH-fluctuations and, thereby, to minimize the risk of immediate precipitation of the europium. The samples were then equilibrated at room temperature for 24 h. After the experiment was finished, the suspension was centrifuged for 30 min and the wet solid was dried at 80 °C for 24 h. The pure Eu precipitate was synthesized by dropwise addition of the acidic 0.1 M Eu3+ stock solution to a 0.5 M NaCl solution at pH 12. The suspension was centrifuged for 30 min and afterwards washed twice with MilliQ water. The white precipitate was dried at 80 °C for 24 h.

For the Eu-doped experiments, 120 mg of the doped zirconia powder was submerged in 40 mL 0.5 M NaCl and the pH was adjusted to 12 with 2 M NaOH. This resulted in a 1:3 ratio volume (mL)/mass (mg). The suspension was transferred into a 50 mL plastic Greiner tube and stored in an oven at 80 °C for the desired amount of time. After the thermal treatment, the sample was taken out of the oven and cooled down to room temperature. The sample was then centrifuged for 30 min and washed with 20 mL MilliQ water twice, following the procedure discussed above. The solid was dried at 80 °C for 24 h and then ground in an agate mortar into a fine powder for further analysis.

For the Cm-containing samples, 100 mg of the powder was submerged in 13.3 mL 0.5 M NaCl with pH values of 5 and 12, resulting in a 1:7.7 ratio volume (mL)/mass (mg). The hydrous zirconia sample at pH 5 was solely used for comparison to adsorption investigations, to elucidate the role of surface adsorption in the removal of f-elements from solution in the zirconia precipitation step. Time-dependent hydrothermal studies were only conducted for the samples prepared at pH 12. Afterwards, the samples underwent the same treatment as the Eu-doped samples. All Cm-involving tasks were performed insight a glovebox with negative pressure in N2 atmosphere and carbonate-free solutions.

To follow the crystallization of the ZrO2 solid phase as a function of the hydrothermal synthesis time and to identify the structural modifications (monoclinic, tetragonal, cubic) of the crystallized ZrO2, Powder X-ray diffraction (PXRD) analyses were performed on a Rigaku MiniFlex 600 using a Cu Kα X-ray source (40 kV/15 mA operation for X-ray generation). Samples were finely ground in an agate mortar. The samples containing inactive Eu3+ were measured on a silicon wafer in the 2θ range between 5° and 90° with a step size of 0.02° 2Θ/min while spinning the Si wafer. The radioactive Cm-containing samples were prepared inside a glovebox under N2 atmosphere using an air-tight sample holder consisting of a Kapton-domed silicon wafer. The Cm-containing samples were measured in a range between 5° and 60° 2Θ with a step size of 0.02° 2\(\Theta\)/min without spinning to avoid possible contamination of the sample holder. The crystallinity was only determined for the inactive samples by the addition of 20 wt% monocrystalline silicon as an internal standard to each sample, followed by thorough grinding of the sample mixture into a homogenous fine powder. For the Rietveld refinement, references of the monoclinic (space group:\(P{2}_{1}/c\)) and tetragonal (\(P{4}_{2}/nmc\)) zirconia modification as well as silicon (\(Fd\overline{3}m\)) were used in the PDXL2 software (Rigaku Corporation). More specifically, the database cards 04-004-4339 (m-ZrO2), 04-005-4479 (t-ZrO2), and 04-014-8844 (Si) were taken for the refinement. The cubic zirconia modification (\(Fm\overline{3}m\)) was not used for the refinement due to the peak width overlapping with the tetragonal modification. However, as we cannot exclude its presence in the synthesized solid phases, the quantity of the refined tetragonal phase is referred to as c + t-ZrO2 throughout the text. The background, peak shape, and preferred orientation effects were corrected or accounted for using B-spline, Split pseudo-Voigt, and March–Dollase models in the refinement.

The crystallite size was determined using the Debye–Scherrer Eq. (1):

Here, L represents the crystallite size perpendicular to the lattice plane of the Bragg reflection in nm. K is a shape factor, which is dependent of the shape of the crystals and the absolute value can only be determined empirically. In the current work, a shape factor of \(K=0.9\) was used, in agreement with previous studies22,42. The diffractograms were measured using a Cu X-ray source (\({\lambda }_{C{u}_{\alpha }}=0.154nm\)). \(\Delta \left(2\Theta \right)\) is the full width at half maximum (FWHM) of the diffraction maximum with Θ being the Bragg angle.

To examine the morphology of the ZrO2 solid phases following aqueous synthesis and the distribution of the dopants within the solid ZrO2 structure, high-angle annular dark-field scanning transmission electron microscopy (HAADF-STEM) imaging combined with spectrum imaging analysis based on energy-dispersive X-ray spectroscopy (EDXS) were performed at an accelerating voltage of 200 kV with a Talos F200X microscope equipped with an X-FEG electron source and a Super-X EDX detector system (FEI). Prior to STEM analysis, the specimen mounted in a high-visibility low-background holder was placed for 2 s into a Model 1020 Plasma Cleaner (Fischione) to remove possible contaminations. The TEM specimens were prepared by creating a very diluted suspension of the sample in MilliQ water and pipetting it onto a carbon-coated TEM grid and letting the water evaporate under ambient atmosphere.

The Eu3+ and Cm3+ environments in the solid phases were investigated via luminescence spectroscopy. Two different spectroscopic approaches were used as explained below.

All synthesized samples were subjected to luminescence spectroscopic investigations at room temperature, using a non-selective excitation wavelength of 394 nm or 396.6 nm, capable of exciting all non-equivalent Eu3+ species or Cm3+ species, respectively. These non-selective luminescence spectroscopic experiments were performed using a pulsed dye-laser (NarrowScan, Radiant Dyes) with a 1:1 dye mixture of Exalite 389 and Exalite 398, coupled to a Nd:YAG (Continuum, Surelite) pump laser. The Eu3+ luminescence emission spectra were recorded in a range between 575 and 645 nm, 1 µs after the exciting laser pulse, while a spectral range between 560 and 700 nm was used for Cm3+, 1 µs after the exciting laser pulse. The spectra were recorded using the Andor SOLIS software. The laser pulse energy was measured to be between 2 and 3 mJ. Luminescence emission spectra were detected using an optical multichannel analyzer (Shamrock 303i) with 300 or 600 lines/mm grating and an iCCD camera (iStar, Andor). Delayed spectra for the identification of the Eu3+ or Cm3+ species with different lifetimes were recorded by delaying the signal detection after the laser pulse up to 15 ms (Eu) or 800 μs (Cm). Luminescence lifetime measurements were performed by collecting multiple spectra and stepwise increasing the time delay of the signal detection and the laser pulse between 20 and 500 µs.

To determine the number of water molecules in the first ligand sphere the equations of Horrocks et al.52 (2) and Kimura et al.53 (3) can be used:

Based on the results of the non-selective luminescence spectroscopic investigations, the chosen samples were subjected to additional spectroscopic experiments with selective, tunable laser excitations in a He-cooled cryostat < 10 K (Janis and Sumitomo, SHI cryogenics group). The samples were placed in a copper sample holder covered by quartz glass. Therefore, a Nd:YAG pump laser (Continuum SureLite II) was coupled to a dye system (NarrowScan, Radiant Dyes) containing Rhodamine 6G/Rhodamine B (560–585 nm), DCM (600–660 nm), and Sulforhodamine (585–605 nm) dye. An optical multichannel analyzer (Shamrock 303i) with gratings of 150, 600, and 1200 lines/mm and an iCCD camera (iStar, Andor) were used. The laser energy was measured by the help of an energy meter (Newport 1918-R) which gave the mean energy while collecting one spectrum which was later used for normalizing the resulting excitation spectrum. The exact wavelength was measured with a wavelength meter (High Finesse WS-5). All spectra were recorded via the Andor SOLIS software.

The datasets used and/or analyzed during the current study available from the corresponding authors on reasonable request.

Gompper, K. Zur Abtrennung Langlebiger Radionuklide 152–167 (Radioaktivität und Kernenergie, 2001).

Google Scholar

Raison, P. E. & Haire, R. G. Zirconia-based materials for transmutation of americium and curium: Cubic stabilized zirconia and zirconium oxide pyrochlores. Prog. Nucl. Energy 38, 251–254 (2001).

Article CAS Google Scholar

Svitlyk, V., Weiss, S. & Hennig, C. Immobilization of radiotoxic elements with Y-stabilized zirconia: The thorium case. J. Am. Ceram. Soc. 105, 5975–5983 (2022).

Article CAS Google Scholar

Voice, P. J., Kowalewski, M. & Eriksson, K. A. Quantifying the timing and rate of crustal evolution: Global compilation of radiometrically dated detrital zircon grains. J. Geol. 119, 109–126 (2011).

Article ADS CAS Google Scholar

Stubican, V. S., Hink, R. C. & Ray, S. P. Phase equilibria and ordering in the system ZrO2-Y2O3. J. Am. Ceram. Soc. 61, 17–21 (1978).

Article CAS Google Scholar

Veevers, J. J. Beach sand of SE Australia traced by zircon ages through ordovician turbidites and S-type granites of the Lachlan Orogen to Africa/Antarctica: A review. Aust. J. Earth Sci. 62, 385–408 (2015).

Article ADS CAS Google Scholar

Cappia, F. et al. Detailed characterization of a PWR fuel rod at high burnup in support of LOCA testing. J. Nucl. Mater. 569, 153881 (2022).

Article CAS Google Scholar

Kelly, J. R. & Denry, I. Stabilized zirconia as a structural ceramic: An overview. Dent. Mater. 24, 289–298 (2008).

Article CAS PubMed Google Scholar

Li, K. et al. One-step preparation of CaO-doped partially stabilized zirconia from fused zirconia. Ceram. Int. 46, 6484–6490 (2020).

Article CAS Google Scholar

Heuer, C., Storti, E., Graule, T. & Aneziris, C. G. Electrospinning of Y2O3- and MgO-stabilized zirconia nanofibers and characterization of the evolving phase composition and morphology during thermal treatment. Ceram. Int. 46, 12001–12008 (2020).

Article CAS Google Scholar

Guo, F. & Xiao, P. Effect of Fe2O3 doping on sintering of yttria-stabilized zirconia. J. Eur. Ceram. Soc. 32, 4157–4164 (2012).

Article CAS Google Scholar

Ray, J. C., Pati, R. K. & Pramanik, P. Chemical synthesis and structural characterization of nanocrystalline powders of pure zirconia and yttria stabilized zirconia (YSZ). J. Eur. Ceram. Soc. 20, 1289–1295 (2000).

Article CAS Google Scholar

Zeng, Z., Liu, Y., Qian, F., Guo, J. & Wu, M. Role of a low level of La2O3 dopant on the tetragonal to monoclinic phase transformation of ceria-yttria co-stabilized zirconia. J. Eur. Ceram. Soc. 39, 4338–4346 (2019).

Article CAS Google Scholar

Kountouros, P. & Petzow, G. Defect chemistry, phase stability and properties of zirconia polycrystals. In Science and Technology of Zirconia V (eds Badwal, S. P. S. et al.) 30–48 (Technomic Publishing Company, Inc., 1993).

Google Scholar

Svitlyk, V., Weiss, S. & Hennig, C. Stability of doped zirconia under extreme conditions: Toward long-term and secure storage of radioactive waste. J. Am. Ceram. Soc. 105, 7831–7839 (2022).

Article CAS Google Scholar

Lv, H. et al. Preparation and properties of black Ti-doped zirconia ceramics. J. Mater. Res. Technol. 9, 6201–6208 (2020).

Article CAS Google Scholar

Štefanić, G., Musić, S. & Molčanov, K. The crystallization process of HfO2 and ZrO2 under hydrothermal conditions. J. Alloys Compd. 387, 300–307 (2005).

Article Google Scholar

Garvie, R. C. The occurrence of metastable tetragonal zirconia as a crystallite size effect. J. Phys. Chem. 69, 1238–1243 (1965).

Article CAS Google Scholar

Shukla, S. & Seal, S. Mechanisms of room temperature metastable tetragonal phase stabilisation in zirconia. Int. Mater. Rev. 50, 45–64 (2005).

Article CAS Google Scholar

Tsukada, T., Venigalla, S., Morrone, A. A. & Adair, J. H. Low-temperature hydrothermal synthesis of yttrium-doped zirconia powders. J. Am. Ceram. Soc. 82, 1169–1174 (1999).

Article CAS Google Scholar

Babić, V., Ivanda, M. & Štefanić, G. Phase development in the metastable solid solutions of ZrO2-YO1.5 system. J. Mol. Struct. 1185, 310–322 (2019).

Article ADS Google Scholar

Štefanić, G., Popović, S. & Musić, S. Influence of pH on the hydrothermal crystallization kinetics and crystal structure of ZrO2. Thermochim. Acta 303, 31–39 (1997).

Article Google Scholar

Samper, J., Naves, A., Montenegro, L. & Mon, A. Reactive transport modelling of the long-term interactions of corrosion products and compacted bentonite in a HLW repository in granite: Uncertainties and relevance for performance assessment. Appl. Geochem. 7, 42–51 (2006).

Google Scholar

Bildstein, O., Trotignon, L., Perronnet, M. & Jullien, M. Modelling iron–clay interactions in deep geological disposal conditions. Phys. Chem. Earth 31, 618–625 (2006).

Article ADS Google Scholar

Tyupina, E. K., Kozlov, P. P. & Krupskaya, V. V. Application of cement-based materials as a component of an engineered barrier system at geological disposal facilities for radioactive waste—a review. Energies 16, 605 (2023).

Article CAS Google Scholar

King, F. Nuclear waste canister materials: Corrosion behavior and long-term performance in geological repository systems. In Geological Repository Systems for Safe Disposal of Spent Nuclear Fuels and Radioactive Waste 2nd edn (eds Apted, M. J. & Ahn, J.) 365–408 (Woodhead Publishing, 2017).

Chapter Google Scholar

Delage, P., Cui, Y. J. & Tang, A. M. Clays in radioactive waste disposal. J. Rock Mech. Geotech. Eng. 2, 111–123 (2010).

Article Google Scholar

Holliday, K. et al. TRLFS of Eu3+ and Cm3+ doped La2Zr2O7: A comparison of defect fluorite to pyrochlore structures. J. Nucl. Mater. 433, 479–485 (2013).

Article ADS CAS Google Scholar

Huittinen, N. et al. Probing structural homogeneity of La1-xGdxPO4 monazite-type solid solutions by combined spectroscopic and computational studies. J. Nucl. Mater. 486, 148–157 (2017).

Article ADS CAS Google Scholar

Huang, C., Tang, Z. & Zhang, Z. Differences between zirconium hydroxide (Zr(OH)4·nH2O) and hydrous zirconia (ZrO2·nH2O). J. Am. Ceram. Soc. 84, 1637–1638 (2001).

Article CAS Google Scholar

Smits, K. et al. Europium doped zirconia luminescence. Opt. Mater. 32, 827–831 (2010).

Article ADS CAS Google Scholar

Hemmer, E. et al. Influence of the host phase on the vibrational spectra of europium-doped zirconia prepared by hydrothermal processing. J. Am. Ceram. Soc. 93, 3873–3879 (2010).

Article CAS Google Scholar

Marin, R. et al. Monitoring the t → m martensitic phase transformation by photoluminescence emission in Eu3+-doped zirconia powders. J. Am. Ceram. Soc. 96, 2628–2635 (2013).

Article CAS Google Scholar

Eibl, M. et al. Understanding the local structure of Eu3+- and Y3+-stabilized zirconia: Insights from luminescence and X-ray absorption spectroscopic investigations. J. Mater. Sci. 55, 10095–10120 (2020).

Article ADS CAS Google Scholar

Soares, M. R. N. et al. Structural and optical properties of europium doped zirconia single crystals fibers grown by laser floating zone. J. Appl. Phys. 109, 013516 (2011).

Article ADS Google Scholar

Edelstein, N. M., Klenze, R., Fanghänel, T. & Hubert, S. Optical properties of Cm(III) in crystals and solutions and their application to Cm(III) speciation. Coord. Chem. Rev. 250, 948–973 (2006).

Article CAS Google Scholar

Eibl, M. et al. A spectroscopic study of trivalent cation (Cm3+ and Eu3+) sorption on monoclinic zirconia (ZrO2). Appl. Surf. Sci. 487, 1316–1328 (2019).

Article ADS CAS Google Scholar

Binnemans, K. Interpretation of europium(III) spectra. Coord. Chem. Rev. 295, 1–45 (2015).

Article CAS Google Scholar

Shannon, R. D. Revised effective ionic radii and systematic studies of interatomic distances in halides and chalcogenides. Acta Cryst. 32, 751 (1976).

Article Google Scholar

Ohmichi, T., Fwkushima, S., Maeda, A. & Watanabe, H. On the relation between lattice parameter and O/M ratio for uranium dioxide-trivalent rare earth oxide solid solution. J. Nucl. Mater. 102, 40–46 (1981).

Article ADS CAS Google Scholar

Eibl, M. Actinide interaction with Zr-bearing phases: Spectroscopic investigations of An3+ sorption and incorporation reactions with zirconia. Ph.D. Thesis, Technische Universität Dresden (2020)

Arantes, T. M. et al. Stable colloidal suspensions of nanostructured zirconium oxide synthesized by hydrothermal process. J. Nanopart. Res. 12, 3105–3110 (2010).

Article ADS CAS Google Scholar

Denkewicz, R. P., TenHuisen, K. S. & Adair, J. H. Hydrothermal crystallization kinetics of m-ZrO2 and t-ZrO2. J. Mater. Res. 5, 2698–2705 (1990).

Article ADS CAS Google Scholar

Kobayashi, T., Sasaki, T., Takagi, I. & Moriyama, H. Solubility of zirconium(IV) hydrous oxides. J. Nucl. Sci. Technol. 44, 90–94 (2007).

Article CAS Google Scholar

Clearfield, A. Crystalline hydrous zirconia. Inorg. Chem. 3, 146–148 (1964).

Article CAS Google Scholar

Zhang, F. et al. Effect of cation dopant radius on the hydrothermal stability of tetragonal zirconia: Grain boundary segregation and oxygen vacancy annihilation. Acta Mater. 106, 48–58 (2016).

Article ADS CAS Google Scholar

Zhang, F. et al. 3Y-TZP ceramics with improved hydrothermal degradation resistance and fracture toughness. J. Eur. Ceram. Soc. 34, 2453–2463 (2014).

Article CAS Google Scholar

Garvie, R. C. Stabilization of the tetragonal structure in zirconia microcrystals. J. Phys. Chem. 82, 218–224 (1978).

Article CAS Google Scholar

Fabris, S., Paxton, A. T. & Finnis, M. W. A stabilization mechanism of zirconia based on oxygen vacancies only. Acta Mater. 50, 5171–5178 (2002).

Article ADS CAS Google Scholar

Borik, M. A. et al. Spectroscopy of optical centers of Eu3+ ions in partially stabilized and stabilized zirconium crystals. Opt. Spectrosc. 122, 580–587 (2017).

Article ADS CAS Google Scholar

Yugami, H., Koike, A., Ishigame, M. & Suemoto, T. Relationship between local structures and ionic conductivity in ZrO2-Y2O3 studied by site-selective spectroscopy. Phys. Rev. B 44, 9214–9222 (1991).

Article ADS CAS Google Scholar

De William, W. H. & Sudnick, D. R. Lanthanide ion probes of structure in biology. Laser-induced luminescence decay constants provide a direct measure of the number of metal-coordinated water molecules. J. Am. Chem. Soc. 101, 334–340 (1979).

Article Google Scholar

Kimura, T., Choppin, G. R., Kato, Y. & Yoshida, Z. Determination of the hydration number of Cm(III) in various aqueous solutions. Radiochim. Acta 72, 61–64 (1996).

Article CAS Google Scholar

Download references

Special thanks to Verona Kolpe for creating the graphical abstract. Furthermore, the use of the HZDR Ion Beam Center TEM facilities is acknowledged.

This work was supported by the Department of Energy, Office of Nuclear Energy under NEUP Grant DE-NE0009169 and the Bundesministerium für Bildung und Forschung (BMBF) AcE Grant (02NUK060). BMBF Grant No. 03SF0451 in the framework of HEMCP funded the TEM Talos. Open Access funding enabled and organized by Projekt DEAL.

Institute of Resource Ecology, Helmholtz-Zentrum Dresden-Rossendorf, 01328, Dresden, Germany

Lucas Opitz, Salim Shams Aldin Azzam, Sara E. Gilson & Nina Huittinen

Department of Chemistry, University of California, Irvine, Irvine, CA, 92697, USA

Lucas Opitz & Sarah C. Finkeldei

Institute of Ion Beam Physics and Materials Research, Helmholtz-Zentrum Dresden-Rossendorf, 01328, Dresden, Germany

René Hübner

Department of Materials Science and Engineering, University of California, Irvine, Irvine, CA, 92697, USA

Sarah C. Finkeldei

Department of Chemical and Biomolecular Engineering, University of California, Irvine, Irvine, CA, 92697, USA

Sarah C. Finkeldei

Institute of Chemistry and Biochemistry, Freie Universität Berlin, 14195, Berlin, Germany

Nina Huittinen

You can also search for this author in PubMed Google Scholar

L.O. formal analysis, investigation, methodology, data curation, writing, review and editing; R.H. formal analysis, data curation, methodology, review and editing; S.S.A.A. formal analysis, data curation, methodology, review and editing; S.E.G. formal analysis, data curation, methodology, review and editing; S.C.F. funding acquisition, project administration, resources, supervision, writing, review and editing; N.H. conceptualization, funding acquisition, project administration, resources, supervision, writing, review and editing.

Correspondence to Sarah C. Finkeldei or Nina Huittinen.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Opitz, L., Hübner, R., Shams Aldin Azzam, S. et al. Investigations towards incorporation of Eu3+ and Cm3+ during ZrO2 crystallization in aqueous solution. Sci Rep 13, 12276 (2023). https://doi.org/10.1038/s41598-023-39143-0

Download citation

Received: 14 May 2023

Accepted: 20 July 2023

Published: 28 July 2023

DOI: https://doi.org/10.1038/s41598-023-39143-0

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.