Open Access Article
Viktoria Baumann
*ab,
Karin Popa
b,
Marco Colognab,
Murielle Riveneta and
Olaf Walter
*b
aUniv. Lille, CNRS, Centrale Lille, Univ. Artois, UMR 8181 – UCCS – Unité de Catalyse et Chimie du Solide, F-59000 Lille, France. E-mail: viktoria.baumann@centralelille.fr
bEuropean Commission, Joint Research Centre, Karlsruhe, Germany. E-mail: olaf.walter@ec.europa.eu
First published on 22nd February 2023
We report on the crystallite growth of nanometric NpO2 and UO2 powders. The AnO2 nanoparticles (An = U and Np) were synthesized by hydrothermal decomposition of the corresponding actinide(IV) oxalates. NpO2 powder was isothermally annealed between 950 °C and 1150 °C and UO2 between 650 °C and 1000 °C. The crystallite growth was then followed by high-temperature X-ray diffraction (HT-XRD). The activation energies for the growth of crystallites of UO2 and NpO2 were determined to be 264(26) kJ mol−1 and 442(32) kJ mol−1, respectively, with a growth exponent n = 4. The value of the exponent n and the low activation energy suggest that the crystalline growth is rate-controlled by the mobility of the pores, which migrate by atomic diffusion along the pore surfaces. We could thus estimate the cation self-diffusion coefficient along the surface in UO2, NpO2 and PuO2. While data for surface diffusion coefficients for NpO2 and PuO2 are lacking in the literature, the comparison with literature data for UO2 supports further the hypothesis of a surface diffusion controlled growth mechanism.
The densification rate is a function of particle size and temperature. If the particle size is decreased, the densification rate will increase, and the temperature and duration of the process can be decreased, saving costs and energy.
In the last decade, different syntheses routes of nanocrystalline AnO2 have been presented.5–10 The hydrothermal decomposition of actinide(IV) oxalates to nanocrystalline actinide dioxide powder has shown several advantages11–14 and the potential to decrease drastically the sintering temperature.15,16
In order to control the sintering process and the final microstructure, its associated mechanisms must be well understood. Sintering is accompanied by grain growth and the elimination of pores. Such processes can occur through different mechanisms as surface diffusion, grain boundary diffusion, lattice diffusion or vapor transport.4,17 In crystalline ceramics, grain boundary (GB) and lattice diffusion (L) from grain boundary to pore contribute most to the densification stage, while diffusion from the surface (S) leads to neck-growth but not to densification. The activation energies Q for the diffusion coefficients are typically in the order QL > QGB > QS, and surface diffusion is thus dominating at low temperature, which is one of the reasons why fast-firing techniques employ a rapid heating rate in the low temperatures, to overcome the non-densifying range.18–20
The migration of pores in nuclear fuels is also controlled by a surface diffusion or evaporation–condensation mechanisms.21 In situ-methods can provide valuable information on the mechanisms involved in sintering, grain growth and pore-migration. Environmental scanning electron microscope at high temperature (HT-ESEM) was introduced as an innovative and powerful method to investigate in situ the sintering behaviour of actinide oxide materials.22–25 For example, Clavier et al. used HT-ESEM device to obtain sintering maps of ThO2, as well as an impressive 2-grain scale observation to capture the first stage of sintering.26 Bouëxière et al. measured the crystallite growth of PuO2 nanocrystals in situ with a HT-XRD device, obtaining an activation energy for crystallite growth of 351(5) kJ mol−1.27 In this work we applied the same method to UO2 and NpO2 nanocrystals. The samples were isothermally annealed for 30 hours in the range of 950 °C to 1150 °C for NpO2, and 650 °C to 1000 °C for UO2. During this process, XRD patterns were recorded so that particle growth could be determined over time at the respective temperatures. We calculated activation energies for UO2 and NpO2 and elaborated on the crystallite growth exponent.
![]() | ||
| Fig. 1 TEM images of NpO2 (left) and UO2 (right) obtained by hydrothermal decomposition of the corresponding actinide oxalates. | ||
The neptunium sample was first heated at 950 °C with a heating rate of 10 °C min−1 and held at this temperature for 30 hours, and then annealed at 1050 °C and 1150 °C following the same procedure. 50 diffractograms were recorded at each temperature in the range of 45° < 2Θ < 61°, with each measurement lasting 36 minutes. A similar procedure was performed with uranium sample. However, because of the broader peak at 46.8°, the three XRD patterns were recorded in the range of 43° < 2Θ < 60°. In contrast to NpO2, fresh UO2 nanocrystals were used for each isothermal annealing. In this case, after reaching the desired temperature (which was maintained for 30 hours), the HT-XRD device was cooled to room temperature and the powder was then replaced with new sample before being heated to the next higher temperature. The isothermal measurement for UO2 was performed at 650 °C, 700 °C, 800 °C, 900 °C and 1000 °C.
![]() | (1) |
| Gn = G0n + kt | (2) |
![]() | (3) |
![]() | (4) |
The subsequent plot of ln(k) as a function of 1/T gives the value of −Q/R as the slope of the straight line, which, multiplied by R, corresponds to the activation energy Q of the mechanism for grain growth.
The powder formed after hydrothermal conversion of actinide oxalates had a fluorite-type cubic structure and crystallized in the space group Fm
m (225). A spherical nearly shape was found for both crystallites from the transmission electron micrographs (Fig. 1). The XRD data of the isotherms measured at 950 °C (for NpO2) and at 650 °C (for UO2) are given on the right side of Fig. 2 as a function of time. The time t0 corresponds to the time after the acquisition of the first XRD spectrum (36 minutes) and consequently t1 = 72 min. The patterns of both samples showed a sharpening of the peaks with temperature and time, caused by crystallite growth. The crystallite size, as a function of time for the first 6 hours are presented in Fig. 3 and in Table 1, the full range (30 h) is shown in Fig. S1.†
![]() | ||
| Fig. 3 The average crystallite size of NpO2 (left) and UO2 (right) as a function of time for various temperature for the first 6 hours. The full time interval is shown in Fig. S1.† | ||
| Time [min] | NpO2 | UO2 | ||||||
|---|---|---|---|---|---|---|---|---|
| Crystallite size [nm] | Crystallite size [nm] | |||||||
| 950 °C | 1050 °C | 1150 °C | 650 °C | 700 °C | 800 °C | 900 °C | 1000 °C | |
| 36 | 20.1(0.9) | 41.8(1.7) | 105.9(1.1) | 6.2(0.2) | 6.4(0.1) | 9.0(1.4) | 25.2(2.5) | 34.4(3.0) |
| 72 | 22.9(0.5) | 46.8(0.2) | 112.9(6.5) | 6.2(0.2) | 6.8(0.3) | 14.4(1.1) | 29.1(0.8) | 40.5(1.0) |
| 108 | 24.2(0.6) | 48.7(2.5) | 119.0(3.3) | 6.7(0.5) | 7.3(0.3) | 17.0(0.5) | 31.1(0.2) | 43.1(2.2) |
| 144 | 25.4(0.6) | 52.0(1.6) | 121.6(8.5) | 6.5(0.2) | 7.6(0.5) | 18.3(0.2) | 31.9(0.4) | 41.8(2.1) |
| 180 | 26.8(0.7) | 55.1(1.6) | 127.4(3.5) | 6.9(0.4) | 7.9(0.6) | 18.1(1.9) | 33.4(0.3) | 45.3(0.4) |
| 216 | 27.5(0.3) | 56.3(2.2) | 130.4(6.0) | 6.9(0.4) | 8.1(0.5) | 20.0(0.4) | 33.0(0.5) | 46.1(1.0) |
| 252 | 27.9(0.3) | 59.0(1.0) | 133.6(7.9) | 7.2(0.6) | 8.4(0.7) | 20.9(0.2) | 34.3(0.5) | 46.1(1.0) |
| 288 | 27.7(0.9) | 59.5(2.2) | 137.5(4.3) | 7.1(0.4) | 8.3(0.4) | 21.2(0.4) | 36.1(0.3) | 47.0(1.8) |
| 324 | 29.1(0.8) | 60.7(3.1) | 144.4(4.2) | 7.3(0.4) | 8.5(0.3) | 21.3(0.8) | 36.1(0.3) | 47.8(0.5) |
| 360 | 28.3(1.3) | 64.6(1.8) | 141.4(9.6) | 7.4(0.4) | 8.6(0.3) | 21.9(0.2) | 37.1(0.7) | 47.2(3.5) |
After the initial rapid growth phase in the first hours, growth with slower kinetics was observed at most annealing temperatures, so that the data after 6 hours were not included in the analysis of the rate constant. Miao et al. also observed a similar phenomenon for the annealing of UO2.30 The grain growth exponent n and the growth constant k were determined from eqn (2) applied to the first 6 hours of annealing, by plotting Gn − G0n vs. time at constant temperature and performing a linear regression with a fixed intercept at 0. The exponent n was restricted in the range n = 2 to 4. The value of n = 4 was chosen because it gave the best linear fit of the experimental points as measured by the adjusted R2 values (Fig. S2 and S3†). The grain growth constant k is the slope of the straight line, and is listed in Table 2 together with the results of the linear fit. The activation energy was then obtained from eqn (4) as the slope of the linear fit of the natural logarithm of k as a function of the reciprocal temperature (1/T) (Fig. 4). In the case of UO2 we neglected the data at T = 1000 °C in the calculation of Q because the linear fit was significantly worse than those obtained for the other temperatures (Table 2), which is as well true for all other exponents n. The activation energies for n in the range n = 2 to 4 are given in Fig. S4† and the fitted crystallite size extrapolated to the full experimental range of 30 h in Fig. S5.†
| Temperature [°C] | D0 [nm] | n | k [nm4 h−1] | R2 | |
|---|---|---|---|---|---|
| NpO2 | 950 | 20.1 | 4 | 127 243(4952) |
0.99 |
| 1050 | 41.8 | 4 | 2 319 034(51 642) |
0.99 | |
| 1150 | 105.9 | 4 | 57 871 100(1 661 392) |
0.99 | |
| UO2 | 650 | 6.2 | 4 | 288(13) | 0.98 |
| 700 | 6.4 | 4 | 788(28) | 0.99 | |
| 800 | 9.0 | 4 | 45 438(177) |
0.99 | |
| 900 | 25.2 | 4 | 288 485(1130) |
0.99 | |
| 1000 | 34.4 | 4 | 835 046(64 084) |
0.94 |
![]() | ||
| Fig. 4 Arrhenius plot for the grain growth constant k using the exponent n = 4 for UO2 (black), PuO2 (ref. 27) (blue) and NpO2 (red). | ||
As a result, we obtained activation energies for crystallite growth of 264(26) kJ mol−1 and 442(32) kJ mol−1 for nanocrystalline UO2 and NpO2, respectively. Fig. 4 shows also data for PuO2 (ref. 27) with a recalculated activation energy of 349(6) kJ mol−1.
Here the exponent n = 4 gave the best linear fit for both UO2 and NpO2 in agreement with what observed in PuO2,27 suggesting that in our conditions the growth of AnO2 (An = U, Np, Pu) nanopowder is controlled by the mobility of the pores, which are migrating via a surface diffusion mechanism. This is not surprising, as pore migration by surface diffusion is favoured at small grain size and low temperatures as the conditions investigated here, whereas pore migration by evaporation–condensation typically dominates at larger pore sizes and temperatures.21,31
The calculated activation energies with n = 4 are 264(26) kJ mol−1 for U in UO2, 442(32) kJ mol−1 for Np in NpO2, and 349(6) kJ mol−1 for Pu in PuO2. Thus the activation energy for the growth of AnO2 nanocrystallites (possibly equivalent to the activation energy for cation surface diffusion) increases in the order Np < Pu < U.
The comparison of our results with the activation energy from the literature is not straightforward: data on activation energy for surface diffusion or for grain growth in PuO2 and NpO2 are scarce or lacking, while data on UO2 are existing, but obtained mostly at higher temperatures and very likely larger pore sizes, where different diffusion mechanisms become dominant.
The activation energy for grain growth of UO2 nanocrystallites controlled by pore migration found here (264 kJ mol−1) is in the broad range of activation energies for grain growth reported in the literature, which is very scattered in the 100–600 kJ mol−1 range.21,33–36 A direct comparison is however not sensible because of the different mechanisms taking place for grain growth at different temperature and sizes (grain boundary, volume or evaporation–condensation) and the grain growth exponent is often assumed or measured between 2 and 4. It is important to note that the choice of the exponent n is crucial in determining the activation energy Q: for example, by varying from n = 2 to 4 in our analysis, the associated activation energies can be twice as high (Fig. S4†). This value of 264 kJ mol−1 is instead to be compared with the activation energy for surface diffusion, which was determined from old experiments and more recently molecular dynamic (MD) simulations. Matzke reviewed nine experiments and proposed an activation energy for surface diffusion of 454 kJ mol−1 in the 1200–1700 °C range.32 This surprisingly high value could be a consequence of the contribution of concurrent mechanisms at such temperatures, as grain boundary, volume and evaporation–condensation. Indeed Zhou and Olander performed more sophisticated experiments by isolating the contribution of evaporation–condensation and obtained a much lower value of 300(60) kJ mol−1 in the 1760–2100 °C range, which is in the same range of what found here.37 More recent molecular dynamics simulations for diffusion of U on the surface of nanocrystals or nanopores consistently confirm a value in the 260–320 kJ mol−1 range.38–40
![]() | (5) |
, with D0s the diffusion pre-exponential factor and Qs the activation energy for surface diffusion, δs the thickness of the surface diffusion layer (taken equal to the lattice parameter a, 0.54 nm),38 Ω the atomic volume (4.09 × 10−29 m3), kB the Boltzmann constant, and r the pore radius. Assuming coarsening by grain growth and pore coalescence (r ∼ G), eqn (5) becomes:
![]() | (6) |
![]() | (7) |
![]() | (8) |
| Actinide dioxide | Temperature range (°C) | D0s (cm2 s−1) | Qs (kJ mol−1) | Method | Reference |
|---|---|---|---|---|---|
| UO2 | 650–900 | 8 × 10−2 | 264 | Nanocrystals growth | This work |
| NpO2 | 950–1150 | 4 × 105 | 442 | Nanocrystals growth | This work |
| PuO2 | 820–1000 | 1 × 102 | 349 | Nanocrystals growth | 27 |
| UO2 | 1200–1700 | 5.00 × 105 | 454 | Review (species: UO2, UO3) | 32 |
| UO2 | 1760–2100 | 5.00 × 106 | 301 | Tracer diffusion (specie: UO2) | 37 |
| UO2 | 1760–2100 | <103 | 301 | Tracer diffusion (specie: U4+) | 37 |
| UO2 | 2427–2827 | 4.49 × 10−2 | 257 | MD simulation, nanopores | 39 |
| UO2 | 1747–2597 | 1.91 × 100 | 319 | MD simulation, nanocrystals | 40 |
| UO2 | 1927–2907 | 4.80 × 10−1 | 302 | MD simulation, nanocrystals | 40 |
| UO2 | 2351–2907 | 4.40 × 100 | 270 | MD simulation, nanocrystals | 38 |
The data for UO2 compares very well with the molecular dynamics simulations. The agreement with older experimental data is less good. One of the reason could be that competing diffusion mechanisms are involved at high temperature, as previously mentioned. The large scattering in published experimental data may also be due to the reported significant difficulties in the theoretical interpretation of tracer diffusion experiments.32,37 Such models were based on the assumption that the migrating species are UO2 and UO3 molecules with rotational degrees of freedom. If U4+ ions are assumed as the migrating species, the pre-exponential factor D0 is reduced by more than 3 orders of magnitude (see the arrow in Fig. 5).37
Footnote |
| † Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d3ra00487b |
| This journal is © The Royal Society of Chemistry 2023 |