Anomalies of Brillouin Light Scattering in Selected Perovskite Relaxor Ferroelectric Crystals

Compositionally disordered perovskite compounds have been one of the exotic topics of research during the past several years. Colossal piezoelectric and electrostrictive effects have been observed in disordered perovskite ferroelectric materials. The key ingredient in the physical behavior of disordered perovskites is the nucleation and growth of the local dipolar regions called polar nanoregions (PNRs). PNRs begin to nucleate far above the temperature of the dielectric maximum Tm and exhibit varied relaxation behavior with temperature. The evidence for the existence of various stages in the relaxation dynamics of PNRs was revealed through the study of the temperature evolution of optical phonons by Raman scattering. The quasi-static regime of PNRs is characterized by the strong coupling between the local polarization and strain with the local structural phase transition and the critical slowing of the relaxation time. Strong anomalies in the frequency and the width of the acoustic phonons, and emergence of the central peak in the quasi-static region of the relaxation dynamics of PNRs have been observed through Brillouin scattering studies. In this review, we discuss the anomalies observed in Brillouin scattering in selected disordered perovskite ferroelectrics crystals such as Pb(Mg1/3Ta2/3)O3, Pb(Sc1/2Ta1/2)O3, 0.65PIN-0.35PT and Sr0.97Ca0.03TiO3 to understand dynamical behavior of PNRs.


Introduction
The ABO 3 perovskite-type ferroelectric oxides have been extensively studied over several decades both from the fundamental and applications point of view. The prototype ferroelectrics such as BaTiO 3 and KNbO 3 undergo a sequential structural phase transition from initial cubic phase to final rhombohedral phase, involving softening of zone-center optical phonon modes in each phase. On the other hand, PbTiO 3 exhibits a single-phase transition from cubic to tetragonal ferroelectric phase and is regarded as a text book example of the soft-mode-driven phase transition. A ferroelastic transition at T tr = 110 K driven by the softening of R 25 zone boundary phonon mode was observed in SrTiO 3 , while KTaO 3 remains paraelectric down to the lowest temperatures. In both SrTiO 3 and KTaO 3 , the low temperature dielectric constant increases upon cooling but becomes temperature independent (<10 K for SrTiO 3 and KTaO 3 ) due to quantum fluctuations. In both these systems, the quantum fluctuations suppress the ferroelectric instability and they enter a phase called quantum paraelectrics [1,2].
The substitution of ions of different valences and polarizabilities at both A and B sites in the parent perovskite compound results in the compositional disorder. This structural disorder gives rise to the formation of local dipolar impurities and defects that profoundly affect the static and dynamic properties of these materials [2]. The off centering of the A and/or B ion in the high polarizable perovskite host results in the development of net accompanied by the change in structure from cubic to rhombohedral at the temperature T c below T m [35][36][37]. The relaxor behavior becomes suppressed and the transition temperature decreases with increase in the degree of B-site ordering [37]. The nuclear magnetic resonance (NMR) studies revealed a substantial difference in the polarization dynamics between the disordered and the ordered regions of Pb(B 1/2 B" 1/2 )O 3 relaxor group [39,40]. NMR studies also revealed that the long-range polar order occurs only in fully B-site ordered region of the crystal [39,40].
In RFEs, a significant change in the dynamical behaviour of PNRs takes place in the quasi-static regime between T* and T f . This temperature range gives rise to the formation of static PNRs due to the strongly correlated shifts of A and/or B ions and associated quenched random dipolar fields.The local dipolar field results in the variation in the local strain field and both become coupled through quadratic electrostrictive effect. Brillouin light scattering spectroscopy is a valuable technique to study the coupling between the strain and polarization and relaxation dynamics of PNRs in relaxor ferroelectrics. The characteristic frequency range associated with the relaxation dynamics of PNRs is rightly located in the low frequency part of the Brillouin scattering, which is normally inaccessible through conventional Raman and thermal neutron scattering. The anomalous behavior of elastic constants in the quasi-static region in PMN and PMT was studied through Brillouin light scattering spectroscopy. Below T*, a broad central peak associated with the polarization fluctuations of PNRs was observed in PMN and PMT [41]. Both in PMN and PMT, the relaxation time of PNRs showed a characteristic slowing behavior with a marked anomaly below T f . The anomalous behavior of acoustic phonons and the central peak has been studied extensively in various Pb-based relaxor systems by Brillouin light scattering spectroscopy [42][43][44][45][46][47][48][49][50][51][52][53][54]. A significant softening of LA and TA phonon modes accompanied by a broad central peak was observed below T* in PZN-0.045PT and PZN-0.09PT [44]. In the Brillouin light scattering studies of PZN-xPT with x = 0.07, 0.1 and 0.12, a minimum was observed in the temperature derivative of the width of the LA phonon, dΓ/dT at T* [45]. Brillouin light scattering from acoustic phonons was investigated for a wide range of compositions in PMN-xPT solid solution [47]. As in the case of PZN-xPT, below T*, a significant softening of frequency of the LA phonon followed by increase in the width was observed. The activation energy of flipping of PNRs calculated using a modified superparaelectric model showed a dramatic increase across the composition of MPB. In 0.65PIN-0.35PT, anomalous behavior of the LA phonon mode at the cubictetragonal structural phase transition and two-component central peaks were studied using Brillouin light scattering spectroscopy [50]. The width of the LA phonon showed Landau-Khalatnikov-type behavior in the vicinity of the structural phase transition in a highly ordered PSN [52]. This review focuses on the Brillouin light scattering studies of dynamical aspects of PNRs in some selected RFEs carried out by the authors in the past. In the following, the results obtained in the Brillouin light scattering studies on RFEs such as PMT, PST, 0.65PIN-0.35PT and SCT-0.03 are reviewed in detail.

Brillouin Scattering
In Brillouin scattering, the incident photon undergoes inelastic scattering by the near-Brillouin zone-center acoustic phonons (q ≈ 0). The shift in the frequency of the scattered photon ∆ν B , known as the Brillouin shift is given by [53]: where ± in the above Equation (1) represents the Stokes and anti-Stokes Brillouin scattering. λ 0 , n, V and θ correspond to the wavelength of the incident photon, refractive index, phase velocity of the acoustic phonon and the scattering angle between the incident and scattered photons, respectively. The intensity of the scattered light I(q,ω) at the given wave vector q Materials 2023, 16, 605 5 of 24 and the frequency ω is related to the imaginary part of the complex dynamic susceptibility through the fluctuation-dissipation theorem, as given by [55]:

Brillouin Scattering
In Brillouin scattering, the incident photon undergoes inelastic scattering by the near-Brillouin zone-center acoustic phonons (q ≈ 0). The shift in the frequency of the scattered photon Δν B , known as the Brillouin shift is given by [53]: where ± in the above Equation (1) represents the Stokes and anti-Stokes Brillouin scattering. λ 0 , n, V and θ correspond to the wavelength of the incident photon, refractive index, phase velocity of the acoustic phonon and the scattering angle between the incident and scattered photons, respectively. The intensity of the scattered light I(q,ω) at the given wave vector q and the frequency ω is related to the imaginary part of the complex dynamic susceptibility through the fluctuation-dissipation theorem, as given by [55]: where n(ω), the Bose-Einstein thermal population factor, is given by: where is ħ the Planck's constant, kB the Boltzmann constant and T the temperature. For zone-center acoustic phonons, the dynamic susceptibility χ * (q,ω) is given by the Lorentzian response function: where f 0 is the oscillator strength, and ω o and γ o are the frequency and damping factor of the acoustic phonon. The high-contrast multi-pass Fabry-Perot interferometer (FPI) developed by Sandercock [56] is well suited to study Brillouin scattering in solids. A 3 + 3 pass tandem FPI provides adequate contrast and resolution to measure the frequency and the width of the acoustic phonons in solids in various scattering geometries. Since 3 + 3 pass FPI covers a large frequency range, 0.03-30 cm −1 , not only the acoustic phonons but also the features of the central peak can be studied.

Materials and Methods
Single crystals of PMT, 0.65PIN-0.35PT and PST were grown by the flux method as described elsewhere [50,57,58]. Single crystals of SCT-0.03 were grown by flame fusion technique [59]. The degree of B-site ordering, S parameter, for PST was determined using powder X-ray diffraction. Brillouin light scattering measurements for all the crystals were performed using a high-contrast 3+3 pass Sandercock tandem Fabry-Perot Interferometer in backscattering geometry. Brillouin spectra were excited using a single longitudinal mode diode-pumped solid-state laser operating at 532 nm. Measurements were carried out in different free spectral ranges (FSR) in the temperature range 870-100 K using a Linkam (THMS 600) stage. The frequency and the width of acoustic phonons were deduced by fitting the experimental spectra with the Voigt function. During fitting, the Gaussian part of the Voigt function was kept constant to account for the instrumental broadening. The instrumental broadening was obtained by fitting the transmitted peak of the interferometer with the Gaussian function. The full width at half maximum of the π [n(ω) + 1] Imχ * (q, ω ) (2) where n(ω), the Bose-Einstein thermal population factor, is given by:

Brillouin Scattering
In Brillouin scattering, the incident photon undergoes inelastic scattering by the near-Brillouin zone-center acoustic phonons (q ≈ 0). The shift in the frequency of the scattered photon Δν B , known as the Brillouin shift is given by [53]: where ± in the above Equation (1) represents the Stokes and anti-Stokes Brillouin scattering.
λ 0 , n, V and θ correspond to the wavelength of the incident photon, refractive index, phase velocity of the acoustic phonon and the scattering angle between the incident and scattered photons, respectively. The intensity of the scattered light I(q,ω) at the given wave vector q and the frequency ω is related to the imaginary part of the complex dynamic susceptibility through the fluctuation-dissipation theorem, as given by [55]: where n(ω), the Bose-Einstein thermal population factor, is given by: where is ħ the Planck's constant, kB the Boltzmann constant and T the temperature. For zone-center acoustic phonons, the dynamic susceptibility χ * (q,ω) is given by the Lorentzian response function: where f 0 is the oscillator strength, and ω o and γ o are the frequency and damping factor of the acoustic phonon.
The high-contrast multi-pass Fabry-Perot interferometer (FPI) developed by Sandercock [56] is well suited to study Brillouin scattering in solids. A 3 + 3 pass tandem FPI provides adequate contrast and resolution to measure the frequency and the width of the acoustic phonons in solids in various scattering geometries. Since 3 + 3 pass FPI covers a large frequency range, 0.03-30 cm −1 , not only the acoustic phonons but also the features of the central peak can be studied.

Materials and Methods
Single crystals of PMT, 0.65PIN-0.35PT and PST were grown by the flux method as described elsewhere [50,57,58]. Single crystals of SCT-0.03 were grown by flame fusion technique [59]. The degree of B-site ordering, S parameter, for PST was determined using powder X-ray diffraction. Brillouin light scattering measurements for all the crystals were performed using a high-contrast 3+3 pass Sandercock tandem Fabry-Perot Interferometer in backscattering geometry. Brillouin spectra were excited using a single longitudinal mode diode-pumped solid-state laser operating at 532 nm. Measurements were carried out in different free spectral ranges (FSR) in the temperature range 870-100 K using a Linkam (THMS 600) stage. The frequency and the width of acoustic phonons were deduced by fitting the experimental spectra with the Voigt function. During fitting, the Gaussian part of the Voigt function was kept constant to account for the instrumental broadening. The instrumental broadening was obtained by fitting the transmitted peak of the interferometer with the Gaussian function. The full width at half maximum of the

Brillouin Scattering
In Brillouin scattering, the incident photon undergoes inelastic scattering by the near-Brillouin zone-center acoustic phonons (q ≈ 0). The shift in the frequency of the scattered photon Δν B , known as the Brillouin shift is given by [53]: where ± in the above Equation (1) represents the Stokes and anti-Stokes Brillouin scattering.
λ 0 , n, V and θ correspond to the wavelength of the incident photon, refractive index, phase velocity of the acoustic phonon and the scattering angle between the incident and scattered photons, respectively. The intensity of the scattered light I(q,ω) at the given wave vector q and the frequency ω is related to the imaginary part of the complex dynamic susceptibility through the fluctuation-dissipation theorem, as given by [55]: where n(ω), the Bose-Einstein thermal population factor, is given by: where is ħ the Planck's constant, kB the Boltzmann constant and T the temperature. For zone-center acoustic phonons, the dynamic susceptibility χ * (q,ω) is given by the Lorentzian response function: where f 0 is the oscillator strength, and ω o and γ o are the frequency and damping factor of the acoustic phonon.
The high-contrast multi-pass Fabry-Perot interferometer (FPI) developed by Sandercock [56] is well suited to study Brillouin scattering in solids. A 3 + 3 pass tandem FPI provides adequate contrast and resolution to measure the frequency and the width of the acoustic phonons in solids in various scattering geometries. Since 3 + 3 pass FPI covers a large frequency range, 0.03-30 cm −1 , not only the acoustic phonons but also the features of the central peak can be studied.

Materials and Methods
Single crystals of PMT, 0.65PIN-0.35PT and PST were grown by the flux method as described elsewhere [50,57,58]. Single crystals of SCT-0.03 were grown by flame fusion technique [59]. The degree of B-site ordering, S parameter, for PST was determined using powder X-ray diffraction. Brillouin light scattering measurements for all the crystals were performed using a high-contrast 3+3 pass Sandercock tandem Fabry-Perot Interferometer in backscattering geometry. Brillouin spectra were excited using a single longitudinal mode diode-pumped solid-state laser operating at 532 nm. Measurements were carried out in different free spectral ranges (FSR) in the temperature range 870-100 K using a Linkam (THMS 600) stage. The frequency and the width of acoustic phonons were deduced by fitting the experimental spectra with the Voigt function. During fitting, the Gaussian part of the Voigt function was kept constant to account for the instrumental broadening. The instrumental broadening was obtained by fitting the transmitted peak of the interferometer with the Gaussian function. The full width at half maximum of the is the Planck's constant, k B the Boltzmann constant and T the temperature. For zone-center acoustic phonons, the dynamic susceptibility χ * (q,ω) is given by the Lorentzian response function: where f 0 is the oscillator strength, and ω o and γ o are the frequency and damping factor of the acoustic phonon. The high-contrast multi-pass Fabry-Perot interferometer (FPI) developed by Sandercock [56] is well suited to study Brillouin scattering in solids. A 3 + 3 pass tandem FPI provides adequate contrast and resolution to measure the frequency and the width of the acoustic phonons in solids in various scattering geometries. Since 3 + 3 pass FPI covers a large frequency range, 0.03-30 cm −1 , not only the acoustic phonons but also the features of the central peak can be studied.

Materials and Methods
Single crystals of PMT, 0.65PIN-0.35PT and PST were grown by the flux method as described elsewhere [50,57,58]. Single crystals of SCT-0.03 were grown by flame fusion technique [59]. The degree of B-site ordering, S parameter, for PST was determined using powder X-ray diffraction. Brillouin light scattering measurements for all the crystals were performed using a high-contrast 3 + 3 pass Sandercock tandem Fabry-Perot Interferometer in backscattering geometry. Brillouin spectra were excited using a single longitudinal mode diode-pumped solid-state laser operating at 532 nm. Measurements were carried out in different free spectral ranges (FSR) in the temperature range 870-100 K using a Linkam (THMS 600) stage. The frequency and the width of acoustic phonons were deduced by fitting the experimental spectra with the Voigt function. During fitting, the Gaussian part of the Voigt function was kept constant to account for the instrumental broadening. The instrumental broadening was obtained by fitting the transmitted peak of the interferometer with the Gaussian function. The full width at half maximum of the Gaussian fit is the measure of the instrumental broadening. The central peak was fitted using a Lorentzian function with ω o = 0.

Elastic Anomalies and Polarization Dynamics in PMT
The Brillouin spectra of PMT measured in the FSR range of 75 GHz at different temperatures are shown in Figure 1. A single LA phonon mode with a frequency around 42 GHz can be seen. Since the macroscopic structure of PMT is cubic, the selection rule permits only the excitation of LA phonon mode in the backscattering geometry. The corresponding elastic constant C 11 is given by the following expressions [60]: where ν LA , Γ, V 11 and α in Equation (5) are the frequency, width, velocity and attenuation of the LA phonon, and ρ and n are the density and the refractive index of PMT and λ o is the wavelength of the incident photon. Using the parameters of the LA phonon, C 11 and α were calculated using the known values of ρ = 9.65 gm/cm 3 and the refractive index n = 2.475 for PMT [60]. The temperature dependence of the refractive index was not considered in the calculation of C 11 , as it does not significantly affect its magnitude. Figure 2 shows the temperature dependence of C 11 and α for PMT. Gaussian fit is the measure of the instrumental broadening. The central peak was fitted using a Lorentzian function with  o = 0.

4.1.Elastic Anomalies and Polarization Dynamics in PMT
The Brillouin spectra of PMT measured in the FSR range of 75 GHz at different temperatures are shown in Figure 1. A single LA phonon mode with a frequency around 42 GHz can be seen. Since the macroscopic structure of PMT is cubic, the selection rule permits only the excitation of LA phonon mode in the backscattering geometry. The corresponding elastic constant C 11 is given by the following expressions [60]: where  LA , , V 11 and α in Equation (5) are the frequency, width, velocity and attenuation of the LA phonon, and  and n are the density and the refractive index of PMT and  o is the wavelength of the incident photon. Using the parameters of the LA phonon, C 11 and α were calculated using the known values of  = 9.65 gm/cm 3 and the refractive index n = 2.475 for PMT [60]. The temperature dependence of the refractive index was not considered in the calculation of C 11 , as it does not significantly affect its magnitude. Figure 2 shows the temperature dependence of C 11 and α for PMT. The elastic constant C 11 shows a broad minimum and α shows a broad maximum in the region of the dielectric anomaly at T m . Our data are in good agreement with those of the earlier studies on the elastic anomalies observed in PMT [60]. The minima of C 11 agree well with the corresponding dielectric anomaly of the PMT [42]. In the present review, a mean field model that was previously used to study the coupling between the polarization and the strain in KTN is employed [9]. The temperature dependence of C 11 according The elastic constant C 11 shows a broad minimum and α shows a broad maximum in the region of the dielectric anomaly at T m . Our data are in good agreement with those of the earlier studies on the elastic anomalies observed in PMT [60]. The minima of C 11 agree well with the corresponding dielectric anomaly of the PMT [42]. In the present review, a mean field model that was previously used to study the coupling between the polarization and the strain in KTN is employed [9]. The temperature dependence of C 11 according to this model is given by the following expression: The first two terms of Equation (6) describe the normal anharmonic process, and the third term accounts for the critical softening of C 11 due to the critical fluctuations. For three-dimensional (3D) fluctuations, the critical exponent p = 0.5, for 2D and 1D fluctuations; it is 1 and 1.5 [61]. Since PNRs begin to appear around T B , the temperature dependence of C 11 was fitted with this equation for temperatures higher than T m . As the structure of PMT is cubic, it is reasonable to assume p to be 0.5. C 11 was therefore fitted with the critical exponent p fixed at 0.5 and all the other parameters in Equation (6) are allowed to vary freely [9]. The parameters obtained from the fitting are: A 1 = 223.45 GPa, A 2 = 0.024 GPa/m, A 3 = 36.5 GPa and T o = 463 K. Figure 3 shows the result of the fitting for PMT.  The experimental C 11 shows good agreement with the mean-field perature above T o + 40K. C 11 shows a marked deviation from the meanthis temperature and exhibits reduced softening. The attenuation  a precisely at the temperature where C 11 deviates from the mean field m ature at which C 11 begins to show reduced softening can be identified teristic percolation temperature T * [26]. T * obtained for PMT is 503 K agreement with that of the X-ray diffraction data [30]. The temperature between T B and T * is determined by the third term of Equation (6), whic quasi-dynamic fluctuations inpolarization. The experimental C 11 shows good agreement with the mean-field model in the temperature above T o + 40 K. C 11 shows a marked deviation from the mean-field model below this temperature and exhibits reduced softening. The attenuation α also begins to rise precisely at the temperature where C 11 deviates from the mean field model. The temperature at which C 11 begins to show reduced softening can be identified with the characteristic percolation temperature T* [26]. T* obtained for PMT is 503 K, which is well in agreement with that of the X-ray diffraction data [30]. The temperature dependence of C 11 between T B and T* is determined by the third term of Equation (6), which accounts for the quasi-dynamic fluctuations inpolarization.
The reduced softening of C 11 and increased α below T * indicates the appearance of excess strain in the crystal. This additional strain ∆S arises from the electrostrictive coupling between the polarization and strain and is expressed as ∆S = (γ 11 + 2γ 12 ) P 2 , where γ 11 and γ 12 are the electrostriction coefficients [2]. Equation (6) does not account for the excess contribution from the electrostriction γ to the elastic constant C 11 . The reduced softening of C 11 below T* can be qualitatively explained using Kubo's formalism [62]. According to this formalism, the polarization fluctuations δP contain two components with very distinct dynamics: Here, δP d is the quasi-dynamic part, i.e., the third term of Equation (6), and δP s corresponds to much slower fluctuations (quasi-static) with the development of long-range correlations between the PNRs. The quasi-static part δP s is expressed in terms of the Edwards-Anderson order parameter q EA = V −1 δP s 2 dV, where V is the volume over which long-range correlations exist [62,63]. The change in the elastic constant ∆C in the quasi-static region can be written as [63]: where χ and χ are the real and imaginary parts of the low-frequency dielectric susceptibility. In the disordered Rb 1−x (NH 4 ) x H 2 PO 4 , the quasi-static term of Equation (8) makes a dominant contribution to the anomalous region of the frequency and width of LA phonon mode [62,63]. A detailed study of the contribution ofthe quasi-static term ∆C s on the temperature dependence of C 11 for Pb-based relaxors is yet to be carried out. However, a qualitative idea of the temperature dependence of ∆C s can be obtained from the study of q EA . The direct determination of q EA by 93 Nb NMR for PMN shows a strong increase with a decrease in temperature below T* [64]. The deviation of the elastic constant C 11 from the mean field behavior below T* therefore evidently arises from the additional term ∆C s of the quasi-static polarization within PNRs. In fact, in PMN, the reduced softening of C 11 almost coincides with the rise inthe local polarization P d below T* [65]. Therefore, in PMT, the reduced softening of C 11 and the increased attenuation in the quasi-static region are evidently associated with the appearance of additional strain due to the strong coupling between the polarization and the strain. In PMT, a broad central peak with marked temperature dependence appears below T*. The typical central peak spectra measured in the FSR range of 500 GHz at different temperatures are displayed in Figure 4a-d. The intensity of the central peak can be described as a Debye-type relaxation according to the following equation [66,67]: where τ is the relaxation time.In the strong damping limit, (γ o >>ω o ), the Lorentzian function given by Equation (4) reduces to the Debye-type function with ω o = 0 [55]. The width of the central peak obtained from the fit of the Lorentzian function is then proportional to 1/τ. The width of the central peak obtained from the fit is shown in Figure 5. With the decrease in temperature below T*, the width decreases with a change in slope between T m and 160 K. An anomalous decrease in the width of the central peak can be observed below 160 K. This temperature corresponds to the freezing temperature T f for PMT [68]. The relaxation time of the central peak corresponds to the short relaxation time τ s obtained from the broadband dielectric measurements [68]. and 160 K. An anomalous decrease in the width of the central peak can be low 160 K. This temperature corresponds to the freezing temperature T f for P relaxation time of the central peak corresponds to the short relaxation tim from the broadband dielectric measurements [68].  To obtain a qualitative understanding of the relaxation behavior, the temperature dependence of relaxation time was analyzed using the simple Arrhenius law  =  o exp[E a /kT], as shown in Figure 6. The parameters obtained from the fit are  o = 2.0 ± 0.14× 10 −12 S, E a = 220 ± 9 K for PMT. Though the fitting parameters are reasonable, the To obtain a qualitative understanding of the relaxation behavior, the temperature dependence of relaxation time was analyzed using the simple Arrhenius law τ = τ o exp[E a /kT], as shown in Figure 6. The parameters obtained from the fit are τ o = 2.0 ± 0.14 × 10 −12 S, E a = 220 ± 9 K for PMT. Though the fitting parameters are reasonable, the Arrhenius behavior does not agree well with the experimental data. The observed non-Arrhenius behavior suggests that the relaxation dynamics are not governed by dipolar fluctuations alone. Broadband dielectric studies in Pb-based relaxor systems suggest that the relaxation dynamics of PNRs are not governed by a single relaxation mechanism [68][69][70]. According to these studies, the dominant mechanism of relaxation between T* and T f is the dipole flipping, while that below T f is the volume fluctuations (breathing) of PNRs [68][69][70]. The strong increase in the relaxation time below T f suggests that the relaxation mechanism above T f is different from that below T f . Therefore, the observed critical slowing of the relaxation time of the central peak between T* and T f can be attributed to the dipole flipping and below T f to the volume fluctuations of PNRs. This interpretation seems to be quite consistent with the existence of different stages of relaxation dynamics of PNRs as revealed by Raman scattering studies [26]. The increase in correlations between PNRs in the quasi-static regime between T * and T f increase the relaxation time. In the temperature range below T f , due to the freezing, the dipolar flipping in PNRs is no longer possible. Therefore, the observed anomalous increase in the relaxation time below T f could be due to the volume fluctuations of PNRs. Figure 7 shows the powder X-ray diffraction pattern of PST. The broad superlattice reflection due to the B-site ordering is observed at 2 = 19 o . This interpretation seems to be quite consistent with the existence of different stages of relaxation dynamics of PNRs as revealed by Raman scattering studies [26]. The increase in correlations between PNRs in the quasi-static regime between T* and T f increase the relaxation time. In the temperature range below T f , due to the freezing, the dipolar flipping in PNRs is no longer possible. Therefore, the observed anomalous increase in the relaxation time below T f could be due to the volume fluctuations of PNRs. Figure 7 shows the powder X-ray diffraction pattern of PST. The broad superlattice reflection due to the B-site ordering is observed at 2θ = 19 o . crease in the relaxation time below T f could be due to the volume fluctuations of PNR Figure 7 shows the powder X-ray diffraction pattern of PST. The broad superlat reflection due to the B-site ordering is observed at 2 = 19 o .  is the theoretical intensity ratio for the completely ordered PST. The theoretical intensity ratio I 111 /I 200 for the completely ordered PST is estimated to be 1.33 [71]. The S parameter for the as-grown PST single crystals is found to be S = 0.55. This indicates that the as-grown PST single crystal is partially disordered. Figure 8 shows the low-frequency dielectric dispersion of partially ordered PST. Since it is known that an increase in B-site ordering suppresses the relaxor behavior, the partial ordering in PST results in the weak relaxor-like behavior with T m = 297 K. is the theoretical intensity ratio for the completely ordered PST. The theoretical intensity ratio I 111 /I 200 for the completely ordered PST is estimated to be 1.33 [71]. The S parameter for the as-grown PST single crystals is found to be S = 0. 55. This indicates that the as-grown PST single crystal is partially disordered. Figure 8 shows the low-frequency dielectric dispersion of partially ordered PST. Since it is known that an increase in B-site ordering suppresses the relaxor behavior, the partial ordering in PST results in the weak relaxor-like behavior with T m = 297 K. The typical Brillouin spectra of PST measured in the FSR range of 75 GHz at different temperatures are displayed in the Figure 9. In the cubic phase above 295 K, only the LA phonon mode could be observed at the backscattering geometry. Since PST exhibits a spontaneous phase transition from cubic to rhombohedral structure [72], a transverse acoustic (TA) phonon mode at 27 GHz, allowed by symmetry, appears below 295 K in the rhombohedral phase. However, the intensity of the TA phonon was too weak to obtain a reliable value for the frequency and width. Hence, the temperature dependence of TA The typical Brillouin spectra of PST measured in the FSR range of 75 GHz at different temperatures are displayed in the Figure 9. In the cubic phase above 295 K, only the LA phonon mode could be observed at the backscattering geometry. Since PST exhibits a spontaneous phase transition from cubic to rhombohedral structure [72], a transverse acoustic (TA) phonon mode at 27 GHz, allowed by symmetry, appears below 295 K in the rhombohedral phase. However, the intensity of the TA phonon was too weak to obtain a reliable value for the frequency and width. Hence, the temperature dependence of TA mode is not discussed further.   Interaction between the acoustic phonon and polarization affects the static and dynamic parts of the complex elastic constant C. Depending on the form of the interaction energy, the interaction can be classified as either linear or nonlinear. For bilinear coupling, the interaction energy F c = SP, where S is the strain, P is the polarization and  is the piezoelectric coefficient. F c = SP 2 for electrostrictive coupling, where  is the electro-  Figure 10 displays the temperature variation in frequency ω LA and the width Γ LA of the LA phonon. Above 600 K, the frequency and the width of the LA phonon are almost independent of temperature. The LA phonon mode begins to deviate from the normal hightemperature behavior below 600 K and exhibits marked softening with a sharp minimum at T c = 295 K, the transition temperature. The width of the phonon mode begins to show an anomalous increase at 450 K with a sharp maximum at T c .    Interaction between the acoustic phonon and polarization affects the static and dynamic parts of the complex elastic constant C. Depending on the form of the interaction energy, the interaction can be classified as either linear or nonlinear. For bilinear coupling, the interaction energy F c = βSP, where S is the strain, P is the polarization and β is the piezoelectric coefficient. F c = γSP 2 for electrostrictive coupling, where γ is the electrostric-tive coefficient. The coupling between the strain and polarization affects both the static and dynamic parts of the complex elastic constant C, as given by [73][74][75][76]:

Anomalous Behavior of Longitudinal Acoustic Phonons in PST
In Equation (10), C ∞ is the elastic constant far from the transition temperature and ∆C is the change in elastic constant due to the coupling between strain and polarization. For bilinear coupling, the elastic constant softens as /T − T c / −1 on both sides of the transition. For electrostrictive coupling, the elastic constant abruptly decreases to a lower value at T c [73,76]. Brillouin scattering measures the dynamic part of C according to Equation (11): where ρ is the bulk density of the crystal and V is the corresponding velocity of the acoustic phonon. Assuming a Debye-type relaxation, the dynamic behavior of C is given by the following form [74,75]: where τ is the relaxation time. The frequency of the acoustic phonon is related to the real part of the elastic constant, while the width is related to the imaginary part of the above equation, as given by: The relaxation mechanism is also known as Landau-Khalatnikov (LK) mechanism. The line shape of the width of the LA phonon, as shown in Figure 10, is markedly asymmetric about T c , as expected from the LK mechanism.
Using Equations (5), (13) and (14), the relaxation time can be expressed as: where υ(Γ) is the frequency (width) of the acoustic phonon at a given temperature and υ ∞ (Γ ∞ ) is the corresponding value of the acoustic phonon far above the transition temperature. In the relaxor-like system Rb 1−x (NH 4 ) x H 2 PO 4 with x= 0.35, the temperature dependence of the width of the LA phonon in the quasi-static region was shown to be equivalent toan LK-type mechanism [62]. Similarly, in the relaxor phase of the PST, the anomalous increase in the width of the LA phonon below T*is evidently connected with the development of the quasi-static part δP s . Therefore, a linear coupling between the polarization and strain should exist in the quasi-static region below T* [62]. Moreover, according to Young and Scott, Equation (12) applies to any form of coupling, such as bilinear or electrostrictive coupling and the coupling between the acoustic phonon mode and other excitations, which are altogether not connected with the order parameter [76]. This relaxation mechanism was used to describe the relaxation behavior of the acoustic phonon mode across commensurate-incommensurate phase transition in Ba 2 NaNb 5 O 15 where the exact form of the coupling was not known [76,77]. Based on the above studies in Rb 1−x (NH 4 ) x H 2 PO 4 and Ba 2 NaNb 5 O 15 , we calculated the relaxation time both in the ferroelectric phase below T c and in the quasi-static region below T* (in the temperature range between T* and T c ) using the LK relaxation mechanism. Figure 11 shows the temperature variation inthe relaxation time τ LA of the LA phonon. range between T * and T c ) using the LK relaxation mechanism. Figure 11 shows the temperature variation inthe relaxation time  LA of the LA phonon. The relaxation time exhibits a critical slowing on approachingT c from the ferroelectric and the relaxor phase. As shown in Figure 11, the temperature variation inthe relaxation rate exhibits approximately linear behavior where 1/ LA satisfies the following relation: for typical order-disorder ferroelectrics such as potassium dihydrogen phosphate KH 2 PO 4 (KDP) and triglycine sulfate (NH 2 CH 2 COOH) 3 H 2 SO 4 (TGS); o is 0.13 ps for KDP [78] and 0.10 ps for TGS [79]. Since the structural phase transition in PST does not involve the condensation of soft mode, the transition can be of the order-disorder-type. The local structure of PNRs in Pb-based relaxors is known to be rhombohedral [7] and sufficiently strong interaction between the quasi-static PNRs in the temperature range below T * may eventually lead to the formation of long-range ferroelectric order with the macroscopic rhombohedral structure below T c .

Brillouin Spectroscopy of Cubic-Tetragonal Structural Phase Transition in a 0.65PIN-0.35PT Single Crystal
PIN-PT solid solution, like PMN-PT and PZN-PT, shows the sequence of structural transition with temperature and PT content and has MPB for x = 0.37 [80]. PIN-PT composition near MPB shows a structural phase transition from cubic to tetragonal (T ct ) and tetragonal to rhombohedral (T tr ) phases [43,81]. While the dynamical behavior of PNRs has been widely studied in the PMN-PT and PZN-PT systems, no such studies have been The relaxation time exhibits a critical slowing on approaching T c from the ferroelectric and the relaxor phase. As shown in Figure 11, the temperature variation inthe relaxation rate exhibits approximately linear behavior where 1/τ LA satisfies the following relation: The solid line is the linear fit to the experimental data. Best fit of the experimental data with Equation (16) yields τ o = 0.22 ps in the ferroelectric phase and 0.19 ps in the relaxor phase. The relaxation time τ o for PST is of a similar order of magnitude as observed for typical order-disorder ferroelectrics such as potassium dihydrogen phosphate KH 2 PO 4 (KDP) and triglycine sulfate (NH 2 CH 2 COOH) 3 H 2 SO 4 (TGS); τ o is 0.13 ps for KDP [78] and 0.10 ps for TGS [79]. Since the structural phase transition in PST does not involve the condensation of soft mode, the transition can be of the order-disorder-type. The local structure of PNRs in Pb-based relaxors is known to be rhombohedral [7] and sufficiently strong interaction between the quasi-static PNRs in the temperature range below T* may eventually lead to the formation of long-range ferroelectric order with the macroscopic rhombohedral structure below T c .

Brillouin Spectroscopy of Cubic-Tetragonal Structural Phase Transition in a 0.65PIN-0.35PT Single Crystal
PIN-PT solid solution, like PMN-PT and PZN-PT, shows the sequence of structural transition with temperature and PT content and has MPB for x = 0.37 [80]. PIN-PT composition near MPB shows a structural phase transition from cubic to tetragonal (T ct ) and tetragonal to rhombohedral (T tr ) phases [43,81]. While the dynamical behavior of PNRs has been widely studied in the PMN-PT and PZN-PT systems, no such studies have been carried out for PIN-PT solid solutions. In this review, the Brillouin scattering studies of the relaxation dynamics of PNRs in 0.65PIN-0.35PT single crystals are discussed.
The Brillouin spectra measured at various temperatures in the FSR range of 75 GHz are shown in Figure 12. The Brillouin doublet of the LA phonon can be seen clearly. A faint and broad TA phonon mode can also be observed in the cubic phase. This could be due to the presence of tetragonal-structure PNRs. The TA mode, however, could not be followed for all the temperatures and hence is not discussed further. Figure 13 shows the frequency and the width Γ LA of the LA phonon as a function of temperature. Above 700 K, the minor temperature dependence in the frequency and the width of the LA phonon is attributed to the normal anharmonic process. Below 700 K, the frequency of the LA phonon shows substantial softening, followed by a strong increase in its width on approaching cubic to tetragonal phase transition temperature T ct = 540 K. Also clear are the LA anomalies that can be seen at the tetragonal-rhombohedral transition temperature T tr = 460 K.
due to the presence of tetragonal-structure PNRs. The TA mode, however, could followed for all the temperatures and hence is not discussed further. Figure 13 sh frequency and the width  LA of the LA phonon as a function of temperature. Abov the minor temperature dependence in the frequency and the width of the LA ph attributed to the normal anharmonic process. Below 700 K, the frequency of phonon shows substantial softening, followed by a strong increase in its width proaching cubic to tetragonal phase transition temperature T ct = 540 K. Also clear LA anomalies that can be seen at the tetragonal-rhombohedral transition tempera = 460 K.  The observed anomalous behavior in the frequency and the width of the LA below 700 K clearly suggests that T B~ 700 K. As discussed in the Introduction, the electrostrictive coupling between the local polarization and strain in PNRs lead anomalous behavior of the LA phonon below T B . As can be seen from Figure 12, a CP (NCP) begins to appear at about 620 K and its intensity grows further subst upon cooling. Figure 14a-d show the Brillouin spectra measured in the FSR rang faint and broad TA phonon mode can also be observed in the cubic phase. This could be due to the presence of tetragonal-structure PNRs. The TA mode, however, could not be followed for all the temperatures and hence is not discussed further. Figure 13 shows the frequency and the width  LA of the LA phonon as a function of temperature. Above 700 K, the minor temperature dependence in the frequency and the width of the LA phonon is attributed to the normal anharmonic process. Below 700 K, the frequency of the LA phonon shows substantial softening, followed by a strong increase in its width on approaching cubic to tetragonal phase transition temperature T ct = 540 K. Also clear are the LA anomalies that can be seen at the tetragonal-rhombohedral transition temperature T tr = 460 K.  The observed anomalous behavior in the frequency and the width of the LA phonon below 700 K clearly suggests that T B~ 700 K. As discussed in the Introduction, the strong electrostrictive coupling between the local polarization and strain in PNRs leads to the anomalous behavior of the LA phonon below T B . As can be seen from Figure 12, a narrow CP (NCP) begins to appear at about 620 K and its intensity grows further substantially upon cooling. Figure 14a-d show the Brillouin spectra measured in the FSR range of 400 The observed anomalous behavior in the frequency and the width of the LA phonon below 700 K clearly suggests that T B~7 00 K. As discussed in the Introduction, the strong electrostrictive coupling between the local polarization and strain in PNRs leads to the anomalous behavior of the LA phonon below T B . As can be seen from Figure 12, a narrow CP (NCP) begins to appear at about 620 K and its intensity grows further substantially upon cooling. Figure 14a-d show the Brillouin spectra measured in the FSR range of 400 GHz. A broad CP begins (BCP) to appear below T B and its intensity increases with a further decrease in temperature. The width of both NCP and BCP decreases upon cooling toward T ct , reflecting the growth of the PNRs and the corresponding slowing of the relaxation time. The relaxation time of the LA phonon was calculated using Equation (15). We use the highest temperature value of the frequency and width of the LA mode for ν ∞ and Γ ∞ . In Figure 15, the calculated value of the relaxation time from the above expression is shown along with the relaxation times obtained for narrow and broad central peaks above T ct . toward T ct , reflecting the growth of the PNRs and the corresponding slowing of the relaxation time. The relaxation time of the LA phonon was calculated using Equation (15).
We use the highest temperature value of the frequency and width of the LA mode for   and   . In Figure 15, the calculated value of the relaxation time from the above expression is shown along with the relaxation times obtained for narrow and broad central peaks above T ct .   sion is shown along with the relaxation times obtained for narrow and broad central peaks above T ct .   The relaxation time for the LA phonon in the ferroelectric phase was not calculated due to the lack of sufficient data points. Below 600 K, the relaxation time of the local strain, τ LA qualitatively agrees with the polarization fluctuations (as determined by BCP) of PNRs. Though the broad central peak appears below T B , it effectively becomes coupled with the LA phonon only from 600 K. This result implies that the relaxation process of PNRs below 600 K is determined by the polarization fluctuations coupled with the local strain fluctuations. A similar result was obtained for 0.71PNN-0.29PT [82]. In this solid solution system, broad and narrow central peaks were observed around T* = 500 K. The relaxation time of the LA phonon agrees well with that of the broad central peak. Considering the similarity of relaxation dynamics with 0.71PNN-0.29PT, we suggest that for 0.65PIN-0.35PT, the temperature 600 K at which the narrow central peak appears can be taken as T*. At this point, we would like to remark on the two-component (narrow and broad) central peaks observed in 0.65PIN-0.35PT. Two-component central peaks appear to be the common feature in all the PT-substituted relaxor solid solutions [44,74,82]. We suggest the possible origin of two-component central peaks in the relaxor-PT solid solution. Owing to the high polarizability of the Ti 4+ ion, the substitution of PT may lead to the formation of PNRs with a tetragonal structure in addition to the already existing PNRs with rhombohedral structure. With decrease in temperature, the size of tetragonal PNRs may grow faster and eventually become larger than the size of the rhombohedral one. In this case, the activation energy required for the flipping of dipoles along equivalent directions in tetragonal PNRs becomes higher. The higher activation energy results in a longer relaxation time for tetragonal PNRs. Therefore, we attribute the narrow component of the central peak to the relaxation of tetragonal PNRs and the broad component to the rhombohedral PNRs. It may be noted that in both PST and PSN, increased correlations between rhombohedral PNRs anddecrease in temperature results in long-range ferroelectric order with the macroscopic rhombohedral structure below T c . In both PSN and PST, only a broad central peak corresponding to rhombohedral PNRs was observed below T* [46,51]. Therefore, we suggest that in 0.65PIN-0.35PT, the phase transition from cubic relaxor phase to the tetragonal ferroelectric phase could be due to the nucleation and growth of tetragonal PNRs, whose relaxation dynamics are different from that of the rhombohedral PNRs.

Brillouin Light Scattering Studies of Acoustic Phonon Anomalies in SCT-0.03 Quantum Ferroelectrics
The Brillouin spectra of SCT-0.03 in some temperature ranges are displayed in Figure 16a,b. A longitudinal acoustic (LA) phonon peak around 75 GHz and a transverse acoustic (TA) mode at 45 GHz can be seen. As shown in Figure 16b, the TA mode splits into two components (TA1 and TA2) below 280 K (indicated by the arrow mark). The temperature variation in the frequency and the width of the LA phonon is shown in Figure 17. In the temperature range 500-323 K, a linear increase in the frequency of the LA phonon with the near constancy of its width is ascribed to the normal anharmonic process. Below 323 K, the LA phonon shows anomalous behavior with a sharp minimum in the frequency followed by a sharp maximum in width at the ferroelastic phase transition temperature T ct = 280 K [83]. Similar elastic anomalies are observed in nondoped SrTiO 3 crystal by picosecond ultrasound spectroscopy [84]. The remarkable increase in T ct by Ca doping can be caused by the smaller ionic radius of Ca ions, which induces the antiferro distortive rigid rotation of the oxygen octahedral [85]. A similar increase in T ct was also reported in La-doped SrTiO 3 [86].    The width of the LA phonon is markedly asymmetric in the vicinity of the tra temperature as expected from the LK mechanism. In the X-ray and the neutron tion studies of SCT-0.04, the cubic to tetragonal phase transition with the tra temperature of 210 K was reported [20]. The phase transition in SCT-0.04 was obse be a ferroelastic-type as it was not associated with any dielectric anomalies. The the observed anomalous behavior of the LA phonon at 280 K in SCT-0.03 in the p study is attributed to the ferroelastic-type structural phase transition from cubi tragonal phase. The higher transition temperature in SCT-0.03 could be due to the Ca content. The observed anomalies of the LA phonon in the phase transition reg The width of the LA phonon is markedly asymmetric in the vicinity of the transition temperature as expected from the LK mechanism. In the X-ray and the neutron diffraction studies of SCT-0.04, the cubic to tetragonal phase transition with the transition temperature of 210 K was reported [20]. The phase transition in SCT-0.04 was observed to be a ferroelastic-type as it was not associated with any dielectric anomalies. Therefore, the observed anomalous behavior of the LA phonon at 280 K in SCT-0.03 in the present study is attributed to the ferroelastic-type structural phase transition from cubic to tetragonal phase. The higher transition temperature in SCT-0.03 could be due to the lower Ca content. The observed anomalies of the LA phonon in the phase transition region are caused by the interaction between the phonon mode and the order parameter. It can be seen from Figure 17 that the order parameter interaction with the LA phonon does not agree with the bilinear coupling that requires /T −T ct / −1 -type temperature dependence on either side of the phase transition. Slonczewski and Thomas reported that the ferroelastic phase transition at 110 K in SrTiO 3 can be described by coupling between the strain S and the order parameter Q in the form F c = δSQ 2 , where the order parameter Q is the soft mode coordinate and δ is the coupling coefficient [87]. The ferroelastic transition in SrTiO 3 could be well described by Landau theory, taking into account the interaction between the strain and the soft optical phonon mode.We suggest that in SCT also, the observed acoustic anomaly could be due to the coupling between the LA phonon and the soft mode. The temperature dependence of the relaxation time τ LA of the LA phonon, calculated using Equation (15), is plotted in Figure 18. The relaxation time of the LA phonon exhibits critical slowing when approaching T ct from both sides of the transition. The value of τ o obtained from Equation (16) is 0.077 ps below T ct and 0.057 ps above T ct . This is about half of the value of KDP and TGS and is very similar to that obtained in PbHPO 4 [75]. This suggests the weaker coupling of the LA phonon with the soft mode of the paraelectric host.
The temperature dependence of the frequency and width of the TA1 and TA2 phonons are shown in Figures 19 and 20. The frequency and the width of the TA1 phonon mode exhibit similar behavior as that of the LA phonon mode. However, a large scattering of the data observed in the width of the TA1 phonon mode limits further detailed discussion on its temperature dependence. The TA2 mode appears just near T ct . The experimental data, however, could be unambiguously fitted only below 265 K. The TA2 mode shows a significant softening with the decrease in temperature in the ferroelastic phase. In the Brillouin scattering study of phase transition in Cs 2 CdBr 4 , the splitting of TA mode into two components was observed on the phase transition from orthorhombic to monoclinic phase [88]. One of the components was purely TA mode (C 66 ) corresponding to the parent orthorhombic phase, while the other was ascribed to the quasi-transverse mode γ 4 (b), propagating along the b axis of the crystal. The frequency of the γ 4 (b) mode was observed to decrease with decrease in temperature in the monoclinic phase. In SCT-0.03, the splitting of the TA phonon (TA1 and TA2 in the tetragonal phase) at T ct and the temperature dependence of TA2 mode in the ferroelastic phase appears to be quite similar to that observed in Cs 2 CdBr 4 . We therefore tentatively ascribe the TA2 phonon as the quasitransverse mode propagating along the c axis of SCT-0.03 in the tetragonal phase. However, the correct description of TA2 mode requires detailed Brillouin light scattering studies in various polarization geometries.
caused by the interaction between the phonon mode and the order parameter. It can be seen from Figure 17 that the order parameter interaction with the LA phonon does not agree with the bilinear coupling that requires |T −T ct | -1 -type temperature dependence on either side of the phase transition. Slonczewski and Thomas reported that the ferroelastic phase transition at 110 K in SrTiO 3 can be described by coupling between the strain S and the order parameter Q in the form Fc = SQ 2 , where the order parameter Q is the soft mode coordinate and  is the coupling coefficient [87]. The ferroelastic transition in SrTiO 3 could be well described by Landau theory, taking into account the interaction between the strain and the soft optical phonon mode.We suggest that in SCT also, the observed acoustic anomaly could be due to the coupling between the LA phonon and the soft mode. The temperature dependence of the relaxation time  LA of the LA phonon, calculated using Equation (15), is plotted in Figure 18. The relaxation time of the LA phonon exhibits critical slowing when approaching T ct from both sides of the transition.
The value of  o obtained from Equation (16) is 0.077 ps below T ct and 0.057 ps above T ct . This is about half of the value of KDP and TGS and is very similar to that obtained in PbHPO 4 [75]. This suggests the weaker coupling of the LA phonon with the soft mode of the paraelectric host. The temperature dependence of the frequency and width of the TA1 and TA2 phonons are shown in Figures 19 and 20. The frequency and the width of the TA1 phonon mode exhibit similar behavior as that of the LA phonon mode. However, a large scattering of the data observed in the width of the TA1 phonon mode limits further detailed discussion on its temperature dependence. The TA2 mode appears just near T ct . The experimental data, however, could be unambiguously fitted only below 265 K. The TA2 mode shows a significant softening with the decrease in temperature in the ferroelastic phase. In the Brillouin scattering study of phase transition in Cs 2 CdBr 4 , the splitting of TA mode into two components was observed on the phase transition from orthorhombic to monoclinic phase [88]. One of the components was purely TA mode (C 66 ) corresponding to the parent orthorhombic phase, while the other was ascribed to the quasi-transverse mode  4 (b), propagating along the b axis of the crystal. The frequency of the  4 (b) mode was observed to decrease with decrease in temperature in the monoclinic phase. In SCT-0.03, the splitting of the TA phonon (TA1 and TA2 in the tetragonal phase) at T ct and the temperature dependence of TA2 mode in the ferroelastic phase appears to be quite

Conclusions
Brillouin light scattering studies on the relaxation dynamics of disordered perovskite compounds have been reviewed. In PMT, the elastic constant C 11 clearly reveals the existence of two distinct dynamical regimes of polarization characterized by the Burns

Conclusions
Brillouin light scattering studies on the relaxation dynamics of disordered perovskite compounds have been reviewed. In PMT, the elastic constant C 11 clearly reveals the existence of two distinct dynamical regimes of polarization characterized by the Burns temperature T B and the intermediate temperature T*. In the quasi-dynamic phase between T B and T*, the elastic constant obeys the mean field model and softens as (T−T o ) 0.5 . The reduced softening of the elastic constant C 11 and the increased attenuation in the quasistatic region below T*are attributed to the strong coupling between the local polarization and strain in PNRs. The temperature dependence of the relaxation time of the central peak was found to be non-Arrhenius. The observed critical slowing of the relaxation time between T* and T f is attributed to dipole flipping and the volume fluctuations of PNRs below T f . In the partially disordered PST, a sharp anomaly was observed in the frequency and the width of the LA phonon upon the transition from cubic relaxor to rhombohedral ferroelectric phase at 295 K. The calculated relaxation time of the LA phonon suggested the order-disorder nature of the structural phase transition. In 0.65PIN-0.35PT single crystal, anomalous behavior in the frequency and the width of the LA phonon was observed at cubic-tetragonal and tetragonal-rhombohedral phase transition temperatures. Twocomponent quasi-elastic (narrow and broad) central peaks were suggested to originate from two types of PNRs with different local structures. While the rhombohedral-type PNRs were pertinent to the parent relaxor PIN, the tetragonal-type PNRs areinduced by PT substitution. The broad central peak was suggested to associate with the relaxation of rhombohedral PNRs, while that of the narrow central peak to the relaxation of tetragonal PNRs. The relaxation time of the LA phonon was found to be qualitatively similar to that of the broad quasi-elastic central peak. The phase transition from cubic to tetragonal structure was suggested to originate from nucleation and growth of tetragonal PNRs. In SCT-0.03, the LA phonon showed critical anomalies at the structural phase transition temperature T ct = 280 K. The Landau-Khalatnikov-type asymmetry was observed in the width of the LA phonon. The relaxation time of the LA phonon showed a critical slowing in the vicinity of the transition.The observed acoustic anomaly was suggested to be due to the coupling between the LA phonon and the soft mode of SrTiO 3 . In the ferroelastic phase, the quasitransverse TA (TA2) phonon mode was found to exhibit softening behavior with decrease in temperature.
Funding: This research received no external funding.

Informed Consent Statement: Not applicable.
Data Availability Statement: The data sets generated during and/or analyzed during the current study are available from the corresponding author on request.