3.1.1. Structural and textural parameters.
Figure 1 shows the X-ray diffractograms of the four catalysts prepared. The three Pr-based formulations chosen (based on previous results) for praseodymium-based catalysts exhibit a single cubic fluorite structure [
14]. The four typical reflections corresponding to the (111), (200), (220), and (311) planes are observed at 28.55, 33.05, 47.50, and 56.40° of 2θ values, respectively. On the other hand, pure ceria exhibits the expected pattern as well.
Additionally, the lattice parameter (
a) can be obtained for each catalyst by using Bragg’s law (and determining the d parameter), while Scherrer’s equation is employed for the average crystal size estimation. All these results are summarized in
Table 1. The values slightly increase by considering a progressive substitution of cerium by praseodymium into the cationic sublattice providing evidence of a higher presence of the Pr
3+ cation reaching the pure praseodymium oxide formulation limit [
7]. On the other hand, quite similar average crystal sizes can be observed for the four formulations. Conversely, a lowering in the BET surface area is seen inside the series and is quite noticeable for pure praseodymia (6 m
2/g). A previous work [
14] reported a negative effect of very high Pr loadings onto ceria (concerning BET surface area and pore development) when the Direct Calcination method was employed. In fact, recent works by Frizon et al [
20] and Fahed et al [
21] reported that the insertion of high Pr loadings in mixed oxides (specifically CZ mixed oxides) significantly decreased the specific surface area. The observation of the non-correlation between S
BET and average crystal size (in spherical approximation; r(nm) = 3/(ρ · S
BET)), where ρ is the compound density) was attributed to the fact that the particles were not spherical and the aggregation phenomena occurred limiting the pore size (see
Table 1). High Pr loading can decrease the inter-granular porosity in line with the drop of the porous volume observed from the data in
Table 1.
Raman spectroscopy is a highly suitable technique for the structural characterization of cerium-based oxides, providing additional structural information beyond that offered by the XRD technique. The XRD patterns provide information mostly about the positions of large cations in the crystal lattice (such as cerium). Still, this technique offers very low sensitivity to light elements (such as oxygen) [
12]. Raman spectra of ceria revealed the
F2g band at 464 cm
-1, which might be attributed to a symmetrical breathing mode of the oxygen ions combined with the tetravalent cations [
8,
9]. The three Pr-containing catalysts show the F
2g vibration mode of the fluorite structure (see
Figure 2), where the band position is moved to lower wavenumber values by the addition of Pr content onto ceria, as listed in
Table 1, where this decrease can be due to the existence of the trivalent cation (rather Pr
3+ than Ce
3+), i.e., a weaker force constant of the corresponding cation-anion bonds. This interpretation is also supported by the observed mode at a wavenumber of about 570 cm
-1, which is attributed to a defect mode caused by oxygen vacancies as mentioned in a previous publication [
24]. The relatively high intensity of this mode under atmospheric conditions indicates the presence of the reduced cations (rather being Pr
3+ than Ce
3+). An additional mode around 195 cm
-1, whose evolution can be correlated with that of the band at 570 cm
-1, can be attributed to the other asymmetric vibration caused by the formation of oxygen vacancies [
25]. The I
560/I
F2g ratio obtained from the normalized spectra correlates with the Pr content, as can be seen in
Table 1. This ratio can be considered as an indicator (or descriptor) of the abundance of oxygen vacancies caused by the substitution of Ce atoms by Pr atoms in the mixed oxides until achieving pure praseodymium oxide [
25].
In order to explore if the catalysts present varied morphologies dependent on the molar composition of cerium and praseodymium, the high-resolution TEM images of the catalysts investigated are illustrated in
Figure 3. The presence of aggregates of multiple nanoparticles can be generally observed: a dense substance made of nanocrystalline domains that have been sintered exhibits a significant number of likely disordered grain boundaries. Although the individual particles are not visible in
Figure 3(a
1-d
1), (this might be due to the route of synthesis used, as reported in a previous study [
26]), there is strong evidence of crystallinity within the domains, as seen by the visible lattice fringes (
Figure 3 b
2-d
2) and even clearer in
Figure 3(a
1-a
4), corresponding to the pure ceria catalyst (analyzed for comparison purposes). The interplanar distances of the lattice fringes observed are compatible with the fluorite structure of the catalysts, showing (111) fringes (close to 0.31 nm) [22,23] and (220) facets with an interplanar distance of 0.19 nm as shown in
Figure 3b
2 and c
2 [
29], which confirms the XRD results, where the reflections corresponding to both planes are the most intense ones detected. Actually, the d value obtained for the (111) facets, measured for the three catalysts, exhibits a slight and progressive increase in the corresponding value (0.30, 0.301, and 0.31 nm, respectively) as the Pr loading enhances from Ce
0.3Pr
0.7O
2 to PrO
2, in complete agreement with
a (lattice parameter) values measured by XRD where this value is also increased. It should be noted that the crystal structure of ceria and praseodymia show large similarities [
30]. Therefore, phase segregation is difficult to infer (if any), and according to these results, there is no evidence of such phenomena [
16].
The representative HAADF images
Figure 3 (a
3-d
3) of CeO
2, PrO
2, Ce
0.2Pr
0.8O
2, and Ce
0.3Pr
0.7O
2, respectively show the appearance of spherical agglomerates in the mixed oxide, these images corresponded well with the energy dispersive spectroscopy (EDS) element mapping images of Ce (a
4-c
5 and d
5), Pr ( b
4-c
4 and d
4) and O ( a
5 -b
5) showing relatively high-density distributions, where there is a homogenous distribution of cerium and praseodymium on the nanometric scale over the Ce
0.2Pr
0.8O
2 and Ce
0.3Pr
0.7O
2 catalysts, which is also consistent with the TEM and XRD results.
Surface-sensitive X-ray photoelectron spectroscopy (XPS) has been used to study the surface chemical state of mixed oxide particles.
Table 2 shows some representative parameters obtained from the XPS spectra analysis (see
Figure S1 in the Supplementary Information).
The obtained results suggest that the Ce/Pr surface atomic ratios were consistently lower than theoretical or nominal values in ceria-praseodymia samples. This indicates that there was a higher concentration of praseodymium on the particle surface (or the periphery of the particles) than that of cerium (thus exhibiting Pr surface enrichment) [
12]. When conducting estimations, the Ce 3d and Pr 3d levels were utilized, both of which do not exhibit significantly different binding energies. The Ce 3d level has a binding energy of approximately 870 eV, while the Pr 3d level has a binding energy of around 970 eV. These values indicate the amount of energy required to remove an electron from the respective 3d orbitals of cerium (Ce) and praseodymium (Pr) atoms [
31]. Overall, the segregation of dopant cations at the surface and domain boundaries of ceria is an important aspect to consider when studying or utilizing ceria-based materials, as it can have significant implications for their properties and performance in various applications [21,28,29].
Pr cations (Pr
3+ and Pr
4+) are incorporated into the cerium oxide lattice (and or partially segregated at the particle’s surface), as suggested by comparison of Ce/Pr surface values with the nominal ones. Therefore, the presence of nano-domains with different Ce/Pr compositions (indistinguishable by XRD and/or HR TEM cannot be excluded for the two mixed oxides). Corresponding estimations are compiled in
Table 3.
The expected and well-described “carbonation” process that occurs on solid surfaces during calcination and storage of samples can pose challenges for the analysis and interpretation of Ce
3+ and Pr
3+ cations existence [
14,
34]. This process leads to the accumulation of carbon on the surface of these oxides, as indicated in
Table 2. The presence of carbon species on the surface can make the accurate determination of the Ce
3+ and Pr
3+ cations, (truly involved in the creation of oxygen vacancies) quite difficult. The estimated values of Pr
3+ (%) are always higher than those of Ce
3+ (%), for the two mixed oxides prepared, which are the expected results since Pr is a more reducible cation (due to the higher reduction potential and more labile Pr-O bond) and is more basic, compared to Ce [
1,
30,
31].
3.1.1. Comparison of the oxygen lability of the catalysts under different environments
The redox properties and the corresponding lability of the oxygen species in these formulations were explored by means of two procedures that can be complementary one each other: H
2-TPR and O
2-TPD (to measure the O
2 evolved under an inert atmosphere). For more details about definitions, procedures, and estimations conducted to determine the amount of oxygen vacancies created after the two different procedures, please see the
Supplementary Information.
O
2-TPD under an inert atmosphere (helium) was used to examine the oxygen “delivery” of the catalysts, and in turn, its lability under helium, by following the rate of oxygen desorption in terms of temperature (and its subsequent global quantification, thus obtaining the value of µmol O
2 emitted/g
cat).
Figure 4(a) displays the O
2-TPD profiles of the catalysts, expressed as this mentioned oxygen emission rate versus temperature, and
Table 4 compiles the corresponding quantifications. Even though CeO
2 is a well-known oxygen storage material, O
2 desorption under an inert atmosphere was hardly observed in this case, in line with previous results [
37], as shown in
Figure 4(a). Conversely, the three other samples (the Pr-containing catalysts) yielded three main and large oxygen desorption peaks from 300℃ to 900℃, whose general pattern is very similar (differing mainly in the relative contribution of every peak), which seems to be a function of the Pr-content. Since the two mixed oxides’ general profile is very similar to that obtained for the pure oxide, the discussion about the oxygen mobility can be assessed in a more certain way due to the presumed absence of separate rich Pr and/or Ce domains in the solids, which could probably modify the oxygen emission pattern, thus making distinctions if referred to that of the pure oxide. It is generally expected that the existence of such large domains would display different oxygen labilities. This indirect evidence would be congruent with the XRD and TEM discussion, where a main crystalline phase with a gradual variation of the lattice parameter was observed when going from the Ce
0.3Pr
0.7O
2 sample to PrO
2.
The sharpest peak and the highest values of O
2 emission rates in the low-temperature range are reached by Ce
0.3Pr
0.7O
2, at 363 °C, (this temperature maximum is very similar for the three catalysts), which can be attributed, in principle, to the most surface/reactive oxygen susceptible to being evolved. The profile and magnitude of this first peak, being the sharpest and most intense one of the whole pattern for the two mixed oxides, reveals an accentuated oxygen lability and, in turn, relevant oxygen mobility at very low temperatures, which can be assumed to be very influenced by the BET surfaces, since the magnitude and maximum oxygen emission rate at this very low temperature follows the same trend than that of the BET surface area and pore volume of the solids. Conversely, the second and third peaks are more accentuated for the pure PrO
2, revealing the maximum bulk oxygen mobility and the best data of global oxygen evolved for this formulation. Therefore, the more the Pr loading, the higher the O
2 evolved amounts, which is in line with the conclusions of other authors [
7,
11]
.
For the deepest insight, some quantitative considerations will be provided now. Considering a mean oxygen density of 21.7 µmol O·m
-2 in ceria-based oxides, according to the assumptions taken by Fahed et al [
21] based on other authors [
38], the oxygen emitted relative to the first peak (see values in parentheses of
Table 4) only concerns to surface oxygen species (or a fraction of them) for the two mixed oxides under study. Conversely, for the case of PrO
2, more than the theoretical oxygen surface population (even though it is small due to its very low surface area) can be desorbed in this low range of temperature, proving evidence of the very high oxygen mobility at the surface/subsurface level for this catalyst.
It is interesting to approach the whole oxygen delivery capacity deduced from these experiments (O
2-TPD). If it is bear in mind that the ideal reduction of PrO
2 to Pr
2O
3 would yield a maximum amount of oxygen evolved of 1445 μmol O
2/g
cat (indicating the maximum oxygen released that could be reached by a bulk oxide of this starting formulation, fully oxidized) [
39], the value of 1045 μmol O
2/g
cat for pure praseodymia (even characterized by a very low surface area) indicates very high lability of oxygen and an accentuated mobility of subsurface/bulk oxygen even under a not reducing atmosphere (He).
H
2-TPR was used to assess the reducibility of the selected catalysts as well as determine the amount of oxygen vacancies susceptible to being created in a different reaction atmosphere.
Figure 4(b) depicts the profiles of H
2 consumption obtained by this procedure. It is interesting to compare as deeply as possible both
Figure 4a and
Figure 4b. Firstly, the profiles under H
2 are characterized by similar patterns (in agreement with the previous discussion) and asymmetrical shapes, the higher the Pr content in the catalysts, the sharper and more intense the peak. This has a certain parallelism with the O
2-TPD profiles, where the three different peaks can be interpreted as the “delivery” of more and more internal lattice oxygen, but appearing, now, as a unique peak (but with weak shoulders, under a more aggressive atmosphere, H
2) [
40]. From a very basic and qualitative perspective, it seems that the three previous peaks, detected under He, have been combined in a single asymmetric peak that decays at an early temperature (around 500°C). The corresponding integration of this peak (in
Figure 4b) for every catalyst confirms the increasing reducibility (ability to form oxygen vacancies) as the Pr loading is progressively enhanced in the catalysts’ formulation [
24].
The delta (δ) parameter has been used to indicate that a single oxygen-deficient non-stoichiometric phase is formed. In this sense, the general formula after the corresponding reduction will be written as MO
2-δ or Ce
1-xPr
xO
2-δ, to express the maximum level of non-stoichiometry that has been reached under the two different atmospheres (He or H
2). For example, for PrO
2, its maximum theoretical level of non-stoichiometry would be δ equal to 0.5 [
39,
41], which corresponds to the reduced formula of PrO
2-0.5 or PrO
1.5 (better written as Pr
2O
3 ), where the whole Pr
+4 cations have been reduced to Pr
+3).
The values of δ or non-stoichiometry reached, together with the final formula after the reduction of the different catalysts, under the two atmospheres, are shown in
Table 4. An increasing trend of δ is observed with the increase in the content in praseodymium, reflecting the better reducibility, and in turn, improved oxygen mobility in both atmospheres, since the oxygen amount evolved is very high (and not only affecting the surface oxygen, as previously commented, and explained below in more detail). Besides, as the Pr content is increased, the capacity of oxygen delivery under the two atmospheres becomes closer and closer, thus reflecting superior oxygen lability as well, which is suggested to be of paramount importance for the soot combustion application.
Figure 5 illustrates the correlations of the amounts of oxygen vacancies created under the two atmospheres (estimated from global O
2 emitted and H
2 consumed amounts shown in
Table 4 and following the assumptions explained in the
Supplementary Information) versus the praseodymium content of the catalysts (ranging from pure ceria to pure praseodymia). It is worth noting that, independently of the different textural properties that characterize the catalysts, the main parameter governing the total oxygen release capacity of the oxides is the praseodymium content (approximately linear correlation see
Figure 5), mainly that obtained by the representation of O
2-TPD data, with a value of r
2 of 0.982). Besides, an interesting observation can be extracted from
Figure 5, as the Pr content is higher and higher, the values obtained under the two different atmospheres become closer and closer (as illustrated in
Figure 5), proving the best oxygen mobility and, in turn, the high ability of stabilizing vacancies for PrO
2 (even under mild conditions).
Interestingly, for the case of pure praseodymia, the maximum theoretical level of reducibility that can be reached (measured from both atmospheres) corresponds with the value of the theoretical non-stoichiometric phase of Bevan’s cluster. It is described in the literature [
42] that fluorite-related solids (MO
2-x) can be anion deficient, with the cation sublattice remaining essentially perfect. The defect cluster would correspond to tightly bound vacancies along a body diagonal [111] (as shown in
Figure S2 of the Supplementary Information). The central cation is six-coordinated and is surrounded by six seven-coordinated cations. Since the geometry of these clusters permits four different equivalent orientations, it is easily accommodated in the fluorite structure and accounts for the wide range of stability of the so-called α-phases (MO
1.71-MO
2). This congruency between the experimental estimation and the theoretical defect cluster mentioned proves the large population of bulk oxygen vacancies accommodated in this solid, thus leading to a bulk defect cluster structure and becoming indirect evidence of the large capacity of stabilizing oxygen vacancies into the lattice after mild or reducing atmospheres.