Open access peer-reviewed chapter

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

Written By

Seweryn Miga, Jan Dec and Wolfgang Kleemann

Submitted: November 8th, 2010 Reviewed: March 19th, 2011 Published: August 23rd, 2011

DOI: 10.5772/20380

Chapter metrics overview

4,056 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 acsusceptometer for simultaneous measurements of the phase resolved complex linear and complex non-linear acsusceptibilities of lossy and dispersive dielectric materials (Miga, Dec & Kleemann, 2007). It allows measurements over a wide range of experimental variables, such as acamplitudes 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 Gof 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, Aand T0 are constants. Band Care 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 Bdetermines the kind of the PT. For positive Ba continuous (second-order) PT occurs. On the other hand, for negative Ba 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 Eand Pdenotes the electric equation of state. After successive differentiation of this equation with respect to Pone 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. Additionally2 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 Psvanishes 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=PsandPS2=(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 Jijare glass-like random intercluster couplings or bonds, hiand Erandom 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 jione can calculate the linear and nonlinear dielectric susceptibilities ias 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 Tcthe 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 Tcis 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 Tca 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 Tf220 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)>0is the reduced temperature and zν=6.6±0.2is the dynamical critical exponent. Essentially the same (static) glass temperature, Tg= (33.5±1.5) K, emerges for the ‛largest’ relaxation time τmas 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.Tfor 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 χ3drops 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 Tmin 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 tw10.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. Tgand Tware 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 Tgand 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.1992Theory of orientational glasses: models, concepts, simulations,Adv. Phys.,41547627
  2. 2. BinderK.YoungA. P.1986Spin glasses: Experimental facts, theoretical concepts, and open questions,Rev. Mod. Phys.,58801976
  3. 3. BobnarV.KutnjakZ.PircR.BlincR.LevstikA.2000Crossover from glassy to inhomogeneous-ferroelectric nonlinear dielectric response in relaxor ferroelectrics,Phys. Rev. Lett.,8458925895
  4. 4. BonneauP.GarnierP.CalvarinG.HussonE.GavarriJ. R.HewatA. W.MorellA.1991X-ray and neutron diffraction studies of the diffuse phase transition in PbMg1/3Nb2/3O3 ceramics,J. Solid State Chem.,91350361
  5. 5. BöttcherC. J. F.1973Theory of Electric Polarization, Elsevier, Amsterdam, Oxford, New York
  6. 6. BurnsG.DacolF. H.1983Crystalline ferroelectrics with glassy polarization behavior,Phys. Rev. B,2825272530
  7. 7. CachR.1988Time changes of the internal bias field in -damaged TGS crystal,Acta Univ. Wrat.LIII,1084125133
  8. 8. ChenA.ZhiY.2004DC electric-field dependence of the dielectric constant inpolar dielectrics: “multi-polarization-mechanism” model,Phys. Rev. B,69174109174101- 174109-8
  9. 9. DebyeP.1929Polar Molecules,The Chemical Catalog Company, Inc., New York
  10. 10. DecJ.MigaS.KleemannW.DkhilB.2008Nonlinear dielectric properties of PMN relaxor crystals within Ginzburg-Landau-Devonshire approximation,Ferroelectrics,363141149
  11. 11. DecJ.MigaS.TrybułaZ.KaszyńskaK.KleemannW.2010Dynamics of Li+ dipoles at very low concentration in quantum paraelectric potassium tantalate,J. Appl. Phys.,107094101 - 094102-8
  12. 12. DevonshireA. F.1949Theory of barium titanate,Part I,Phil. Mag.,4010401063
  13. 13. FisherD. S.1986Scaling and critical slowing down in random-field Ising systems,Phys. Rev. Lett.,56416419
  14. 14. FujimotoM.2003The physics of structural phase transitions(2nd ed.), Springer-Verlag, Berlin, Heidelberg, New York
  15. 15. GinzburgV. L.1945On the dielectric properties of ferroelectric (Seignetteelectric) crystals and barium titanate,Zh. Exp. Theor. Phys.,15739749
  16. 16. GlazounovA. E.TagantsevA. K.2000Phenomenological model of dynamic nonlinear response of relaxor ferroelectrics,Phys. Rev. Lett.8521922195
  17. 17. HembergerJ.RiesH.LoidlA.BöhmerR.1996Static freezing transition at a finite temperature in a quasi-one-dimensional deuteron glass,Phys. Rev. Lett.,7623302333
  18. 18. HöchliU. T.1982Dynamics of freezing electric dipoles,Phys. Rev. Lett.,4814941497
  19. 19. HöchliU. T.KnorrK.LoidlA.1990Orientational glasses,Advan. Phys.,39405615
  20. 20. IkedaS.KominamiH.KoyamaK.WadaY.1987Nonlinear dielectric constant and ferroelectric-to-paraelectric phase transition in copolymers of vinylidene fluoride and trifluoroethylene,J. Appl. Phys.,6233393342
  21. 21. JönssonP. E.2004Superparamagnetism and spin-glass dynamics of interacting magnetic nanoparticle systems,Adv. Chem. Phys.,128191248
  22. 22. JonscherA. K.1983Dielectric Relaxation in Solids, Chelsea Dielectrics, London
  23. 23. KleemannW.KützS.RytzD.1987Cluster glass and domain state properties of KTaO3:Li,Europhys. Lett.,4239245
  24. 24. KleemannW.DecJ.LehnenP.BlincR.ZalarB.PankrathR.2002Uniaxial relaxor ferroelectrics: the ferroic random-field Ising model materialized at last,Europhys. Lett.,571419
  25. 25. KleemannW.ShvartsmanV. V.BedantaS.BorisovP.TkachA.VilarinhoP. M.2008(Sr,Mn)TiO3- a magnetoelectrically coupled multiglass,J. Phys.: Condens. Matter,20434216434211- 434216-6
  26. 26. KleemannW.BedantaS.BorisovP.ShvartsmanV. V.MigaS.DecJ.TkachA.VilarinhoP. M.2009Multiglass order and magnetoelectricity in Mn2+ doped incipient ferroelectrics, Eur. Phys. J. B,71407410
  27. 27. KleemannW.DecJ.MigaS.RytzD.2011Polar states of the impurity system KTaO3:Li,Z. Kristallogr.226145149
  28. 28. KremerF.SchönhalsA. .Eds.2003Broadband Dielectric Spectroscopy, Springer-Verlag, Berlin, Heidelberg, New York
  29. 29. LagutaV. V.KondakovaI. V.BykovI. P.GlinchukM. D.TkachA.VilarinhoP. M.JastrabikL.2007Electron spin resonance investigation of Mn2+ ions and their dynamics in Mn-doped SrTiO3,Phys. Rev. B,76054104054101- 054104-6
  30. 30. LebedevA. I.SluchinskajaI. A.ErkoA.KozlovskiiA. F.2009Direct evidence for off-centering of Mn impurity in SrTiO3,JETP-Lett.,89457467
  31. 31. Leont’evI. N.LeidermanA.TopolovV.YuFesenkoO. E.2003Nonlinear properties of barium titanate in the electric field range 0 ≤ E ≤ 5.5×107 V/m,Phys. Solid State,4511281130
  32. 32. LevinI.KrayzmanV.WoicikJ. C.TkachA.VilarinhoP. M.2010X-ray absorption fine structure studies of Mn coordination in doped perovskite SrTiO3,Appl. Phys. Lett.,96052901 - 052904-3
  33. 33. LinesM. E.GlassA. M.1977Principles and Applications of Ferroelectrics and Related Materials, Oxford University Press, London
  34. 34. MierzwaW.FugielB.ĆwikielK.1998The equation-of-state of triglycine sulphate (TGS) ferroelectric for both phases near the critical point,J. Phys.: Condens. Matter,1088818892
  35. 35. MigaS.DecJ.MolakA.KoralewskiM.2006Temperature dependence of nonlinear susceptibilities near ferroelectric phase transition of a lead germanate single crystal,J. Appl. Phys.,99124107124101- 124107-6
  36. 36. MigaS.DecJ.KleemannW.2007Computer-controlled susceptometer for investigating the linear and non-linear dielectric response,Rev. Sci. Instrum.,78033901 - 033902-7
  37. 37. MigaS.DecJ.2008Non-linear dielectric response of ferroelectric and relaxor materials,Ferroelectrics,367223228
  38. 38. MigaS.DecJ.MolakA.KoralewskiM.2008Barium doping-induced polar nanore-gions in lead germanate single crystal,Phase Trans.,8111331140
  39. 39. MigaS.CzaplaZ.KleemannW.DecJ.2010aNon-linear dielectric response in the vicinity of the ‘inverse melting’ point of Rochelle salt,Ferroelectrics,4007680
  40. 40. MigaS.KleemannW.DecJ.2010bNon-linear dielectric susceptibility near to the field-induced ferroelectric phase transition of K0.937Li0.063TaO3,Ferroelectrics,4003540
  41. 41. MitsuiT.1958Theory of the ferroelectric effect in Rochelle salt,Phys. Rev.,11112591267
  42. 42. OliverJ. R.NeurgaonkarR. R.CrossL. E.1988A thermodynamic phenomenology for ferroelectric tungsten bronze Sr0.6Ba0.4Nb2O6 (SBN:60),J. Appl. Phys.,643747
  43. 43. PircR.TadićB.BlincR.1994Nonlinear susceptibility of orientational glasses,Physica B,193109115
  44. 44. PircR.BlincR.1999Spherical random-bond-random-field model of relaxor ferroelectrics,Phys. Rev. B,601347013478
  45. 45. ShvartsmanV. V.BedantaS.BorisovP.KleemannW.TkachA.VilarinhoP.2008(Sr,Mn)TiO3- a magnetoelectric multiglass,Phys. Rev. Lett.,101165701 - 165704-4
  46. 46. SmolenskiiG. A.IsupovV. A.AgranovskayaA. I.PopovS. N.1960Ferroelectrics with diffuse phase transition (in Russian),Sov. Phys.-Solid State,229062918
  47. 47. StanleyH. E.1971Introduction to Phase Transitions and Critical Phenomena, Clarendon, Oxford
  48. 48. TkachA.VilarinhoP. M.KholkinA. L.2005Structure-microstructure-dielectric tunability relationship in Mn-doped strontium titanate ceramics,Acta Mater.,5350615069
  49. 49. TkachA.VilarinhoP. M.KholkinA. L.2006Dependence of dielectric properties of manganese-doped strontium titanate ceramics on sintering atmosphere,Acta Mater.,5453855391
  50. 50. TkachA.VilarinhoP. M.KholkinA. L.2007Non-linear dc electric-field dependence of the dielectric permittivity and cluster polarization of Sr1−xMnxTiO3 ceramics,J. Appl. Phys.,101084111 - 084110-9
  51. 51. ValasekJ.1920Piezoelectric and allied phenomena in Rochelle salt,Phys. Rev.,15537538
  52. 52. ValasekJ.1921Piezo-electric and allied phenomena in Rochelle salt,Phys. Rev.,17475481
  53. 53. vonHippel. A.1950Ferroelectricity, domain structure, and phase transitions of barium titanate,Rev. Mod. Phys.,22222237
  54. 54. vonHippel. A.1954Dielectrics and Waves, Wiley, New York
  55. 55. VugmeisterB. E.GlinchukM. D.1990Dipole glass and ferroelectricity in random-site electric dipole systems,Rev. Mod. Phys.,629931026
  56. 56. WangY. L.TagantsevA. K.DamjanovicD.SetterN.2006Anharmonicity of BaTiO3 single crystals,Phys. Rev. B,73132103132101- 132103-4
  57. 57. WangY. L.TagantsevA. K.DamjanovicD.SetterN.YarmarkinV. K.SokolovA. I.LukyanchukI. A.2007Landau thermodynamic potential for BaTiO3,J. Appl. Phys.,101104111 - 104115-9
  58. 58. WeiX.YaoX.2006aReversible dielectric nonlinearity and mechanism of electrical tunability for ferroelectric ceramics,Int. J. Mod. Phys. B,2029772998
  59. 59. WeiX.YaoX.2006bAnalysis on dielectric response of polar nanoregions in paraelectric phase of relaxor ferroelectrics,J. Appl. Phys.,100064311 - 064319-6
  60. 60. WestphalV.KleemannW.GlinchukM.1992Diffuse phase transitions and random field-induced domain states of the “relaxor” ferroelectric PbMg1/3Nb2/3O3,Phys. Rev. Lett.,68847950
  61. 61. WickenhöferF.KleemannW.RytzD.1991Dipolar freezing of glassy K1-xLixTaO3, x = 0.011,Ferroelectrics,124237242
  62. 62. ZalarB.LagutaV. V.BlincR.2003NMR evidence for the coexistence of order-disorder and displacive components in barium titanate,Phys. Rev. Lett.,90037601037601- 037601-4

Written By

Seweryn Miga, Jan Dec and Wolfgang Kleemann

Submitted: November 8th, 2010 Reviewed: March 19th, 2011 Published: August 23rd, 2011