Open access

Non-Linear Dielectric Response of Ferroelectrics, Relaxors and Dipolar Glasses

Written By

Seweryn Miga, Jan Dec and Wolfgang Kleemann

Submitted: 08 November 2010 Published: 23 August 2011

DOI: 10.5772/20380

From the Edited Volume

Ferroelectrics - Characterization and Modeling

Edited by Mickaël Lallart

Chapter metrics overview

4,511 Chapter Downloads

View Full Metrics

1. Introduction

The dielectric response of dielectrics with respect to temperature, pressure, frequency and amplitude of the probing electric field is an essential issue in dielectric physics (Jonscher, 1983). This concerns both normal dielectrics and those specified as ferro- or antiferroelectrics (Lines & Glass, 1977). This basic phenomenon of dielectric materials has been extensively studied in the literature both experimentally and theoretically, however, mainly restricting to the linear dielectric response where a linear relationship between polarization, P, and external electric field, E,


is fulfilled. Here ε0 and χ1 stand for the electric permittivity of the free space (dielectric constant) and the linear electric susceptibility of a given dielectric material, respectively. In the framework of this approach the analysis of the dielectric relaxation is possible via the Debye model (including the Cole-Cole equation and related formulae in case of multidispersive behavior) and the relaxation rate can often be described by an Arrhenius relation (Jonscher, 1983; Lines & Glass, 1977; Debye, 1929; von Hippel, 1954; Böttcher, 1973; Kremer & Schönhals, 2003) or by an activated dynamic scaling law (Fisher, 1986; Kleemann et al., 2002) which in the particular case of the exponent θν = 1 converts into the well known Vogel-Fulcher-Tammann relation (Jonscher, 1983; Kremer & Schönhals, 2003). More complex distribution functions of relaxation-times are used as well (Jonscher, 1983; Böttcher, 1973; Kremer & Schönhals, 2003).

At higher electric fields Eq. (1) becomes violated and the relation between polarization and electric field is better represented by a power series of E,


which contains higher order terms with respect to the external electric field, where χi with i > 1 are the non-linear susceptibilities: second-order, third-order and so on (Böttcher, 1973). These non-linear components are referred to as hypersusceptibilities and the non-linear contribution to the polarization response is designated as the hyperpolarization (Jonscher, 1983). While the higher-order susceptibilities contain a wealth of important information (Wei & Yao, 2006a; Wei & Yao, 2006b), the hyperpolarization has nevertheless been less studied due to obstacles like (i) the difficulty of the respective dielectric measurements, because the nonlinear signal is usually some orders of magnitude smaller than the linear response, and (ii) lack of an appropriate theory to deal with dielectric spectra as a function of the electric field (Chen & Zhi, 2004).

In order to overcome the technological challenge we have constructed a fully automatized ac susceptometer for simultaneous measurements of the phase resolved complex linear and complex non-linear ac susceptibilities of lossy and dispersive dielectric materials (Miga, Dec & Kleemann, 2007). It allows measurements over a wide range of experimental variables, such as ac amplitudes up to 40 V, frequencies from 10-2 to 103 Hz, and temperatures from 100 K to 600 K utilizing only current/voltage and analogue/digital converters and a computer. In contrast to the commonly used analysis of the charge accumulated on a standard capacitor in series with the sample our method is based on the analysis of the current flowing directly through the sample. Absence of any capacitive voltage dividers in the measurement circuit eliminates uncontrolled phase shifts. That is why the instrument provides high quality nonlinear susceptibility data and in particular appears as a very convenient tool for discrimination between continuous and discontinuous phase transitions when determining the sign of the real part of the third order dielectric susceptibility.

We have applied this new instrument to various basic ferroelectric scenarios, such as the classic first- and second-order ferroelectric transitions of barium titanate (BaTiO3) (Miga & Dec, 2008), triglycine sulphate (TGS) (Miga & Dec, 2008) and lead germanate (Pb5Ge3O11) (Miga & Dec, 2008), the double anomalous second-order transitions of Rochelle salt (Miga et al., 2010a), the smeared transition of the classic relaxor ferroelectrics lead magno-niobate (PbMg1/3Nb2/3O3, PMN) (Dec et al., 2008) and strontium-barium niobate (Sr0.61Ba0.39Nb2O6, SBN61) (Miga & Dec, 2008), the dipolar glassy and ferroelectric transitions of Li+-doped potassium tantalate (K1-xLixTaO3, KLT) with x = 0.005, 0.011 and 0.063 (Dec et al., 2010; Miga et al., 2010b).


2. Theoretical background of nonlinear dielectric response

Dielectric properties of materials are usually investigated via linear dielectric response. In this case a linear relationship, Eq. (1), between polarization, P, and external electric field, E, is fulfilled. At higher field intensities, the polarization may be a non-linear function of the electric field strength. For not too high electric field strength, one can present the polarization as the power series expansion in the variable E, Eq. (2). For symmetry reasons the second-order dielectric susceptibility is nonzero only for macroscopically noncentrosymmetric systems.

Let us consider the Landau-Ginzburg-Devonshire (LGD) theory of ferroelectric phase transitions (PT) (Ginzburg, 1945; Devonshire, 1949). According to this theory the free energy density G of a ferroelectric material within its paraelectric phase is represented as a power series expansion with respect to the polarization P,


where G0 stands for the free energy at P = 0, A and T0 are constants. B and C are usually smooth functions of temperature (Fujimoto, 2003). In the ferroelectric phase the polarization P, yields the spontaneous one, Ps, as an order parameter. The sign of B determines the kind of the PT. For positive B a continuous (second-order) PT occurs. On the other hand, for negative B a discontinuous (first-order) PT appears. Let us first consider a continuous PT to occur at T0 = Tc. In this case one can drop the last term in Eq. (3) without loss of generality. Thermodynamic equilibrium requires fulfilment of the condition


This relationship between E and P denotes the electric equation of state. After successive differentiation of this equation with respect to P one obtains the following susceptibilities (Ikeda et al., 1987):


The second-order susceptibility is proportional to the polarization P, therefore it changes sign when the polarization changes direction. Additionally 2 vanishes when the polarization of the sample vanishes. That is why 2 is a sensitive probe of the net polarization, but it is unsuitable for determination of the PT order. In contrast to 2 the odd order susceptibilities 1 and 3 depend on the square of the polarization, which makes them insensitive to the orientation of the polarization. Moreover, 1 and 3 do not vanish, even if the polarization is equal to zero. For classic ferroelectrics on heating the spontaneous polarization Ps vanishes at the PT point and is zero within the paraelectric phase, where Eq. (7) simplifies to


Due to positive B, 3 is negative above a continuous PT point. Within the ferroelectric phase P=Ps andPS2=(A/B)(TcT), hence,


Within the ferroelectric phase 3 has a positive sign. Thus the LDG theory predicts a change of sign of 3 at a continuous PT. The scaling theory (Stanley, 1971) predicts that B scales as B=B0|τ|γ2β, whereτ=(TTc)/Tc. and are critical exponents of the linear susceptibility and of the order parameter, respectively. For the Landau universality class, where - 2 = 0 (i. e. = 1, = ½) one expects a temperature independent value of B within the paraelectric phase.

A similar calculation for a discontinuous PTs, where B < 0 and C > 0, yields 3 > 0 at all temperatures, in particular also above PT (Ikeda et al., 1987). Fortunately from an experimental point of view 3 is given by this same equation (8) within paraelectric phase independent of the ferroelectric PT order. The sign of 3 is a sensitive probe for discrimination between continuous and discontinuous ferroelectric PTs.

In all known ferroelectrics the paraelectric phase is located above the stability range of the ferroelectric one. However, sodium potassium tartrate tetrahydrate (Rochelle salt, RS) (Valasek, 1920, 1921) apart from a classic high temperature paraelectric phase has an additional, unusual one located below the ferroelectric phase. Both PTs, between the paraelectric phases and the ferroelectric one have continuous character. In order to predict the sign of 3 within the low-temperature paraelectric phase one can refer to the theory of Mitsui (Mitsui, 1958). Within this theory the electric equation of state for RS is similar to Eq. (4) with a positive coefficient of the cubic term, P3. Therefore a negative sign of 3 is expected within low-temperature paraelectric phase (Miga et al., 2010a).

Another class of dielectrics are relaxor ferroelectrics. They are usually considered as structurally disordered polar materials, which are characterized by the occurrence of polar nanoregions (PNRs) of variant size below the so-called Burns temperature, Td (Burns & Dacol, 1983) far above the ferroelectric Curie temperature, Tc. In contrast to conventional ferroelectrics, relaxors do not exhibit any spontaneous macroscopic symmetry breaking. In addition they are characterized by a large, broad and frequency-dependent peak in the temperature dependence of dielectric susceptibility. According to the spherical random-bond-random-field (SRBRF) model (Pirc & Blinc, 1999) the dipole moments of individual nanopolar clusters interact via a spin-glass-type random exchange coupling, and are subject to quenched random local electric fields. The model Hamiltonian of such a system is formally written as


where Jij are glass-like random intercluster couplings or bonds, hiand E random local and uniform external electric fields, respectively. In the dynamic approximation it is assumed that PNRs reorient by means of stochastic flips described by a relaxation time . This model yields negative 3 with two extremes observed at the freezing temperature, Tf, and at the peak temperature of 1, Tm.

The measured values of 1 and 3 can be used for calculating the so-called scaled non-linear susceptibility, a3, which is given by (Pirc & Blinc, 1999)


Within the paraelectric phase of classic ferroelectrics a3 is equal to the nonlinearity coefficient B, cf. Eqs. (8) and (11). For ferroelectics displaying a continuous PT, a3 = -8B within ferroelectric phase. The SRBRF model yields negative 3 and positive a3 with two extremes observed at the freezing temperature, Tf, and at the peak temperature of 1, Tm. It should be noticed that systems like dipolar glasses a3 are also expected to exhibit a critical singularity at Tf, where 3 (T − Tf )- with the mean field exponent = 1 (Pirc et al., 1994).

A schematic comparison of predictions of the above theories is presented in Fig. 1.

Figure 1.

Schematic presentation of linear and nonlinear responses of classic ferroelectrics with (a) continuous PTs, (b) discontinuous PTs, and (c) relaxor ferroelectrics (see text).


3. Methods of measurement of nonlinear dielectric response

Nonlinear dielectric response can be measured using two different kinds of experimental methods. The first one is based on the investigation of the ac linear dielectric susceptibility as a function of a dc bias field. A schematic presentation of this method is shown in Fig. 2a. One applies to the sample a weak probing ac electric field with fixed amplitude (that warrants a linear response) and a superimposed variable dc bias field, EB. The bias field amplitude may reach values up to 5.5 107 Vm-1 (Leont’ev et al., 2003). In barium titanate a strong dc field induces the PT between paraelectric and ferroelectric phases (Wang et al., 2006). By changing the value of the dc electric field the local slope of the polarization curve is probed in many points. For such a kind of experiment one can use commercially available LRC meters or impedance analyzers, e.g. Agilent E4980A or Solartron 1260. The electric field dependence of the linear susceptibility of ferroelectrics displaying continuous PT fulfils the following relations (Mierzwa et al., 1998):


for the paraelectric phase and


for the ferroelectric one. (E) and (0) are the susceptibilities for bias, E, and zero electric field, respectively. By use of Eqs. (12) and (13) one can calculate the nonlinearity coefficient B. Its knowledge allows us to calculate the third-order nonlinear susceptibility 3. Unfortunately, the above described method has at least one restriction when investigating the nonlinear dielectric response. Namely, during so-called field heating/cooling runs, unwanted poling and remnant polarization of the investigated sample can evolve under a high bias field.

Figure 2.

Presentation of methods of measurement of nonlinear dielectric response using (a) a weak probing ac electric field with fixed amplitude and a superimposed variable dc bias field, and (b) an enhanced ac field. The effects are exaggerated for visualisation.

The second kind of method is free from this restriction. A schematic presentation of this method is shown in Fig. 2b. During the experiment the sample is exposed to an ac probing field with sufficiently large amplitude. Consequently, under this condition the temperature dependences of the linear and nonlinear susceptibilities are determined under zero dc field in heating and cooling runs. That is why the nonlinear susceptibility detected this way can be referred to as a dynamic nonlinear susceptibility related to ac dielectric nonlinearity. Usually the amplitude of this field is much smaller than the bias field strength used in the above described experiment. As a result of the nonlinear P(E) dependence, the output no longer remains harmonic. The distorted signal may be subjected to Fourier analysis revealing all harmonic components in the polarization response. This is the main idea of our nonlinear ac susceptometer (Miga et al., 2007). Using harmonics of displacement current density ji one can calculate the linear and nonlinear dielectric susceptibilities i as follows:


where and E0 are the angular frequency and the amplitude of the applied electric field respectively. In Eq. 14 the terms up to seventh order are involved, hence, susceptibilities up to the fifth order can be regarded as reliable even for strongly nonlinear materials. Neglecting harmonics higher than the order of a considered susceptibility may lead to artificial effects. The next important point is the simultaneous measurement of all displacement current components, which considerably improves the accuracy of the measured susceptibilities (Bobnar et al., 2000). In case that the phase shifts of the displacement current harmonics are known it is possible to calculate the real, χi, and imaginary, χi, parts of all susceptibilities. This is impossible, when measurements are done with a dc bias electric field. It should be noted that up to now commercial instruments for the dynamic measurement of complex nonlinear dielectric susceptibilities are unavailable.


4. Experimental results

4.1. Ferroelectrics displaying continuous PTs

Triglycine sulphate, (NH2CH2COOH)3H2SO4 (TGS) is a model ferroelectric displaying a continuous PT. Sodium potassium tartrate tetrahydrate (Rochelle salt, RS), NaKC4H4O64H2O, and lead germanate, Pb5Ge3O11 (LGO) exhibit continuous PTs, as well. Despite the quite different structures of the crystals and different mechanisms of their PTs, a negative sign of the real part of the third-order nonlinear susceptibility is expected in their paraelectric phases. Fig. 3 shows the temperature dependences of the real parts of the linear, the third order nonlinear dielectric susceptibilities, and the scaled susceptibility, a3. The linear susceptibility of a TGS crystal (Fig. 3a) obeys a Curie-Weiss law within the paraelectric phase very well with a critical exponent =1.0000.006. According to predictions of the phenomenological theory of ferroelectric PT, 3’ changes its sign at the PT point (Fig. 3b).

Figure 3.

Temperature dependences of the real parts of the linear (a) and third order non-linear (b) susceptibilities and a3 coefficient (c) of TGS. The amplitude of the probing ac electric field was 5kVm-1.

While below Tc the third-order susceptibility has positive values, it is negative above Tc. As was mentioned earlier such change in the sign of 3’ is one of the primary features of the continuous ferroelectric PT. The scaled non-linear susceptibility has been evaluated using Eq. (11). The measured values of 1’ and 3’ were used for calculations. The temperature dependence of a3 of TGS crystal is shown in Fig. 3c. Within the paraelectric phase a3 is equal to B and has a positive value. This is an attribute of the continuous PT. Moreover, a3 is practically temperature independent within the temperature range of Tc + 0.2 K < T < Tc + 2 K above the PT point. This behaviour is consistent with predictions of the scaling theory and confirms that TGS belongs to the Landau universality class. Weak temperature dependence of a3 and 1/|T-Tc| dependence of 1’ leads to a 1/|T-Tc|4 anomaly of 3’ (see Eq. 11). Therefore the temperature dependence of 3’ is much sharper than that of 1’. Consequently, 3’ rapidly vanishes in the surrounding of the PT. An increase of a3 as observed within the paraelectric phase close to Tc is probably due to crystal defects. A similar effect was reported for a -ray damaged TGS crystal (Cach, 1988). Then within the ferroelectric phase a3 changes slowly due to the contribution of the domain walls to the resultant dielectric response. In the ferroelectric state the measured values of the linear and non-linear susceptibilities, which were used for calculation, involved the responses not only from the ‘‘pure’’ ferroelectric system but also from the domain walls, which were forced to move under the probing electric field.

Figure 4.

Temperature dependences of the real parts of the linear (a), the second (b) and the third (c) order nonlinear susceptibilities of Rochelle salt. The amplitude of the probing ac electric field was 500 Vm-1.

Fig. 4 shows the temperature dependences of the linear, the second, and the third-order nonlinear susceptibilities of Rochelle salt. The temperature dependences of 1’ and 3’ are, close to the high temperature ferroelectric-paraelectric PT ( 297.4 K), qualitatively similar to those of the TGS crystal. As was mentioned earlier, RS displays an additional low temperature PT. This transition between the ferroelectric and the low temperature paraelectric phase appears at about 254.5 K. At this PT 3’ changes its sign as compared to the high temperature PT (Fig. 4c). As a result of the inverse order of phases, 3’ is negative below and positive above the low temperature PT. In this way the third-order nonlinear susceptibility is negative in both paraelectric phases close to PT points. Fig. 4b shows the temperature dependence of the second-order nonlinear susceptibility. The sign of this susceptibility depends on the net polarization orientation. Therefore it can be changed by polarization of the sample in the opposite direction. The simplest way to change the sign of 2’ is a change of the wires connecting the sample to the measuring setup. So, in contrast to sign of the 3’, the sign of the second order nonlinear susceptibility is not very important. In the case of 2’ most important is its nonzero value and the observed change of sign of this susceptibility at 273 K. This hints at a modification of the domain structure within the ferroelectric phase. The next change of sign at the low temperature PT point is originating from different sources of the net polarization above and below this point. Above 254.5 K this polarization comes from the uncompensated ferroelectric domain structure, whereas within the low-temperature paraelectric phase it originates merely from charges screening the spontaneous polarization within the ferroelectric phase.

Lead germanate (LGO) displays all peculiarities of a continuous ferroelectric PT (see Fig. 5 a, c, e) (Miga et al., 2006). However, a small amount of barium dopant changes this scenario (Miga et al., 2008). Ba2+ ions replacing the host Pb2+ influence the dielectric properties. 2% of barium dopant causes a decrease of the linear susceptibility, broadening of the temperature dependence of 1’, and a decrease of the PT temperature. Despite all these changes the temperature dependences of 1’ for pure and barium doped LGO are qualitatively similar. A completely different situation occurs, when one inspects the nonlinear dielectric response. Small amounts of barium dopants radically change the temperature dependence of the third-order nonlinear susceptibility (Fig.5 b). Similarly to linear one the anomaly of the third-order nonlinear susceptibility shifts towards lower temperature and decreases. The most important difference is the lacking change of sign of 3’. In contrast to pure LGO for barium doped LGO the third-order nonlinear susceptibility is positive in the whole temperature range. Therefore one of the main signatures of classic continuous ferroelectric PTs is absent. Due to the positive value of 3’ the scaled nonlinear susceptibility is negative in the whole temperature range (Fig. 5d). No change of a3 is observed. This example shows the high sensitivity of the nonlinear dielectric response to the character of the PT. In the discussed case a change of character of the PT is due to the presence of barium induced polar nanoregions (PNRs). The occurrence of PNRs results in weak relaxor properties of barium doped LGO.

4.2. Ferroelectrics displaying discontinuous PT

Barium titanate, BaTiO3 (BT), is a model ferroelectric displaying three discontinuous ferroelectric PTs (von Hippel, 1950). Two of them appear – at rising temperatures - between rhombohedral, orthorhombic and tetragonal ferroelectric phases at about 200 K and 280 K respectively. The final discontinuous PT appears between the ferroelectric tetragonal and the cubic paraelectric phase at about 400 K. Fig. 6 shows the temperature dependences of the real parts of the linear (a) and third-order nonlinear (b) susceptibilities, and the a3 coefficient (c) of a BT crystal in the vicinity of the ferroelectric-paraelectric PT. The amplitude of a probing ac electric field was equal to 7.5 kVm-1.

Figure 5.

Temperature dependences of the real part of the linear dielectric susceptibility (a) for LGO:Ba 2% and LGO, real part of third-order nonlinear dielectric susceptibility for (b) LGO: Ba 2% and (c) LGO, and scaled nonlinear susceptibility a3 for (d) LGO:Ba 2%, and (e) LGO crystals. The amplitude of the probing ac electric field was 15 kVm-1.

Figure 6.

Temperature dependences of the real parts of the linear (a) and third order non-linear (b) susceptibilities and of the a3 coefficient (c) of a BaTiO3 crystal. The amplitude of the probing ac electric field was 7.5kVm-1.

Temperature hysteresis is one of the typical features of discontinuous PT. Therefore for viewing this phenomenon Fig. 6 presents results for cooling and heating runs. The third-order nonlinear susceptibility is positive, both within the ferroelectric and the paraelectric phases (Fig. 6b). Just below Tc a rapid decrease of the third-order dielectric susceptibility is observed. For that reason, the value of 3’ is much smaller within the paraelectric phase than in the ferroelectric one, yet it is still positive. The positive sign of 3’ as detected in BT is consistent with predictions of the phenomenological theory of discontinuous ferroelectric PTs. Unlike the scaled nonlinear susceptibilities of crystals displaying a continuous PT, a3 of BaTiO3 is negative in the whole investigated temperature range (Fig. 6c). However, similarly to TGS (Fig. 3) or RS (Fig. 4) a jump-like change of a3 occurs at the PT. Within the paraelectric phase a3 is equal to B (see section 4.1) and displays weak temperature dependence. Wang et al. (Wang et al., 2007) proposed incorporation of higher order terms (up to the eighth power) in the Landau potential for the temperature independence of the experimentally obtained B coefficient. This proposition was given for analyzing the data collected with a bias field up to 1500 kV/m. Fig. 6 shows data collected for a two hundred times smaller electric field, for which terms at such high order are not important. They cannot explain the temperature dependence of the nonlinearity coefficient B. Hence, although the nonlinear properties of BaTiO3 have been investigated already for six decades, the question of the temperature dependence of the nonlinearity parameter B is still open.

4.3. Relaxor ferroelectrics

The classic relaxor ferroelectric lead magno-niobate (PbMg1/3Nb2/3O3, PMN) (Smolenskii et al., 1960) displays an average cubic structure in the whole temperature range (Bonneau et al., 1991). No spontaneous macroscopic symmetry breaking is observed in this relaxor. On the other hand and in contrast to PMN, the relaxor strontium-barium niobate Sr0.61Ba0.39Nb2O6 (SBN61) spontaneously undergoes a structural phase transition. On cooling the vanishing of a mirror plane leads to a lowering of the tetragonal symmetry from 4/mmm to 4mm at a temperature of about 348 K (Oliver et al., 1988).

Figure 7.

Temperature dependences of 1’ (a), 2’ (b), 3’ (c) and a3 (d) of a PMN crystal measured at f = 31, 100, 319, and 1000 Hz. A probing ac electric field with an amplitude of 12 kV/m was applied along the [100] direction.

Figure 8.

Temperature dependences of 1’ (a), 2’ (b), 3’ (c) and a3 (d) of a SBN61 crystal measured at f = 10, 31, 100, 318, and 1000 Hz. A probing ac electric field with an amplitude of 7.5 kV/m was applied along the [001] direction.

Fig. 7 presents the temperature dependences of the real parts of the linear, 1’ (a), second-order, 2’ (b), and third-order, 3’ (c), dielectric susceptibilities and of the scaled nonlinear susceptibility a3 (d) of the PMN single crystal recorded along [100] direction (Dec et al., 2008). The linear susceptibility displays features of a relaxor i.e. a large, broad and frequency-dependent peak in the temperature dependence. Nonzero 2’ as presented in Fig. 7b is an indicator of net polarization of the crystal. The presence of this polarization was independently confirmed by measurements of thermo-stimulated pyroelectric current of an unpoled sample. Integration of this current indicates an approximate value of the average polarization as low as 310-5 C/m2. This polarization is much smaller than the spontaneous polarization of ferroelectrics, but is well detectable. The observed pyroelectric response and nonzero 2’ hints at an incomplete averaging to zero of the total polarization of the PNR subsystem. Fig. 7c shows the temperature dependence of the third-order nonlinear susceptibility. In contrast to the predictions of the SRBRF model, 3’ is positive in the whole temperature range. The positive sign of 3’ may result from a term 180BP21 exceeding unity in Eq. 7. Positive sign of 3’ results in a negative sign of a3. Consequently the sign of a3 differs from that predicted by the SRBRF model. Having in mind that corrections due to the fifth harmonic contribution produce large noise at temperatures below 210 K and above 310 K (Fig. 7d, main panel), less noisy a3 data calculated only from first and third harmonics are presented in the inset to Fig. 7d. Since both 1’ and 3’ do not display any anomalies in the vicinity of the freezing temperature Tf 220 K, also a3 does not exhibit any maximum in contrast to SRBRF model predictions. It continues to increase monotonically by almost two orders of magnitude when decreasing the temperature from 260 to 180 K. This result is independent of the number of harmonics used for experimental data analysis. It is worth stressing that despite the remarkable dispersion of 1’ and 3’, the scaled susceptibility a3 does not display any sizable frequency dependence, even around the temperatures of the susceptibility peaks. It may, hence, be considered as a static quantity (Glazounov & Tagantsev, 2000). According to our result the dispersion of a3 is even weaker than reported previously (Glazounov & Tagantsev, 2000).

Fig. 8 shows results of measurements of the linear and nonlinear dielectric response of SBN61 crystal (Miga & Dec, 2008). The probing electric ac field was applied along [001], which is the direction of the polar axis below Tc. The peak temperatures of the linear susceptibility (Fig. 8a) are a few degrees above the temperature of the structural phase transition. Fig. 8b presents the temperature dependence of the second-order susceptibility. Similarly to results obtained on PMN this susceptibility is non-vanishing. However, the values of 2’ of SBN61 are almost two orders of magnitude larger than those measured on PMN. Therefore 2’ is detectable even forty degrees above Tc. Analysis of the thermo-stimulated current indicates a net polarization of a nominally unpoled SBN61 crystal equal to 510-2 C/m2. The higher value of this polarization (in comparison with PMN) results in higher values of 2’. Fig. 8c shows the temperature dependence of 3’. This susceptibility is positive within both the ferroelectric and the paraelectric phase. Again, the sign of 3’ disagrees with predictions of the SRBRF model. As was discussed above, the positive sign of 3’ presumably originates from the presence of net polarization (detected by 2’, see Fig. 8b). The disagreeing sign of 3’ leads to a disagreement of the sign of a3 as well (Fig. 8d). The scaled nonlinear susceptibility shows an anomaly related to the phase transition, but it does not exhibit any additional peak as predicted by the SRBRF model.

The results obtained for both relaxor ferroelectrics are qualitatively similar. Therefore, they are independent of the presence or absence of a structural phase transition and macroscopic symmetry breaking. As discussed in Section 2 the dielectric properties of relaxors are mainly determined by PNRs, which were detected in both of the above relaxors. Unfortunately, the early version of the SRBRF model (Pirc et al., 1994) predicts a negative sign of the third-order nonlinear dielectric susceptibility, which is not confirmed in experiments. Consequently the sign of the scaled nonlinear susceptibility a3 is incorrect as well. Very probably this unexpected result is due to the fact that the PNRs primary do not flip under the ac electric field, but merely change their shape and shift their centers of gravity in the sense of a breathing mode (Kleemann et al., 2011).

4.4. Dipolar glasses

4.4.1. Orientational glasses

The formation of dipolar glasses in incipient ferroelectrics with perovskite structure, ABO3, such as SrTiO3 and KTaO3, by A-site substitution with small cations at low concentrations has been a fruitful topic since more than 20 years (Vugmeister & Glinchuk, 1990). For a long time probably the best-known example has been the impurity system K1-xLixTaO3 (KLT for short) with x<< 1, whose complex polar behavior is known to be due to the interaction of the (nearly) softened transverse-optic mode of the host-lattice and the impurity dynamics (Höchli et al., 1990). Fig. 9 shows the structure model of A-site substituted Li+ viewing the nearest neighbor environment in the KTaO3 lattice from two different perspectives. At very low concentrations, x ≈ 0.01, it reveals signatures of glasslike behavior (Höchli, 1982; Wickenhöfer et al., 1991), while a ferroelectric ground state with inherent domain structure is encountered at higher concentrations, x ≥ 0.022 (Kleemann et al., 1987).

Figure 9.

Displacement vectors (blue arrows) of an off-center Li+ (Mn2+) ion in A-site doped KTaO3 (SrTiO3) viewed (a) from the center of 12 surrounding oxygen ions and (b) from the corner of the elementary cell.

Only recently a similar system has been discovered with qualitatively new properties. The impurity system Sr1-xMnxTiO3 (SMnT for short; Fig. 9) reveals very similar dipolar glassy properties, but the additional magnetic degrees of freedom of the Mn2+ dopant enable simultaneously a spin glass state (Shvartsman et al., 2008). This unique ‛multiglass’ situation has paved the way to a new materials class: ‛disordered multiferroics’ (Kleemann et al., 2008).

In addition to conventional tests of the glass transition, e. g. by verifying the divergence of the polar relaxation times, the behavior of the nonlinear susceptibility is believed to similarly decisive. As was first acknowledged in spin glass physics (Binder & Young, 1986), but later on also in the field of orientational glasses (Binder & Reger, 1992), criticality at the glass temperature, Tg, is expected to give rise to a divergence of the third-order nonlinear susceptibility,χ33P/E3 , where the polarization P denotes the homogeneous order parameter. Recent experimental attempts will be reported below.

4.4.2. Orientational glass K0.989Li0.011TaO3

The experiments on KLT were performed on a Czochralski grown single crystal sample with x = 0.011 with dimensions 3×2×0.5 mm3 and (100) surfaces (Kleemann et al., 1987). Dipolar relaxation was studied as a function of temperature T via measurements of the complex dielectric susceptibility, = ′ - iχ″ vs. T, by use of different experimental methods adapted to different frequency ranges, 10-3 ≤ f ≤ 106 Hz. They included a Solartron 1260 impedance analyzer with 1296 dielectric interface (Fig. 10a and b) and a digital lock-in analyzer (Wickenhöfer et al., 1991) (Fig. 10c) for linear, and a homemade computer-controlled digital susceptometer (Miga et al., 2007) for non-linear dielectric susceptibility data (Fig. 11) at high precision under relatively low excitation voltages.

Fig. 10a and b show dielectric susceptibility data, χ′(T) and χ″(T), for frequencies 10-3 < f < 106 Hz, which reveal various signatures of glassy behaviour. The peak of χ′(T) in Fig. 10a converges toward a finite glass temperature, Tm → Tg = (33.6 ± 0.1) K (Wickenhöfer et al., 1991) following critical dynamics with a power lawτ=1/(2πf)=τ0εzν, where ε=(Tm/Tg1)>0 is the reduced temperature and zν=6.6±0.2 is the dynamical critical exponent. Essentially the same (static) glass temperature, Tg = (33.5±1.5) K, emerges for the ‛largest’ relaxation time τm as defined by the condition χ(fmin,T)=0.05χmax(T) from frequency spectra of the loss function χ(f) in Fig. 10c (Kleemann et al., 2011). In this figure it is clearly seen, how the center of gravity of the loss spectrum shifts to very low frequencies and its width broadens toward very long relaxation times, τm, as T→ Tg.

Figure 10.

Dielectric susceptibility, χ′(T) (a), χ″(T) (b), and χ″(f) (c) of K0.989Li0.011TaO3 measured at frequencies 10-3 < f < 106 Hz and temperatures 10 < T < 90 K. Decades of f and steps of ∆T = 2.5 K are parameters in (a, b) and (c), respectively.

Figure 11.

Nonlinear susceptibility χ3'(T) of K0.989Li0.011TaO3 measured at frequencies 11 < f < 3000 Hz.

A complementary test of glassy criticality refers to the third-order nonlinear susceptibility,χ33P/E3. It fulfils the expectation of a divergence at Tg (Binder & Reger, 1992) only partially. Fig. 11 shows χ3vs.T for frequencies 11 Hz ≤ f ≤ 3 kHz on cooling with dT/dt = -0.2K/min under an ac field amplitude E0 = 45 kV/m and by analyzing the emerging signal up to the fifth harmonic (Kleemann et al., 2011). While χ3 drops below zero on the high-T edge at low frequencies, all signals are dominated by positive low-T peaks. This unexpected result is very probably due to the fact that the dipolar glass KLT is only insufficiently modeled by a pseudospin system in full analogy to an Ising spin glass (Binder & Reger, 1992). Actually we deal with a frustrated system of nanoclusters (Vugmeister & Glinchuk, 1990), which are subject to complex dipolar interactions and underlie internal dynamical degrees of freedom. As a consequence the nanodipolar clusters do not primarily flip under the external electric ac field, but merely change their shape and shift their centers of gravity in the sense of a breathing mode. Obviously this action becomes more effective the lower the temperature (i.e. the higher the lattice permittivity) and the lower the frequency. Thus the expected negative divergence is replaced by a positive peak nearly coinciding with that of the linear susceptibility, χ1(Fig. 10a).

4.4.3. Multiglass Sr0.98Mn0.02TiO3

The experiments on Sr0.98Mn0.02TiO3 were performed on a ceramic sample prepared by a mixed oxide technology (Tkach et al., 2005). Preponderant incorporation of Mn2+ onto A-sites of the perovskite structure (Fig. 9) was confirmed by energy dispersive X-ray spectra (Tkach et al., 2006), Mn2+ ESR analysis (Laguta et al., 2007), and EXAFS spectroscopy (Lebedev et al., 2009; Levin et al., 2010). Fig. 12 shows the components of the complex dielectric susceptibility χ=χiχrecorded at frequencies 10-3 f 106 Hz and temperatures 10 T 100 K as χ(T) (a) and χ(f) (c). The broad and strongly frequency dependent peaks of both components are related to the dynamics of polar clusters created by off-center displacements of Mn2+ cations (Tkach et al., 2007). The position of the peak temperatures Tm in Fig. 12a is well described by a power law of the respective frequency, f(Tm)(Tm/Tg1)zν, which is a typical manifestation of glassy critical behavior (Shvartsman et al., 2008). Best fits of the experimental data yield the glass temperature Tg = 38.3 0.3 K (Fig. 12b) and the dynamical critical exponent, z = 8.5 0.2, which compares well with that of spin glasses (Jönsson, 2004). Obviously the dynamics of the polar clusters becomes frozen at Tg, where the relaxation time τ=1/(2πf) diverges on a percolating network. This defines a PT from the disordered super-paraelectric to a cluster glass state. Similarities with superspin glass characteristics (Jönsson, 2004) are obvious. As described elsewhere (Shvartsman et al., 2008, Kleemann et al., 2009) this freezing process initiates also the transition of the Mn2+ spin moments into a spin glass state, which is magneto-electrically coupled to the dipolar glass state.

Another striking indicator of the dipolar glass state is the memory effect, which arises after isothermally annealing the sample below Tg. Fig. 12b shows an example of ‛burning a hole’ with depth Δχ (Tw) = χwait – χref(Tw) -6 at the wait temperature Tw = 32.5 K after waiting for tw 10.5 h (arrow). It signifies the asymptotic approach to the ‛stiff’ glassy ground state at Tw, which has the lowest susceptibility with respect to an external homogeneous field. Since the structure of the glassy ground state varies as a function of the temperature, the system is ‛rejuvenating’ at temperatures sufficiently far from Tw (Jönsson, 2004), hence, localizing the ‛burnt hole’ around Tw. This contrasts with the global decrease expected for an ordinarily relaxing metastable system.

Figure 12.

a) χ(T)of Sr0.98Mn0.02TiO3 recorded at Eac= 60 V/m and frequencies f = 10-1, 100, 101, 102, 103, 104, 105, and 0.4106 Hz. (b) Frequency (f) dependence of the peak temperature (Tm) of (T) taken from (a), plotted as ln(f/Hz) vs. Tm, and fitted by a critical power law (solid line), and difference curve Δχ = χwait – χref vs. T obtained at f = 10 Hz and Eac= 60 V/m upon heating after zero-field cooling from 80 K and waiting for 10.5 h at Tw = 32.5 K (wait) or without waiting (ref), respectively. Tg and Tw are marked by arrows. (c) χ(f)as measured at frequencies 10-3 f 106 Hz and temperatures 28.1 T 47.7 K in 2.2 K steps.

Figure 13.

Temperature dependence of the third-order nonlinear dielectric susceptibility of Sr0.98Mn0.02TiO3 measured at f = 37, 111, 333 and 1000 Hz. The dipolar glass freezing temperature Tg ≈ 38 K is indicated by a vertical arrow.

The glass transition may also be judged from spectra χvs. log f (Shvartsman et al., 2008) similarly as shown for K0.989Li0.011TaO3 (Fig. 10c). Since χ(f) measures the distribution function of relaxation times, its extension over more than nine decades of frequencies clearly signifies the glassy nature of the system. At low frequencies, f < 10-1 Hz, the low-f branch of χ(T) is observed to gradually lift up and to become horizontal at T < 38.9 K. This suggests it to extend with finite amplitude to fmin 0, hence, τmax →∞ for the percolating glass cluster. Maybe the ultimate proof of the very existence of a generic dipolar glass is given by the dynamic nonlinear susceptibility χ3=(3P/E3)/(6ε0) (Miga et al., 2007) in Fig. 13. This seems to indicate a divergence as T Tg and f 0 as predicted by theory (Binder & Reger, 1992), but has rarely been evidenced on orientational glasses (Hemberger et al., 1996). At the lowest frequency, f = 37 Hz, a fairly sharp peak is encountered, whose high-T branch might be considered as a critical hyperbola. Excess noise, however, prevents from seriously fitting a critical exponent, which should be close to = 1 as found previously (Hemberger et al., 1996). However, since the nonlinear susceptibility is modified by non-diverging ferroic (viz. ferroelectric) correlations due to progressive cluster formation on cooling, we cannot expect rigorous proportionality to the ‛spin’ glass susceptibility (Binder & Reger, 1992). Very probably the expected divergence is damped out similarly as in the case of KLT (Fig. 11).


5. Conclusion

In subsections 4.1 – 4.4 we presented different groups of polar materials separately and compared qualitatively their dielectric properties with predictions of suitable theories. It is likewise interesting to compare quantitatively different materials. The scaled nonlinear susceptibility, a3, is defined independently of the kind of material and its symmetry. This quantity is a measure of the nonlinearity of the investigated object. For an adequate comparison we have chosen values of a3 within centrosymmetric phases of different materials close to their temperatures of phase transition or of peak positions of the linear susceptibility, respectively. Fig. 14 shows thus collected values of a3. In view of the very large differences between the different a3 values a logarithmic scale is used. Consequently, only the magnitude values, |a3|, are presented. The highest |a3| was found for Rochelle salt, RS. Lower nonlinearity appears in sequence in the ferroelectric crystals TGS, LGO, BT, barium doped LGO, multiglass Sr0.98Mn0.02TiO3, orientational glass K0.989Li0.011TaO3 and relaxor ferroelectrics PMN and SBN 61. Obviously, classic ferroelectrics undergoing continuous PT are characterized by high nonlinearity, while structural disorder and presence of PNRs diminish the nonlinearity. Consequently relaxor ferroeletrics are characterized by the smallest values of a3. In other words, displacive ferroelectrics exhibiting soft-mode softening are most affected by nonlinearity, while typical order-disorder systems do not obtain their ferroelectricity primarily from the nonlinear interionic potential. The mechanism of their phase transition rather reflects the statistics of local hopping modes, in particular when being accompanied by quenched random fields as in relaxor ferroelectrics (Westphal et al., 1992; Kleemann et al., 2002), but probably also in partial order-disorder systems like BaTiO3 (Zalar et al., 2003).

Figure 14.

Absolute values of the scaled nonlinear susceptibility, |a3|, within centrosymmetric phases close to temperatures of phase transitions (for classic ferroelectrics) or peak positions of the linear susceptibility at low-f (for relaxor ferroelectrics and glasses).

The comparison of nonlinear dielectric response of various kinds of polar materials presented in this chapter gives evidence that such kind of measurement is a very sensitive tool for determination of the nature of ferroelectrics. Particularly useful for this purpose are the third-order dielectric susceptibility and the scaled non-linear susceptibility, a3. In our opinion, the second-order dielectric susceptibility is less significant, since it is more characteristic of the sample state than of a particular group of ferroelectrics. If anything, this susceptibility contains information about the distinct polar state of the sample. This information may be used for checking the presence of a center of inversion. The nonlinear dielectric response of classic ferroelectric crystals displaying continuous or discontinuous phase transition stays in a good agreement with predictions of the thermodynamic theory of ferroelectric phase transition. Predictions of the scaling theory for TGS crystal are also successfully verified experimentally. The situation is much more complex in disordered systems like ferroelectric relaxors or dipolar glasses. Respective theories properly explaining the observed features have still to be developed and tested via dynamic nonlinear dielectric response.



The authors are grateful to D. Rytz, A. Tkach, and P.M. Vilarinho for providing samples. WK thanks the Foundation for Polish Science (FNP), Warsaw, for an Alexander von Humboldt Honorary research grant.


  1. 1. BinderK.RegerJ. D. 1992 Theory of orientational glasses: models, concepts, simulations, Adv. Phys., 41 547627
  2. 2. BinderK.YoungA. P. 1986 Spin glasses: Experimental facts, theoretical concepts, and open questions, Rev. Mod. Phys., 58 801976
  3. 3. BobnarV.KutnjakZ.PircR.BlincR.LevstikA. 2000 Crossover from glassy to inhomogeneous-ferroelectric nonlinear dielectric response in relaxor ferroelectrics, Phys. Rev. Lett., 84 58925895
  4. 4. BonneauP.GarnierP.CalvarinG.HussonE.GavarriJ. R.HewatA. W.MorellA. 1991 X-ray and neutron diffraction studies of the diffuse phase transition in PbMg1/3Nb2/3O3 ceramics, J. Solid State Chem., 91 350361
  5. 5. BöttcherC. J. F. 1973 Theory of Electric Polarization, Elsevier, Amsterdam, Oxford, New York
  6. 6. BurnsG.DacolF. H. 1983 Crystalline ferroelectrics with glassy polarization behavior, Phys. Rev. B, 28 25272530
  7. 7. CachR. 1988 Time changes of the internal bias field in -damaged TGS crystal, Acta Univ. Wrat. LIII, 1084 125133
  8. 8. ChenA.ZhiY. 2004 DC electric-field dependence of the dielectric constant in polar dielectrics: “multi-polarization-mechanism” model, Phys. Rev. B, 69 174109174101 - 174109-8
  9. 9. DebyeP. 1929 Polar Molecules, The Chemical Catalog Company, Inc., New York
  10. 10. DecJ.MigaS.KleemannW.DkhilB. 2008 Nonlinear dielectric properties of PMN relaxor crystals within Ginzburg-Landau-Devonshire approximation, Ferroelectrics, 363 141149
  11. 11. DecJ.MigaS.TrybułaZ.KaszyńskaK.KleemannW. 2010 Dynamics of Li+ dipoles at very low concentration in quantum paraelectric potassium tantalate, J. Appl. Phys., 107 094101 - 094102-8
  12. 12. DevonshireA. F. 1949 Theory of barium titanate, Part I, Phil. Mag., 40 10401063
  13. 13. FisherD. S. 1986 Scaling and critical slowing down in random-field Ising systems, Phys. Rev. Lett., 56 416419
  14. 14. FujimotoM. 2003 The physics of structural phase transitions (2nd ed.), Springer-Verlag, Berlin, Heidelberg, New York
  15. 15. GinzburgV. L. 1945 On the dielectric properties of ferroelectric (Seignetteelectric) crystals and barium titanate, Zh. Exp. Theor. Phys., 15 739749
  16. 16. GlazounovA. E.TagantsevA. K. 2000 Phenomenological model of dynamic nonlinear response of relaxor ferroelectrics, Phys. Rev. Lett. 85 21922195
  17. 17. HembergerJ.RiesH.LoidlA.BöhmerR. 1996 Static freezing transition at a finite temperature in a quasi-one-dimensional deuteron glass, Phys. Rev. Lett., 76 23302333
  18. 18. HöchliU. T. 1982 Dynamics of freezing electric dipoles, Phys. Rev. Lett., 48 14941497
  19. 19. HöchliU. T.KnorrK.LoidlA. 1990 Orientational glasses, Advan. Phys., 39 405615
  20. 20. IkedaS.KominamiH.KoyamaK.WadaY. 1987 Nonlinear dielectric constant and ferroelectric-to-paraelectric phase transition in copolymers of vinylidene fluoride and trifluoroethylene, J. Appl. Phys., 62 33393342
  21. 21. JönssonP. E. 2004 Superparamagnetism and spin-glass dynamics of interacting magnetic nanoparticle systems, Adv. Chem. Phys., 128 191248
  22. 22. JonscherA. K. 1983 Dielectric Relaxation in Solids, Chelsea Dielectrics, London
  23. 23. KleemannW.KützS.RytzD. 1987 Cluster glass and domain state properties of KTaO3:Li, Europhys. Lett., 4 239245
  24. 24. KleemannW.DecJ.LehnenP.BlincR.ZalarB.PankrathR. 2002 Uniaxial relaxor ferroelectrics: the ferroic random-field Ising model materialized at last, Europhys. Lett., 57 1419
  25. 25. KleemannW.ShvartsmanV. V.BedantaS.BorisovP.TkachA.VilarinhoP. M. 2008 (Sr,Mn)TiO3- a magnetoelectrically coupled multiglass, J. Phys.: Condens. Matter, 20 434216434211 - 434216-6
  26. 26. KleemannW.BedantaS.BorisovP.ShvartsmanV. V.MigaS.DecJ.TkachA.VilarinhoP. M. 2009 Multiglass order and magnetoelectricity in Mn2+ doped incipient ferroelectrics, Eur. Phys. J. B, 71 407410
  27. 27. KleemannW.DecJ.MigaS.RytzD. 2011 Polar states of the impurity system KTaO3:Li, Z. Kristallogr. 226 145149
  28. 28. KremerF.SchönhalsA. .Eds. 2003 Broadband Dielectric Spectroscopy, Springer-Verlag, Berlin, Heidelberg, New York
  29. 29. LagutaV. V.KondakovaI. V.BykovI. P.GlinchukM. D.TkachA.VilarinhoP. M.JastrabikL. 2007 Electron spin resonance investigation of Mn2+ ions and their dynamics in Mn-doped SrTiO3, Phys. Rev. B, 76 054104054101 - 054104-6
  30. 30. LebedevA. I.SluchinskajaI. A.ErkoA.KozlovskiiA. F. 2009 Direct evidence for off-centering of Mn impurity in SrTiO3, JETP-Lett., 89 457467
  31. 31. Leont’evI. N.LeidermanA.TopolovV.YuFesenkoO. E. 2003 Nonlinear properties of barium titanate in the electric field range 0 ≤ E ≤ 5.5×107 V/m, Phys. Solid State, 45 11281130
  32. 32. LevinI.KrayzmanV.WoicikJ. C.TkachA.VilarinhoP. M. 2010 X-ray absorption fine structure studies of Mn coordination in doped perovskite SrTiO3, Appl. Phys. Lett., 96 052901 - 052904-3
  33. 33. LinesM. E.GlassA. M. 1977 Principles and Applications of Ferroelectrics and Related Materials, Oxford University Press, London
  34. 34. MierzwaW.FugielB.ĆwikielK. 1998 The equation-of-state of triglycine sulphate (TGS) ferroelectric for both phases near the critical point, J. Phys.: Condens. Matter, 10 88818892
  35. 35. MigaS.DecJ.MolakA.KoralewskiM. 2006 Temperature dependence of nonlinear susceptibilities near ferroelectric phase transition of a lead germanate single crystal, J. Appl. Phys., 99 124107124101 - 124107-6
  36. 36. MigaS.DecJ.KleemannW. 2007 Computer-controlled susceptometer for investigating the linear and non-linear dielectric response, Rev. Sci. Instrum., 78 033901 - 033902-7
  37. 37. MigaS.DecJ. 2008 Non-linear dielectric response of ferroelectric and relaxor materials, Ferroelectrics, 367 223228
  38. 38. MigaS.DecJ.MolakA.KoralewskiM. 2008 Barium doping-induced polar nanore-gions in lead germanate single crystal, Phase Trans., 81 11331140
  39. 39. MigaS.CzaplaZ.KleemannW.DecJ. 2010a Non-linear dielectric response in the vicinity of the ‘inverse melting’ point of Rochelle salt, Ferroelectrics, 400 7680
  40. 40. MigaS.KleemannW.DecJ. 2010b Non-linear dielectric susceptibility near to the field-induced ferroelectric phase transition of K0.937Li0.063TaO3, Ferroelectrics, 400 3540
  41. 41. MitsuiT. 1958 Theory of the ferroelectric effect in Rochelle salt, Phys. Rev., 111 12591267
  42. 42. OliverJ. R.NeurgaonkarR. R.CrossL. E. 1988 A thermodynamic phenomenology for ferroelectric tungsten bronze Sr0.6Ba0.4Nb2O6 (SBN:60), J. Appl. Phys., 64 3747
  43. 43. PircR.TadićB.BlincR. 1994 Nonlinear susceptibility of orientational glasses, Physica B, 193 109115
  44. 44. PircR.BlincR. 1999 Spherical random-bond-random-field model of relaxor ferroelectrics, Phys. Rev. B, 60 1347013478
  45. 45. ShvartsmanV. V.BedantaS.BorisovP.KleemannW.TkachA.VilarinhoP. 2008 (Sr,Mn)TiO3- a magnetoelectric multiglass, Phys. Rev. Lett., 101 165701 - 165704-4
  46. 46. SmolenskiiG. A.IsupovV. A.AgranovskayaA. I.PopovS. N. 1960 Ferroelectrics with diffuse phase transition (in Russian), Sov. Phys.-Solid State, 2 29062918
  47. 47. StanleyH. E. 1971 Introduction to Phase Transitions and Critical Phenomena, Clarendon, Oxford
  48. 48. TkachA.VilarinhoP. M.KholkinA. L. 2005 Structure-microstructure-dielectric tunability relationship in Mn-doped strontium titanate ceramics, Acta Mater., 53 50615069
  49. 49. TkachA.VilarinhoP. M.KholkinA. L. 2006 Dependence of dielectric properties of manganese-doped strontium titanate ceramics on sintering atmosphere, Acta Mater., 54 53855391
  50. 50. TkachA.VilarinhoP. M.KholkinA. L. 2007 Non-linear dc electric-field dependence of the dielectric permittivity and cluster polarization of Sr1−xMnxTiO3 ceramics, J. Appl. Phys., 101 084111 - 084110-9
  51. 51. ValasekJ. 1920 Piezoelectric and allied phenomena in Rochelle salt, Phys. Rev., 15 537538
  52. 52. ValasekJ. 1921 Piezo-electric and allied phenomena in Rochelle salt, Phys. Rev., 17 475481
  53. 53. vonHippel. A. 1950 Ferroelectricity, domain structure, and phase transitions of barium titanate, Rev. Mod. Phys., 22 222237
  54. 54. vonHippel. A. 1954 Dielectrics and Waves, Wiley, New York
  55. 55. VugmeisterB. E.GlinchukM. D. 1990 Dipole glass and ferroelectricity in random-site electric dipole systems, Rev. Mod. Phys., 62 9931026
  56. 56. WangY. L.TagantsevA. K.DamjanovicD.SetterN. 2006 Anharmonicity of BaTiO3 single crystals, Phys. Rev. B, 73 132103132101 - 132103-4
  57. 57. WangY. L.TagantsevA. K.DamjanovicD.SetterN.YarmarkinV. K.SokolovA. I.LukyanchukI. A. 2007 Landau thermodynamic potential for BaTiO3, J. Appl. Phys., 101 104111 - 104115-9
  58. 58. WeiX.YaoX. 2006a Reversible dielectric nonlinearity and mechanism of electrical tunability for ferroelectric ceramics, Int. J. Mod. Phys. B, 20 29772998
  59. 59. WeiX.YaoX. 2006b Analysis on dielectric response of polar nanoregions in paraelectric phase of relaxor ferroelectrics, J. Appl. Phys., 100 064311 - 064319-6
  60. 60. WestphalV.KleemannW.GlinchukM. 1992 Diffuse phase transitions and random field-induced domain states of the “relaxor” ferroelectric PbMg1/3Nb2/3O3, Phys. Rev. Lett., 68 847950
  61. 61. WickenhöferF.KleemannW.RytzD. 1991 Dipolar freezing of glassy K1-xLixTaO3, x = 0.011, Ferroelectrics, 124 237242
  62. 62. ZalarB.LagutaV. V.BlincR. 2003 NMR evidence for the coexistence of order-disorder and displacive components in barium titanate, Phys. Rev. Lett., 90 037601037601 - 037601-4

Written By

Seweryn Miga, Jan Dec and Wolfgang Kleemann

Submitted: 08 November 2010 Published: 23 August 2011