Controllable sites and high-capacity immobilization of uranium in Nd₂Zr₂O₇ pyrochlore

As potential nuclear waste host matrices, two series of uranium-doped Nd₂Zr₂O₇ nanoparticles were successfully synthesized using an optimized molten salt method in an air atmosphere. Our combined X-ray diffraction, Raman and X-ray absorption fine-structure (XAFS) spectroscopy studies reveal that uranium ions can precisely substitute the Nd site to form an Nd_2–xU_xZr₂O_7+δ (0 ≤ x ≤ 0.2) system and the Zr site to form an Nd₂Zr_2–yU_yO_7+δ (0 ≤ y ≤ 0.4) system without any impurity phase. With increasing U concentration, there is a phase transition from pyrochlore (Fd3m) to defect fluorite (Fm3m) structures in both series of U-doped Nd₂Zr₂O₇. The XAFS analysis indicates that uranium exists in the form of high-valent U⁶⁺ in all samples. To balance the extra charge for substituting Nd³⁺ or Zr⁴⁺ by U⁶⁺, additional oxygen is introduced accompanied by a large structural distortion; however, the Nd₂Zr_1.6U_0.4O_7+δ sample with high U loading (20 mol%) still maintains a regular fluorite structure, indicating the good solubility of the Nd₂Zr₂O₇ host for uranium. This study is, to the best of our knowledge, the first systematic study on U-incorporated Nd₂Zr₂O₇ synthesized via the molten salt method and provides convincing evidence for the feasibility of accurately immobilizing U at specific sites.

1. Introduction

The disposal of actinide-rich nuclear waste is a key issue for the development of sustainable nuclear energy (Ewing, 2011; Abdou et al., 2018; Finkeldei et al., 2020; Liu et al., 2021). The search for high-capacity, low-cost and radiation-tolerant materials that can be used as host materials for nuclear wastes has been a major research interest for a long time. Pyrochlore oxides have received considerable attention in immobilizing actinides owing to their desirable physical-chemical properties, including their high thermal stability, high chemical durability, and strong resistance to radiation damage (Abdou et al., 2018; Fuentes et al., 2018; Finkeldei et al., 2020). Pyrochlore, ideally A₂B₂O₇, can accommodate a variety of chemical compositions, in which typically a trivalent cation occupies the A sites coordinated by eight oxygen atoms and tetravalent cation occupies the B sites coordinated by six oxygen atoms (Marlton et al., 2021; Talanov & Talanov, 2021). Among them, zirconate pyrochlore is the most important class of host materials due to its remarkable resistance to amorphization under ion beam irradiation (Lu, Shu, Chen et al., 2018; Lu, Shu, Shao et al., 2018; Lu et al., 2019; Shu et al., 2017). Previous studies have demonstrated that actinides can be effectively incorporated into the zirconate pyrochlore matrix at both the A and B sites (Blackburn et al., 2021; Jafar et al., 2015; Lu et al., 2017; Shu et al., 2016; Wang et al., 2021), but many impurities, such as ZrO₂, U₃O₈, ThO₂ and other complex uranium oxides, were introduced at high actinide concentrations, and the stability of the pyrochlore was reduced. Moreover, the immobilization capacity, local structure and valence state of the incorporated actinides (i.e. their speciation) have not been fully explored in previous reports. Realizing the rational manipulation of the phase structure and substituting sites in zirconate pyrochlore, as well as developing low-cost and sustainable synthetic methods, remains a tremendous challenge, which has practical significance in the application of immobilizing actinides.

The conventional method for the preparation of zirconate pyrochlore is a solid-state reaction, which normally results in inhomogeneity and impurities in the final phases (Shu et al., 2016). To improve the chemical homogeneity and phase accuracy of the system, several wet-chemistry routes have been employed for the preparation of these zirconium-based pyrochlore materials, which involves the mixing of reactant precursors and fast mass transport at the molecular level, resulting in better compositional homogeneity and controlled stoichiometric ratio of the obtained samples. However, subsequent sintering is generally carried out at temperatures above 1400°C for more than 48 h, which is of high cost and not suitable for large-scale waste disposal (Lu, Shu, Chen et al., 2018). Principally, molten salt synthesis is a potential methodology that can lower the reaction temperature due to faster mass transport. In previous work (Mao et al., 2009; Gupta & Mao, 2021), La₂Zr₂O₇ pyrochlore ceramic was synthesized via the molten salt method by heating mixtures at 650°C for 6 h. However, until now, only limited work has been reported on actinide-incorporated pyrochlore via the molten salt method (Abdou et al., 2018, 2019; Wang et al., 2020). Neodymium zirconate has also been proposed as a potential ceramic matrix for the immobilization of nuclear waste; therefore, in this work, we chose Nd₂Zr₂O₇ as a suitable host material for studying uranium immobilization via a molten salt process. We successfully synthesized two series of Nd₂Zr₂O₇ (NZO) nanoparticles doped with different U concentrations by an optimized molten salt method and achieved precise control of doping at either Nd or Zr sites to form Nd_2–xU_xZr₂O_7+δ or Nd₂Zr_2–yU_yO_7+δ, respectively. X-ray diffraction, Raman and transmission electron microscopy were used to characterize the crystalline structure, phase transition and morphology in U-doped NZO nanoparticles. XAFS was used to obtain detailed information on the oxidation state and local environments of metal ions.

2. Experimental methods and characterization

2.1. Reagents

The reagents used in this experiment include Nd(NO₃)₃·6H₂O (99.9%), ZrO(NO₃)₂·2H₂O (99.5%), UO₂(NO₃)₂·6H₂O, NH₄OH (28%), NaNO₃ (99.9%) and KNO₃ (99.9%). All reagents are used as received without further purification.

2.2. Material synthesis

All samples were synthesized by an innovative molten salt method (Mao et al., 2009). In the first step, the single-source composite precursor material was synthesized by the coprecipitation method. The original reactants contained Nd(NO₃)₃·6H₂O, ZrO(NO₃)₂·2H₂O and UO₂(NO₃)₂·6H₂O. The reactants were measured stoichiometrically as (5−α) mmol Nd(NO₃)₃·6H₂O, 5 mmol ZrO(NO₃)₂·2H₂O and α mmol UO₂(NO₃)₂·6H₂O corresponding to A site uranium doping, and 5 mmol Nd(NO₃)₃·6H₂O, (5−β) mmol ZrO(NO₃)₂·2H₂O and β mmol UO₂(NO₃)₂·6H₂O corresponding to B site uranium doping. All reagents were dissolved in 200 ml deionized water and stirred at 400 rev min⁻¹. Then, 200 ml of dilute ammonia [NH₄OH:H₂O = 1:9 (v/v)] solution was dropped into the solution, forming a flocculent and exhibiting stratification. The stabilized precipitate was centrifuged out and dried in an oven at 80°C for 24 h to obtain a composite homogeneous composite precursor (2−x)Nd(OH)₃·2ZrO(OH)₂·xUO₂(OH)₂·nH₂O and 2Nd(OH)₃·(2−y)ZrO(OH)₂·yUO₂(OH)₂·nH₂O]. In the second step, NZO and uranium-doped NZO were synthesized by the molten salt method. 0.35 g of the precursor, 2.55 g of NaNO₃ and 3.03 g of KNO₃ were placed into an agate mortar and ground for 30 minutes. The fine mixture was then placed in a corundum crucible and annealed at 750°C for 6 h with a heating rate of 5°C min⁻¹. After natural cooling, the mixture was soaked and washed with deionized water five times to remove the residual molten salt. Finally, the sample was dried at 80°C for 12 h to obtain the product. By controlling the stoichiometric ratio of the initial reactants, NZO and U-substituted NZO materials were obtained, including compounds NZO, Nd_1.9U_0.1Zr₂O_7+δ (A−0.1), Nd_1.8U_0.2Zr₂O_7+δ (A−0.2), Nd_1.6U_0.4Zr₂O_7+δ (A − 0.4), Nd₂Zr_1.9U_0.1O_7+δ (B−0.1), Nd₂Zr_1.8U_0.2O_7+δ (B−0.2) and Nd₂Zr_1.6U_0.4O_7+δ (B−0.4).

2.3. Characterization

The crystal structure was analysed by powder X-ray diffraction spectroscopy (XRD) on a Bruker D8 advanced AXS diffractometer with Cu Kα (λ = 1.54178 Å) irradiation and micro-Raman spectroscopy on a Horiba LabRAM HR Evolution with laser wavelength of 633 nm. The lattice parameters were obtained by Rietveld refinement using GSAS II software (version 4300) (Toby & Von Dreele, 2013). The ideal pyrochlore structure model was used, assuming the 8a and 48f sites fully occupied. The morphology was characterized by scanning electron microscopy (SEM; Zeiss LEO 1530VP) and transmission electron microscopy (TEM; FEI Tecnai G2 F20 S-TWIN). The X-ray absorption spectra experiments were performed at BL14W1 in the Shanghai Synchrotron Radiation Facility (SSRF) and BL16A1 in the National Synchrotron Radiation Research Center (NSRRC). The storage ring of SSRF was operated at 3.5 GeV with a maximum current of 260 mA. Using a Si (111) double-crystal monochromator, the data collection was carried out in transmission mode for the Zr K-edge with an ionization chamber and in fluorescence mode for the U L₃-edge with a Lytle detector. The storage ring of NSRRC was operated at 1.5 GeV with a maximum current of 360 mA. Using a double-crystal monochromator, the sample was placed in an argon atmosphere, and the data collection was carried out in fluorescence mode for the Nd L₃-edge with a Lytle detector. All data were collected at ambient temperature. The XAFS data were analysed using the program Demeter (Ravel & Newville, 2005). For extended X-ray absorption fine-structure (EXAFS) data fitting of Zr-K, the Fourier transform window of the k space was of the sine type in the range 3.0–12.0 Å, and the sine type of the R space was in the range 1.0–4.0 Å. The amplitude reduction factor #!S₀ ² was fixed as 0.75 during EXAFS fitting, based on the fitting result for the known structure of Zr foil. For EXAFS data fitting of U-L₃, the Fourier transform window for k space was of the sine type in the range 2.0–8.5 Å and sine type for R space in the range 1.0–5.0 Å. The amplitude reduction factor #!S₀ ² was set to 0.9 (Zhang et al., 2018).

3. Results and discussion

3.1. X-ray diffraction

The crystal structures of pure Nd₂Zr₂O₇ (NZO) and U-doped NZO materials synthesized by the optimized molten salt method were characterized by powder X-ray diffraction. As displayed in Fig. 1(a), the diffraction peaks of the as-prepared NZO sample are in good agreement with the reflections of a typical pyrochlore structure (space group ; PDF: 78-1617). However, the superlattice peaks of (111), (311), (331) and (511) for the NZO sample are relatively weak, indicating relatively poor crystallinity and the possible existence of defective fluorite structures (space group: ) (Hagiwara et al., 2019). Rietveld refinement of the XRD pattern demonstrates the purity of the NZO sample with a pyrochlore structure [Fig. 1(b)], and the parameters are summarized in Table S1 of the supporting information.

Figure 1 Structural characterization. (a) XRD profiles of Nd2Zr2O7, Nd2–xUxZr2O7+δ (A−0.1, A−0.2, A−0.4) and Nd2Zr2–yUyO7+δ (B−0.1, B−0.2, B−0.4). (b) Rietveld refinement of XRD profile for Nd2Zr2O7 (NZO). (c) Evolution of the fitted lattice parameters with the U concentrations. Raman spectra of (d) Nd2–xUxZr2O7+δ samples and (e) Nd2Zr2–yUyO7+δ samples.

The U-doped NZO samples have a similar XRD pattern as the pyrochlore structure, supporting the introduction of uranium into A or B sites to generate Nd_2–xU_xZr₂O_7+δ or Nd₂Zr_2–yU_yO_7+δ, respectively (Gupta et al., 2016). Compared with that of NZO, the diffraction patterns of all U-doped NZO samples show increasing broadening with increasing uranium substitution level. Such broadening XRD patterns could be ascribed to an increase in lattice distortions induced by the different effective ionic radii of the U (0.73 Å for six-fold coordination and 0.86 Å for eight-fold coordination), Nd (1.11 Å for eight-fold coordination) and Zr (0.72 Å for six-fold coordination) ions (Wang et al., 2021; Ziolkowski, 1985). This may even lead to possible phase transitions in the sample with high U content. In the Nd_1.6U_0.4Zr₂O_7+δ (A−0.4) sample, as marked in Fig. 1(a), second-phase UO₃ (PDF: 45-0856) was observed, and the diffraction peaks in the high-angle area split, indicating the existence of other oxides, such as Nd_0.5Zr_0.5O_1.75 (PDF: 78-1289) or Nd_0.2Zr_0.8O_1.9 (PDF: 78-1301). As shown in Fig. 1(c), we clearly observed an opposite trend in the lattice parameters, which shrink in Nd_2–xU_xZr₂O_7+δ samples but expand in Nd₂Zr_2–yU_yO_7+δ samples with increasing uranium content. Considering that the effective radius of U⁶⁺ ions is between that value of Nd³⁺ and Zr⁴⁺ ions, the shrinkage or expansion of lattice parameters is easy to understand, which supports that uranium ions can precisely replace either Nd or Zr sites in NZO hosts, as expected. Quantitative structural parameters are summarized in Table S1, and the x–O_48f positional parameters in Nd_2–xU_xZr₂O_7+δ and Nd₂Zr_2–yU_yO_7+δ increase regularly with increasing uranium content (except for A−0.4), implying that the crystal structure undergoes a phase transition from the pyrochlore phase to the defect fluorite phase due to the distortion of the ZrO₆ octahedron (Wang et al., 2018). However, A−0.4 does not conform to the rule due to the appearance of an impurity phase and phase separation, indicating that the introduction of uranium at the A site by the molten salt method has reached the limit. Clearly, it is much easier for uranium ions to accommodate at Zr sites than at Nd sites. For this type of substitution, no significant phase impurities can be observed. The transition from ordered to dis­ordered pyrochlore phase and defect fluorite phase cannot be well distinguished by XRD alone. However, these results can be supplemented by Raman spectroscopy.

3.2. Raman spectroscopy

The sensitivity of Raman spectroscopy to oxygen displacements makes it widely used to distinguish pyrochlore and fluorite structures (Mączka et al., 2009; Pokhrel et al., 2016; Sanjuán et al., 2011). In fluorite (), only one Raman active mode is expected. In a perfect pyrochlore (), there are six Raman active vibrational modes: T = A_1g + E_g + 4F_2g (Turner et al., 2017), in which A_1g belongs to O–B–O bending vibration (∼520 cm⁻¹), E_g is mainly attributed to the O–B–O bending vibration and B–O stretching vibration (∼300 cm⁻¹), and F_2g is the O–B–O₆ octahedral vibration and O–A–O_48f bending vibration (∼300 cm⁻¹, 400 cm⁻¹, 500 cm⁻¹, 600 cm⁻¹) (Niu et al., 2019; Park et al., 2018; Zhang et al., 2010). As shown in Figs. 1(d) and 1(e), the Raman spectra of NZO were in good agreement with the peaks of the pyrochlore phase. With increasing U doping amount, the intensity of the Raman spectra in Nd_2–xU_xZr₂O_7+δ and Nd₂Zr_2–yU_yO_7+δ begins to weaken and shifts to low frequency, which is consistent with Maya Abdou's result (Abdou et al., 2018). These changes reflect the increase in the structural disorder degree, especially for anions (Zhang et al., 2014). In particular, the peak intensity decreases to its lowest at A−0.4 and B−0.4, suggesting that the pyrochlore phase was gradually transformed into a disordered fluorite phase. In addition, A−0.1 and B−0.1 exhibit three additional Raman active sites (∼720 cm⁻¹, 770 cm⁻¹, 820 cm⁻¹), which may be related to the distorted B–O₆ octahedron induced by the incorporated U cations (Garg et al., 2008; Tang et al., 2020). Similar modes were found in other uranium compounds (Gregg et al., 2013; Manara & Renker, 2003; Palacios & Taylor, 2000; Abdou et al., 2018; Tang et al., 2020), such as U₃O₈, UO₃ and UO₂²⁺ clusters. In U₃O₈ and UO₃, 700 cm⁻¹ is attributed to the stretching of U-O, and 760 cm⁻¹ is related to the symmetric stretching of O=U=O, despite the slightly lower frequency offset. Through semiquantitative analysis of the Raman peak intensity at 700–800 cm⁻¹, we can obtain the corresponding relationship between the doping level of uranium and the structural transformation. The intensity of the peaks at 700–800 cm⁻¹ increases significantly with increasing U substitution, which indicates that uranium ions can be successfully introduced into the parent phase and form a stable conformation. On the other hand, compared with that of A−0.2, the three Raman active sites (∼720 cm⁻¹, 770 cm⁻¹, 820 cm⁻¹) are well maintained in the B−0.2 sample, implying that the structural incorporation of uranium ions at the B site is better than that at the A site. When the U concentration reaches 20 mol% of the substituted site, both A−0.4 and B−0.4 show significant structural evolution, which is consistent with the XRD results (Gupta et al., 2016).

3.3. Morphological characterization

The morphology of all samples was characterized by SEM and TEM. The elemental distributions (Figure S2) and molar ratios (Table S2) of all samples were detected by SEM-EDS, which indicated good elemental uniformity in the parent material. As shown in Fig. 2(a), the NZO sample has a polyhedral structure with dimensions of ∼20 nm. This good uniformity is attributed to the fast mass transfer of the reaction components at the atomic scale, which benefits from the optimized molten salt method (Dargaud, Cormier et al., 2010; Liu et al., 2013).

Figure 2 TEM images and HR-TEM results of U-doped NZO samples. Panels (a, e) refer to the NZO samples; (b, f), (c, g), (d, h) refer to the A−0.1, A−0.2 and A−0.4 samples, respectively; and (i, l), (j, m) and (k, n) refer to the B−0.1, B−0.2 and B−0.4 samples, respectively; additionally, the illustration shows the corresponding SAED result.

Compared with the morphology of the NZO parent phase, the characteristics of the U-doped NZO samples show obvious changes with increasing U concentration. The particle sizes of A−0.1 and B−0.1 are significantly larger than that of NZO. Such differences reflect that the incorporation of uranium affects the initial growth mode, and a higher content of U may further lead to disordered structures and phase transitions, as confirmed by XRD and Raman results. The crystal grains of A−0.2 evolve into a spherical shape, while B−0.2 still maintains a uniform particle size distribution with a cubic structure. When the uranium content continues to increase, the particles in A−0.4 begin to aggregate and grow into irregular masses, while the particles in B−0.4 have a good dispersion despite evolving into irregular shapes, supporting that U doping at A sites is much more different. Importantly, the inset in Fig. 2(d) shows a bimodal size distribution, indicating a kind of phase separation in A−0.4, consistent with the XRD result. Well defined lattice fringes were detected by high-resolution transmission electron microscopy (HRTEM), and the lattice fringe spacing corresponding to the (222) planes of Nd_2–xU_xZr₂O_7+δ and Nd₂Zr_2–yU_yO_7+δ have opposite change trends, which is consistent with the XRD results. The inset selected-area electron diffraction (SAED) patterns of all samples show clear polycrystalline diffraction rings, corresponding to the three strong peaks of the pyrochlore or defect fluorite phase. In addition, obvious defects can be observed in the HRTEM for the A−0.4 and B−0.4 samples [Figs. 2(h), 2(n)].

3.4. Local electronic and crystal structure

X-ray absorption fine structure (XAFS) is highly sensitive to the local electronic and crystal structures and therefore has been widely used in actinide compounds (Sun et al., 2019; Zhang et al., 2018; Zhang, Su et al., 2016; Zhang, Zhou et al., 2016; Hu et al., 1998). Fig. 3(a) shows the X-ray absorption near-edge structure (XANES) spectra of the U L₃-edge for all U-doped NZO samples, together with spectra of UO₂ and γ-UO₃ as a U⁴⁺ reference and a U⁶⁺ reference, respectively. The energy position of U L₃-edge XANES for U-doped NZO samples is basically the same as that of the γ-UO₃ standard, indicating a U⁶⁺ state in all samples. However, to achieve charge balance, the introduction of high-valent U for substituting Nd³⁺ or Zr⁴⁺ sites should induce a large number of cation vacancies or additional oxygen, along with higher structural disorder (Abdou et al., 2018; Finkeldei et al., 2020). It is important to note that the XANES spectral features can be affected by the valence state, local symmetry and structural disorder. Previous works reported that the energy position of spectral features is affected by structural disorder with shifts of approximately ±0.5–1.0 eV (Conradson et al., 2003, 2004). The high structural disorder around high-valent U may explain the wide main peak of U L₃-edge XANES and the reason that no significant differences were observed between the Nd_2–xU_xZr₂O_7+δ and Nd₂Zr_2–yU_yO_7+δ samples in Fig. 3(a). To detect the bonding between uranium and oxygen ligands in U-doped NZO, we performed an extended X-ray absorption fine-structure (EXAFS) analysis of the U L₃-edge by using Fourier transform (FT). Figure S3 shows the k³-weighted χ(k) signal in the U-doped NZO sample and its corresponding U L₃-edge FT. The broad peak at ∼1.5 Å in the R space arising from the U–O coordination shell indicates that there is a high oxygen disorder around the U site. Table S3 shows that this disorder (σ²) reaches its highest value at A−0.4 and B−0.4. At the same time, the coordination number decreases as the uranium concentration increases. It is worth noting that, for the U–O path, Nd_2–xU_xZr₂O_7+δ has a shorter bond length than Nd₂Zr_2–yU_yO_7+δ. This difference implies that the molten salt method realizes the doping of uranium at the A and B sites of the NZO matrix material and produces a slight difference in the oxygen coordination of uranium.

Figure 3 XAFS data of the NZO and U-doped NZO samples. (a) U-L3, (b) Nd-L3 and (c) Zr K-edge XANES spectra of U-doped NZO samples, (d) Fourier transform (FT) of Zr K-edge EXAFS data and their corresponding fits in R space.

To assess the structural changes of Nd and Zr affected by dopant U ions in U-doped NZO samples, Nd-L₃ and Zr-K XANES measurements were also performed. The Nd L₃-edge XANES spectra of all samples exhibit a similar shape, as shown in Fig. 3(b), and the white line energy position is very close to that of the Nd(NO₃)₃ reference, demonstrating the 3+ valence state of the Nd ions. In the case of an increase in the valence by one, the white line maximum would shift by 7–8 eV to higher energies (Lissner et al., 1994; Hu et al., 1997). In a typical pyrochlore, Nd coordinates with eight oxygens to form a hexahedral structure, and it contains two different crystallographic oxygen sites: O_8b and O_48f. Among them, O_48f has special properties concerning the adjustment of oxygen disorder and structural distortion in pyrochlore (Hagiwara et al., 2019; Martin et al., 2009). Previous studies pointed out that the introduction of high-valent ions will cause the 8a oxygen vacancies to be partially filled and cause the anion disorder of Gd₂Zr_xTi_2–xO_7+δ pyrochlore (Nästren et al., 2009). In the Nd L₃-edge XANES, the strong white line originates from a dipole-allowed transition from core level 2p to 5d, and the spectral profile is related to the unoccupied Nd 5d states (Rao et al., 1983). As observed in Figure S4, the intensity of the white line at the Nd L₃-edge XANES increases gradually with increasing uranium concentration. This can be mainly assigned to the strong covalent mixing of broad U 6d with O 2p competing with Nd 5d–O 2p mixing. Furthermore, U doping leads to an increase in the symmetry of the Nd₂Zr₂O₇ system and, in turn, a narrower 5d band. Fig. 3(c) shows the Zr K-edge XANES spectra of the NZO and U-doped NZO samples, which can be characterized by an obvious pre-edge structure (labelled by a) and a main peak split into two well resolved features (labelled by b and c). Peaks a and b/c correspond to the electronic transition from the 1s to 4d orbital and the dipole-allowed transition from the 1s to 5p orbital, respectively (Komyoji et al., 1992; Popov et al., 2019). These features of all samples are close to those in the pyrochlore phase rather than in the defect-fluorite phase (Nandi et al., 2020). In Zr K-edge XANES, the pre-edge feature a for all samples is significant, which indicates tetrahedral distortion of the ZrO₆ coordination polyhedron (Dargaud, Calas et al., 2010). The c/b intensity ratios of all samples are similar to that of pyrochlore-type Gd₂Zr₂O₇ with six-fold coordinated Zr⁴⁺. The ratio would rather deviate for Zr^4+δ (NZO) compounds with seven- and eight-fold coordination (Komyoji et al., 1992; Galoisy et al., 1999; Li et al., 1993). Deeper knowledge of the local structure around Zr can be obtained from the quantitative analysis of EXAFS spectra. As shown in Fig. 3(d), all samples show two main coordination peaks, in which the first peak corresponds to six O_48f atoms at approximately 2.2 Å, and the second peak corresponds to six Zr ions and six Nd ions. The quantitative results are shown in Table S4.

The Zr K-edge EXAFS spectra are sensitive to changes in the local environment with increasing uranium content in the Nd_2–xU_xZr₂O_7+δ and Nd₂Zr_2–yU_yO_7+δ samples. A noticeable change in the coordination number (CN) is observed for Zr–O_48f bonds during uranium doping. Compared with that in NZO, the CN values of Zr–O increased slightly in the A−0.1 and B−0.1 samples, which may be caused by the introduction of anions to maintain charge balance and preferentially fill the unoccupied 8a site in the pyrochlore structure (Hess et al., 2002; Nästren et al., 2009). Since Zr is located at the 16c site (adjacent to the 8a site), the addition of Zr–O CN has little effect on the polyhedral structure of Zr–O_48f. However, the CN of Zr–O decreased significantly when keeping increasing the U content in Nd_2–xU_xZr₂O_7+δ and Nd₂Zr_2–yU_yO_7+δ. In this case, more uranium is introduced into pyrochlore, which causes the Zr–O octahedron to be distorted, and the structure has more oxygen disorder characteristics by the movement of O_48f and the formation of anionic Frenkel defects (Popov et al., 2019). In particular, this disorder is maximized at the highest doping concentration, which is consistent with the extremely low characteristic peak of pyrochlore in the Raman results. Through conventional Zr K-edge EXAFS analysis, the pyrochlore and defect fluorite phases cannot be accurately distinguished. For the coordination of Zr–Zr and Zr–Nd, CN decreases continuously with increasing uranium amount in the Nd_2–xU_xZr₂O_7+δ and Nd₂Zr_2–yU_yO_7+δ systems, and the bond lengths of Zr–Zr and Zr–La increase regularly. This indicates that the movement of pyrochlore cation sites and the deformation of the pyrochlore lattice lead to the formation of a disordered system.

4. Conclusion

In this work, we systematically studied the uranium solubility, crystalline phase, valence state and structural evolution in a Nd₂Zr₂O₇ pyrochlore synthesized by a modified molten salt method. According to the XRD results, uranium ions can be precisely confined at either the Nd or Zr sites in Nd₂Zr₂O_7, and the solubility can reach up to a 20 mol% substitution at the Zr sites and a 10 mol% substitution at the Nd sites without any impurity phase. The Raman analysis shows that there is a phase transition from ordered pyrochlore to defective fluorite with increasing U concentration. The TEM images show that all U-doped Nd₂Zr₂O₇ samples have good crystallinity and dispersity except for Nd_1.6U_0.4Zr₂O_7+δ. The U L₃-edge XANES spectra indicate that uranium exists in the form of high-valent U⁶⁺ in all samples. To balance the extra charge from substituting Nd³⁺ or Zr⁴⁺ by high-valent U⁶⁺, additional oxygen is introduced along with a large structural distortion, as revealed by the XAFS analysis. This systematic work sheds light on understanding actinide immobilization, which will accelerate the development of advanced technology in nuclear waste storage materials.