Temperature Dependent Variations of Phonon Interactions in Nanocrystalline Cerium Oxide

The temperature dependent anharmonic behavior of the phonon modes of nanocrystalline CeO 2 was investigated in the temperature range of 80–440K.The anharmonic constants have been derived from the shift in phonon modes fitted to account for the anharmonic contributions as well as the thermal expansion contribution using the high pressure parameters derived from our own high pressure experimental data reported previously. The total anharmonicity has also been estimated from the true anharmonicity as well as quasiharmonic component. In the line-width variation analysis, the cubic anharmonic term was found to dominate the quartic term. Finally, the phonon lifetime also reflected the trend so observed.


Introduction
CeO 2 is a strategic material which has been extensively investigated and has recently been found to be of use in intermediate temperature oxide fuel cells as it easily absorbs and releases oxygen [1].It also shows enhanced ionic conductivity with temperature [2].Hence, the investigation of the temperature dependent behavior of CeO 2 has gained immense interest.Under ambient conditions, it usually crystallizes in the cubic fluorite structure with space group Fm3m [3].Nanocrystalline cerium oxide has applications in the coating technology, for example, wood-coating industry [4].
It is well known that an increase in temperature in any material introduces perturbations in the harmonic potential term, which changes the vibrational properties.Temperature affects the population of the different levels for each normal mode instigating the change of the lattice parameters.Any change of the lattice parameters with temperature is attributed to the anharmonicity of the lattice potential.As an effect of the anharmonic interaction, a nonequilibrium phonon population decays into phonons of lower energy or is scattered by thermal phonons into modes of frequencies shifted from their original positions [5].These anharmonic interactions occur primarily due to the cubic or quartic anharmonicities, resulting in the decaying of an optical phonon into two or three acoustic phonons, respectively, which also decreases their lifetimes.As a consequence, the phonon line position and bandwidth in the Raman spectra of various amorphous or crystalline materials exhibit significant changes with temperature [6][7][8].Further, the finite-size effects in nanocrystals are expected to modify the anharmonicity and the phonon decay times.It has also been reported that the anharmonic constants related to the peak widths and peak positions are higher in nanocrystals than in the bulk, implying a greater degree of anharmonicity in nanocrystals [9].Very few experiments have been carried out to investigate the temperature dependence of Raman mode frequencies and line-widths in nanocrystals [5,[10][11][12].
Still fewer studies have been carried out on the anharmonic behavior of the CeO 2 phonon modes [13][14][15].Popović et al. [13] reported that three contributions are involved in the temperature-induced Raman mode frequency shifts of nano-CeO 2− : the thermal expansion, the anharmonic effect, and the phonon frequency change due to the change of nanoparticle size with temperature (confinement effect).It is natural to expect that in nanomaterials the contribution of four-phonon processes should be larger compared to that of the three-phonon processes [11].However, Spanier et al. [14] reported that phonon coupling is not stronger in the CeO 2−y nanoparticle system than that in the bulk counterpart.Consequently, the aim of this work is to examine the temperatureinduced changes in the vibrational properties of CeO 2 , to identify the contribution of volume expansion and anharmonic intermode coupling to the Raman frequency shift, line-width, and the phonon lifetime in CeO 2 .
In this context, we report the temperature dependent Raman behavior of nanocrystalline CeO 2 and its anharmonic constants derived using our previously published high pressure data [16].We have analyzed the softening of the most dominant mode ( 2 ) in light of the theory of anharmonic phonon-phonon interaction and thermal expansion of the lattice.

Characterization at Ambient Conditions
The CeO 2 powder used in the present study is a commercially prepared high purity powder from Johnson Matthey Rare Earth Products, UK.The powder was used in the present study in as-received state and was not subjected to any pretreatment and, as per the manufacturer's claims, the common metal impurities are at a ppm level.The characterization of this CeO 2 nanopowder under ambient conditions via X-ray diffraction, atomic force microscopy, and Raman spectroscopy has been discussed in detail elsewhere [17].In brief, the crystalline phase at ambient was found to be cubic.The crystallite size of this cubic CeO 2 was found to be 36 nm calculated using Debye-Scherrer relation and about 32 nm from atomic force microscopy.The predominant Raman peak of CeO 2 was obtained at 465 cm −1 , assigned as  2 mode [17].In Ceria, the first order Raman active mode is the triply degenerate symmetrical stretching vibration of the CeO 8 vibrational unit.In bulk CeO 2 , there is a single allowed Raman active triply degenerate  2 mode centered at about 465 cm −1 .This mode can be viewed as a symmetric breathing vibrational mode of the O ions around each cation [18,19].The CeO 2 belongs to space group Fm3m (225).This structure has two BCC unit cells and has two acoustic modes and 1 optical mode.

Temperature Dependent Raman Studies
The Raman spectra were recorded in the temperature range 80-440 K.The spectra were acquired using Jobin Yvon T64000 triple Raman spectrometer with the 514.5 nm line of the Ar ion laser as the source radiation.The sample temperature was varied using a continuous flow liquid nitrogen cryostat from Janis, model ST-500.
The temperature dependent Raman spectra of CeO 2 are given in Figure 1.The spectra show changes in line position, line-width, and intensity with an increase in temperature from 80 K to 440 K.We have focused on temperature dependence of Raman frequency and line-width of the predominant optical  2 phonon in the studied temperature range.With temperature increase, the  2 mode shifts to lower frequencies and the line-width continuously increases.A vertical line on the  2 mode in Figure 1 facilitates the visualization of the extent of shift in said phonon mode.
Hence, the observed results show temperature dependent phonon softening of  2 mode, that is, red shift with progression of temperature.The possible reasons for red shift can be heating and phonon interactions and so forth.Heat causes expansion of the material which leads to frequency decrease.The frequency softening is significant up to a temperature of ∼380 K. Above this temperature, the phonon softening is seen to slow down.It may be mentioned here that, as per the reported Debye temperature of 355 K [20,21], the phonon softening may be said to be slowing down above the Debye temperature.This may be owing to the fact that, above the Debye temperature, the change in specific heat capacity becomes negligible.Figure 2 shows the  2 peak positions observed at 80 K and 440 K which are centered at 468.2 cm −1 and 463.3 cm −1 , respectively.The frequency of this mode can be clearly seen to shift to lower energy side and the line-width increases with temperature.In the studied temperature range, the total quantitative softening of the  2 mode amounts to about 5 cm −1 .These shifts as well as the line-widths estimated from the "peakfit" software are graphically depicted clearly in Figure 3.It may be interesting to note here that similar to the behavior shown by mode softening the line-width increase also shows a lower slope above about 380 K.
The observed phonon softening of the  2 mode is examined from the theory of anharmonic phonon-phonon interaction and thermal expansion of the lattice.As mentioned, the variation of the frequency of the normal modes with temperature at constant pressure arises from two contributions: pure volume contribution (Δ) latt which results from thermal expansion and pure-temperature contribution (Δ) anh which results from cubic and quartic anharmonicities (self-energy shift).The cubic and quartic anharmonicities are described by factors which include contributions from three-and four-phonon anharmonic processes [11].The quasiharmonic  term (Δ) latt arises from the thermal expansion, while puretemperature contribution (Δ) anh due to cubic and quartic anharmonicities is also termed as intrinsic (true) anharmonicity [9,12].Hence, The quasiharmonic or thermal expansion contribution can be written as where  0 (=470.1 cm −1 ) is the harmonic frequency of  2 mode which is obtained by extrapolating the experimental data down to 0 K.The Grüneisen parameter  (=1.28) of the corresponding mode is estimated from our high pressure Raman scattering data reported elsewhere [16] and  is the linear thermal expansion coefficient.The anharmonic contribution (Δ) anh which arises from the phonon-phonon interaction due to the cubic and quartic terms in the interatomic potential can be written as The first term estimates the coupling of an optical phonon to two low-energy phonons (three-phonon process); the second term is due to the coupling of three phonons (fourphonon process), and  and  are the fitting parameters.We have also calculated the total anharmonicity of  2 mode as a function of temperature using the following equation [22]: The first term of the right hand side describes the true anharmonic contribution and the second quasiharmonic term.As mentioned, the resultant frequency shift may be either positive or negative depending on the relative magnitudes of the anharmonic terms in the interatomic potential.At lower temperatures, the main contribution to the softening is due to the positive quartic anharmonicity.
In Figure 4, we have shown the fitted experimental data points (solid circles) by using (1).The fitting parameters obtained for the present case using the above equations are given in Table 1.From the fitting parameters given in the table, it is clear that coefficients  and  both are negative thus explaining the phonon softening.However, coefficient  dominates over  indicating that three-phonon processes or cubic anharmonicity is contributing significantly to the overall softening observed.It is also seen that the true anharmonic term dominates, albeit slightly, over the quasiharmonic term.
A comparison of these coefficients with those obtained by Dohcevic-Mitrovic et al. [15] and Popović et al. [13] reveals that, in CeO 2−y nanocrystals, although a dominance of fourphonon process was observed in the high temperature range by these researchers, however, we have found significant contribution from three-phonon process to the anharmonicity.Our results differ as we have observed a behavior in nanocrystalline CeO 2 which has been reported for polycrystalline CeO 2 [13].It is understood that, in small particles of pure and doped ceria nanocrystals, when size effects have minor impact on Raman modes, four-phonon anharmonic processes prevail [9,15].However, when nanopowdered particles are grown enough, the size effects provoke changes of the anharmonic interactions whereby three-phonon coupling prevails over the four-phonon anharmonic processes.
The line-width of the  2 mode was also studied as a function of temperature in the 80-440 K range, which, as expected, increased with temperature as shown in Figure 3.The FWHM of the  2 mode increases from 4.85 cm −1 at 80 K to 11.43 cm −1 at 440 K. Hence, the line-width increases with temperature by about 6.58 cm −1 .
The experimental peaks of real materials exhibit an intrinsic width which has contributions from factors such as the finite resolution of the spectrometer and impurities/defects which disturb the translation symmetry of the harmonic crystal [22,23] apart from the anharmonic decay of the phonons.The thermal expansion is a manifestation of the lattice anharmonicity that has its own shifts in the frequency, but not on the line-width [23].Hence, considering that other effects due to impurities, defects, and spectrometer are constants for a particular sample, the line-width variation would reflect the anharmonic decay of phonons.The temperature dependent phonon line-width is then written as [24] where Γ 0 is the harmonic line-width independent of temperature and it is caused by defects, including isotopic mixture; the second part of the equation is the damping part induced due to the three-and four-phonon process cubic and quartic anharmonicities, respectively.Figure 5 shows the fitting to the experimental points using (5).A good agreement can be seen between the calculated curve, given by the solid line, and the experimental points.For the present case, the parameters in (5) are Γ 0 = 2.86 cm −1 and  and  are anharmonic constants where  = 1.65832 and  = 0.32715.It may be mentioned here that coefficients  and  are complementary to coefficients  and  which were used to estimate the contribution of anharmonicity to the phonon frequency shift.The positive coefficients describe the increase in line-width with an increase in temperature.The factors calculated using (5) are listed in Table 2.
In this context, it may be pertinent to mention that the Klemens [25] model for the line-width assumed that the contribution to the line-width arises only from the decay of the optical phonon into two acoustic phonons of the same frequency and opposite momentum [3].Balkanski et al. [11] calculated the change in line-width and line center for silicon using cubic and quartic anharmonicity and found good agreement between theory and experiment.However, it has also been reported that [18] the temperature dependence of the line-width of crystalline CeO 2 can be explained completely with consideration of only the cubic anharmonic term of the crystal potential energy and is governed by the phonon occupation number.In our case, although the cubic anharmonic term is predominant, the contribution from four-phonon process is not negligible.In addition, changes in the line-width of the phonon in Raman and infrared spectra can be used as an indirect measurement for estimating the lifetime of strongly interacting optical phonons.The line-width of a Raman mode is inversely proportional to its lifetime, and anharmonic line-widths (lifetimes) increase (decrease) with increase in temperatures.The estimated lifetime can be calculated by  = 1/Γ, where  is the velocity of light and Γ is the FWHM of the mode [23].Phonon lifetime is 2.1 ps at 80 K and 0.92 ps at 440 K for the  2 phonon of CeO 2 .The variation of the decay time of the  2 phonon with increase in temperature has thus been estimated and is shown in Figure 6.The discontinuity in slope observed at about 380 K is again detectable in the figure.

Conclusions
The behavior of the optical  2 phonon of nanocrystalline cubic CeO 2 as a function of temperature is investigated via Raman spectroscopy.The  2 optical mode of CeO 2 shows phonon softening with increasing temperature up to 440 K.This phenomenon arises from dominating cubic phononphonon interaction and line-width broadening in the temperature range 80 to 440 K.The temperature dependent data analysis reveals that this mode is a combination of quasiharmonic and true anharmonic terms wherein the latter has a higher contribution to the total anharmonicity.The three-phonon process dominated the four-phonon process.The lifetime of the first-order  2 optical mode is found to decrease with increasing anharmonicity.An observable discontinuity in the slopes of the phonon softening and the line-width curves is seen at around 380 K.

Figure 1 :
Figure 1: Temperature dependent Raman spectra of CeO 2 from 80 K to 440 K.

Figure 2 :
Figure 2: Comparison of the line profile of  2 mode as a function of temperature at 80 K and 440 K.

1 )Figure 3 :
Figure 3: Temperature dependent Raman shift and FWHM of the  2 mode of CeO 2 .

− 1 )Figure 4 :
Figure 4: Frequency shift of  2 mode as a function of temperature.The solid line represents the calculated frequency shift using (1).

FWHM (cm − 1 )Figure 5 :
Figure 5: FWHM of  2 mode as a function of temperature.The solid line represents the calculated FWHM using (5).

Table 1 :
Best-fit values for anharmonic constants and true and quasiharmonic contributions in  2 mode of CeO 2 .

Table 2 :
Best-fit values for anharmonic constants for  2 mode of CeO 2 .