Fractal Phenomena in Ferroelectrics

This paper discusses on ferroelectricity from fractal analysis viewpoints. Micro to macro ferroelectric domain growth changes the domain dynamics, leading to the changes in various physical properties such as dielectric permittivity, loss, polarization hysteresis, light scattering and acoustic emission. *Corresponding author: Uchino, K. International Ctr. for Actuators & Transducers, The Pennsylvania State University, University Park, PA 16802, USA, E-mail: kenjiuchino@psu.edu


Introduction
The concept of "fractal" means "self-similarity", that is, the contour/structure of a part of an object is similar to that of the whole object. In other words, the contour/structure seems to be similar, irrelevant to the observing scale. Though the fractal exhibits a rather complicated structure in general, it is expressed quantitatively in a simple way with using the "fractal dimension". Figure 1 shows various growth shapes of the snow crystal for different size ranges [1] . It is noteworthy that apparently-similar 6-fold crystal growth shapeswith a tree-branch pattern is observed irrelevant to the crystal size from μm to mm.The fractal dimension(defined later) of the snowflake is reported about 1.3, much less than the Euclidean dimenstion 2. This sort of geometrical, structural and dimensional specific may provide characteristic physical properties. We will discuss ferroelectricity from a fractal analysis viewpoint in this paper, based on our previously reported experimental papers: • Critical Exponents of the Dielectric Constantsin Diffused-Phase-Transition Crystals: Kenji Uchino and Shoichi-crack propagation, 2. Convolutional effect in non-linear mechanical system, such as chaos, 3. Critical phenomenon in an equilibrium condition, e.g., phase transition, with which this paper is most closely related.

Intuitive origin of fractal
Let us consider the origin of self-similarity of snowflake shapes shown in Figure 1. Referring to Figure 2, when a new cluster (i.e., water particle) approaches to an existing snowflake, it can be deposited easily at the tip of the branch, but may be difficult to reach into the joint part. This intuitive model can explain that the snow crystal grows with keeping the original 6-fold shape. The key is to consider the cluster (new comer)-cluster (ready-existing) interaction. Because ferroelectrics (relaxor ferroelectrics, in particular) generate a macroscopic domain pattern with decreasing temperature through its phase transition (i.e., Curie) temperature originated from micro-micro domain interaction from an elevated temperature, similar fractal behavior can be anticipated. We will discuss the dielectric relaxation and loss in terms of the micro-micro domain interaction in Section IV.

Fractal dimension and scaling
Fractal dimension is defined as follows: There are N particles with each radius ri, which make a cluster. Average cluster radius R can be estimated by If the number of particles N is related with R as D is called "fractal dimension". Let us consider an example pictured in Figure3. We consider a small square particle in a 2-D cluster with a radius R (Euclidean dimension=2). If there is no particle-particle interaction and a circle is fully filled (i.e., Tiling) Figure3(a), since Area=π•R 2 , its fractal dimension is 2. However, like a snowflake, if particles will not reach to the branch joint part Figure3(b), the number N is reduced, leading to a smaller D than 2 (1<D<2).

Figure4:
Characteristic length L of a Koch's snowflake created by an iterative process of equilateral triangles. Self-similarity is maintained from (a), through (b) to (c). Koch fractal dimension D=1. 26. Another Koch curve example is pictured in Figure4. This snowflake-like Figure is generated based on an equilateral triangle. ⅓-length equilateral triangles are aligned iteratively on the large triangle edges. Self-similarity can be maintained with reducing the scale from (a), through (b) to (c). A rough diameter of the Figure L is called "characteristic length". If the contour length P is related with L as D is another fractal dimension. The D for the Koch curve is known to be D=1.26 [9] . You can imagine that if there is a specific self-similarity in the shape/contour, the fractal dimension D may be smaller than the Euclidean dimension d. Though knowing the fractal dimension D does not provide a concrete idea on the contour shape, it will provide at least the suggestion that there may be a sort of cluster-cluster interaction and that we may expect a self-similarity in the domain contour shape.
Supposing that a contour shape is a function f(r) and scaling parameter is λ, this function should keep This is known as "scaling law". By treating Eq. (4) as We obtain an important relationship Leading to the power law for the function g (λ): Here, p can be non-integer.

Domain and Dielectric Properties
The relationship between the domain structure and its dielectric property is considered in this Section. First, we exhibit domain dynamic reorientation mechanisms under electric field observed for typical normal ferroelectric BaTiO 3 (BT) at room temperature Figure 5(a) and relax-orferroelectricPb (Zn 1/3 Nb 2/3 )O 3 (PZN) at 100°C Figure. 5(b) [10] . The electric field was applied along the spontaneous polarization direction(from the bottom to top of the pictures); that is, [100] and [111] of the perovskite axes, respectively in BT and PZN. Typical 45° angled domain walls in BT with respect to the [100] perovskite orientation move perpendicularly to the domain planes independently each other (i.e., tetragonal 90°-domain reorientation), and disappear almost completely with increasing the eternal field much higher than the coercive field. That is, the domain walls move 1-dimensionally. To the contrary, in PZN, "spindle-like" domains (i.e., rhombohedral non-180° (or 73°)-domain walls become slightly "fat" without changing the length much with increasing the external field. The color reversal between black and white corresponds to the 180°-domain reversal (i.e., polarization reversal).Note that the domain walls do not disappear even under a large electric field, but just stay with breathing-like cyclical motion. That is, the domain walls move primarily 2-dimensionally. Figure 6: Dielectric constant vs. temperature curves and domain patternsaPb(Zn 1/3 Nb 2/3 )O 3 single crystal: (a) for an un-poled sample, and (b) for a poled sample [11] .

Domain and dielectric relaxation
More intriguing phenomena are observed in the relation of the dielectric property and domain structure. Refer to Figures 6(a) and 6(b), which show dielectric constant vs. temperature curves and domain patterns in the pure PZN single crystals at room temperature (RT): (a) for an un-poled sample, and (b) for a poled sample [11] . The domain picture at RT for the poled sample clearly exhibits relatively macroscopic domains (spindle-like with 3µm wide and 30~50 µm long), while for the un-poled (or annealed) sample, the picture does not exhibit visible domains, because the domain size seems to be too small (i.e., micro (ornano) domains). In the un-poled sample, significant dielectric "relaxation" (i.e., large frequency dependence of the permittivity) can be observed in the whole temperature range below the Curie region; while in the poled sample, both dielectric constant value and dielectric relaxation diminish significantly in the low temperature range (RTup to 105°C). Furthermore, relatively large dielectric loss tan δ in the un-poled sample reduces dramatically when the macro domains are observed. Just around 105°C, where the macroscopic domain pattern becomes invisible (i.e., micro (nano)-domain), both permittivity and dielectric loss jump up, associated also with dielectric relaxation. We can conclude thus: 1. The dielectric relaxation phenomenon and high dielectric loss tan δ in the un-poled PZN seem to be originated from the micro-domain dynamics. 2. There seems to exist a sort of transition, macro-to-micro domain change, around 105°C in PZN.

Critical Exponent in Relaxor
We discuss on the critical exponent γ in the relation between the dielectric constant and temperature (or hydrostatic pressure) in this section [2] . A high correlation of the γ value with the phase transitiondiffuseness has been found empirically, which is discussed based on the fractal viewpoint.
well-known Curie-Weiss law (i.e., ɛ -ɛ ∞ = C/(T -T 0 ), where ɛ∞is the dielectric constant measure dat an optical-range frequency) is observed in a wide temperature range around the transition point. This is demonstrated in Figure 7(a), which shows the reciprocal dielectric constant (1/ɛ -1/ɛ m ) versus temperature (T -T m ) relation for a Pb(Mg 1/3 Nb 2/3 )O 3 (PMN) single crystal (i.e., convex curve). The notations ɛ m and T m refer to the maximum dielectric constant and the corresponding temperature (i. e., average Curie temperature). Since the dielectric constant ɛ ∞ is considerably smaller than the static ɛ in the samples treated in this paper, it may be neglected in the formula.  [2] Smolenskii [12] and Kirillov et a1. [13] Pointed out that for ferroelectrics with a diffused phase transition (i.e., relaxor ferroelectrics), a quadratic relation between the dielectric constant and temperature (1/ɛ-1/ɛ m )= B(T -T m ) 2 gives better fitting to the experimental data. This "quadratic law" has been derived theoretically on the basis of a microscopic composition-fluctuation model [12,13] . Refer to Figure 8. When we assume that the statistical composition fluctuation in a complex perovskite produces microscopic regions each of which has similar dielectric characteristics, but a slightly different Curie temperature Figure. 8(a),and also that the distribution of the Curie temperature T C is of a Gaussian type: The reciprocal dielectric constant is represented by Where σ is a standard deviation of the Gaussian function [14,15] . It is notable that the Gaussian distribution supposed basically no strong interaction between microscopic regions/ clusters each other. The reciprocal permittivity of each micro-region is supposed to obey the Curie-Weiss law of the first-kind phase transition as follows: ( ) Where C and Δ are the Curie-Weiss constant and the difference between the Curie temperature T C and the Curie-Weiss temperature T 0 Figure.8 When σ is relatively large (i. e., very diffused), this theoretical treatment derives two approximate formulas for the extreme cases; T -T m << σ and T -T m >> σ. In the former case, the reciprocal permittivity is reduced into A series expansion, neglecting terms with (T -T m ) 4 and the higher, gives the quadratic law of On the other hand, in the case of T -T m >>σ, the temperature dependence of the permittivity approaches to the normal Curie-Weiss law.  [14] ; this suggests, strictly speaking, a kind of discrepancy of the above-described Smolenskii's theoretical model.

Crital exponent of permittivity
The discrepancy of the curve fitting in Figures 7(a) and 7(b) may be related with the Gaussian distribution by neglecting the micro domain-micro domain interaction. Actually, we can expect a fractal phenomenon due to the domain-domain interaction. On this background, we tried at first to establish some empirical rules which can represent the real phenomenon much better.
In this paper, we have determined the critical exponent γ precisely in a generalized form of Not only for relaxor ferroelectrics, but also for normal sharp-transition ferroelectrics, so that Eq. (12) may provide the best fit to the experimental results. Moreover, another critical exponent γ* of the permittivity with respect to hydrostatic pressure p, which is given by has also been determined for the same samples.
12PT) have been used for our dielectric measurements. Cut into arectangular plate with the polarization direction normal to thesurface, the crystal was silver-electroded, then electrically poled. The capacitance was precisely measured at 1 kHz and 1 MHz with an automatic capacitance bridge (Hewlett-Packard 4270A) up to about 200°C above the average transition temperature. The differenceof the permittivity between the poled and unpoled samples was hardly observed in the temperature range above the Curie point.
. Note the good linearity for all the samples. The γ values were determined in the temperature range of 5 °C < T -T m < 200°C. [2] Figure 9 shows logarithmic plots of the reciprocal permittivity (1/ɛ -1/ɛ m ) measured at 1 kHz as a function of temperature (T -T m ). The linearity of the curves is remarkably good except in a narrowtemperature range above the peak Tm (T -T m < l °C) for all the samples. The critical exponent γ of the reciprocal permittivity with respect to temperature has been determined from the slope of the straight line in the range of5°C<T -T m < 200°C: γ = 1.08 for BT, 1.17 for KTN, 1.58 for 0.88PZN-0.12PT, 1.64 for PMN and1.76 for PZN. Note that this γ is different from the conventional critical exponent defined for sharp-transition ferroe1ectrics such as BT in a very close vicinity of the Curie temperature (T -T C < 1°C).  [17] . Table 1. summarizes the critical exponent γ and the Curie-like constant C' for all the samples. It is important to mention here that the relaxor ferroelectric sreveal a large relaxation of the permittivity and also that the exponent γ might be changed with the measuring frequency. PMN, for example, shows a gradual decrease in the permittivity ɛ and a monotonous increase of the peak temperature T m with changing the measuring frequency from 1 kHz to 1 MHz, as is inserted in Figure10(a). Figure10(a) and 10(b) show the (1/ɛ -1/ɛ m )versus(T -T m ) plots for PMN and PZN, respectively, measured at 1 kHz and 1 MHz. The critical exponent γ seems to be insensitive to the frequency change at least in the range below 10 MHz, while the constant C' varies with the frequency. . [2] Pressure dependence of permittivity:The effect of hydrostatic pressure on the dielectric constantwas also investigated up to 6 kbar for PNN and KTN single crystals. Pressure was generated in an air-driven intensifier pump (Pressure Products Industries) and transmitted through the pressure medium of Plexol. Applied pressure was monitored by a Heise Gauge with the accuracy of ± 5 bar. Temperature was varied with a lamp heater and a CA thermocouple was set inside the pressure cell. Figure 11: Dependence of the permittivity on hydrostatic pressure for KTN (a) and for PMN (b). [2] Figure 11(a) and 11(b) show isotherms of the permittivity and its reciprocal for KTN [16] and PMN [14] , respectively, measured as a function of hydrostatic pressure. A good linear relation between the reciprocal permittivity and hydrostatic pressure (similar to the Curie-Weiss law) is obtained for KTN above a critical pressure p m which induces the ferroelectric-Para electric transition, whereas gradual increase of the slope is observed for PMN with increasing pressure (i.e., convex curve). Therefore, a new critical exponent γ* of the permittivity with respect to pressure can be introduced by the definition of Eq. (13). Figure 12 (a) and 12(b) show logarithmic plots of the reciprocal permittivity (1/ɛ -1/ɛ m )as a function of hydrostatic pressure (p -p m ). Since the critical pressure pm where the permittivity reveals the maximum could not be observed practically for PNN in our measurement above room temperature, we used a hypothetical p m (T) which is calculated for each observation temperature as Where ∂T m (p)/∂p is a constant of -4.8 [K•kbar -1 ] [14] . The slope gives the exponent γ* for each case: γ* = 1.00 for KTN and 1. 66 for PMN. The critical γ* exponents are also summarized in Table 1, including the value for BT reported by Samara [17] .

Critical Exponent Discussion
It is worth noting the correlation between the criticalexponent γ and the phase transition diffuseness. Refer to Table 1. BT revealing atypical sharp phase transition has a γ value nearly equal to 1: this corresponds to the Curie-Weiss law. The solid solution KTN withan "apparent" sharp transition reveals a slight deviation from theCurie-Weiss law. On the other hand, the γ value is as large as1.7 for PMN and PZN with a very diffused transition; however, thisis relatively smaller than expected for the "quadratic law" by Smolenskii [12] .The solid solution PZN-PT between a relaxor and a normal-transition ferroelectrics reveals an intermediate γ value of 1.58. Thecritical exponent γ of the permittivity, determined in the temperaturerange from several to several hundreds °C above the transition point,may at least represent the degree of the phase transition diffuseness.
It is also interesting that the critical exponent γ* of the permittivity with respect to hydrostatic pressure corresponds exactly to the γ value in both the relaxor and normal ferroelectrics. This correspondence is consistently explained, if we assume that the reciprocal dielectric constant is represented by the following general formula: The transformation of Eq. (15), And Eq. (14), in conjunction with Eq. (13) which is the definition of the critical exponent γ*.

Crital exponent and fractal dimension
The critical exponent around the phase transition has not been calculated theoretically yet except for special and simple models. We review here one of the simplest models, Ising model, for receiving preliminary and intuitive idea in the relationship between the critical exponent and fractal dimension.
Critical phenomenon in ising model: A magnet model for the spin-dimension D s = 1 (for Ising model) on the lattice dimension d = 2 (i.e., plane) is visualized in Figure  13. There are only two spin states (up or down for D s = 1), which are denoted σ X = ± 1, where x = (x, y, z, …) is the lattice point coordinate (only x and y for d = 2). We suppose that the spinspin coupling is isotripic in terms of the space dimension and the total system energy (Hamiltonian) is provided by where <x,y> means the summation for all combination of the adjacent lattice points.When Heisenberg Exchange Coupling: J is larger than zero, σ x =σ y = 1 or -1 (i. e., parallel alignment) provides more stable status, and σ x = 1,σ y = -1or σ x = -1,σ y = 1 gives higher energy status.On the other hand, when J< 0, σ x = 1,σ y = -1or σ x = -1,σ y = 1 gives more stable status (i. e., antiparallel alignment of the adjacent spins). When we consider the effect of an external magnetic field H, extended Hamiltonian is expressed by Summation of the second term is conducted for all lattice points, which promotes the spin alignment along the eternal magnetic field H.
When this spin system is set at a certain temperatureT atmosphere, the spin arrangement is continuously fluctuating due to the thermal fluctuation. In Boltzmann statistical thermodynamics, the probability of a certain status {σ X } under an absolute temperature T can be provided by Here, k is so-called Boltzmann constant, and Z Λ (T,H) is a normalization constant provided by Based on these assumption, spin-spin correlation function <σ x σ y > Λ (T,H) is expressed as When the temperature T is high enough, and the external field H = 0, the two-point correlation factor (regardless of lattice dimension d) is known to follow as <σ 0 σ x > (T,0) means what degree the spin σ 0 at the lattice "origin" point 0 (without losing generality) is aligned with the spin σ x at the lattice point x, and |x| is the distance between the origin and the lattice point x. ξ(T) is called "correlation distance", beyond which <σ 0 σ x > (T,0) becomes very small; in other words, the spins align generally in parallel inside ξ(T) region, but the spin alignment is not maintained beyond ξ(T). Remind the "characteristic distance" of fractal analysis in Section II, that corresponds to ξ(T). By skipping detailed discussion (refer to [8] ), in rough summary as depicted in Figure 13, randomly oriented spins at high temperature (T>T C ) are aligned blow a critical temperature T C (called Curie temperature), exhibiting spontaneous magnetization M(T) with "spontaneous symmetry break". With decreasing temperature, spins align more and more, and finally completely at T = 0 K. η-correlation function: β-Magnetization: Note that the first three exponents can be observed above the critical (Curie) temperature and on the last exponent β is for below T C .
All the critical exponents are not completely independent, but there is an important inter-relationship: If we consider the critical phenomenon at T = T C from the fractal analysis viewpoint, we can evaluate the spin-spin correlation function as Where D is the fractal dimension (not spin dimension Ds). Thus, by comparing Eq. (27) and Eq. (30), we obtain Leading to (from Eq. (29))  [18] . With increasing the spin dimension from 1 to ∞ (Sphere model), γ increases from 1.25 to 2. γ = 2 for the classic model seems to be rather different from the ferroelectrics (γ = 1 for classic ferroelectric phenomenology).
The fractal dimension D increases very slightly (or insensitive) with D s . To the contrary, Table 3 summarized critical exponents of the Ising model (D s = 1) for various lattice dimension d [8] .
First, it is notable that the first row of Table 2 is the same as the second row (d = 3) of Table 3, which is shadowed on the tables. Second, with increasing the lattice dimension d from 2 to 3, then 4 and higher, γ decreases from 1.75 to 1.25, then 1 finally. γ = 1 for 4 and higher dimension calculation matches exactly to the classic phenomenological model, which is popularly known as the "Curie-Weiss law". Third, with increasing the lattice dimension d from 2 to 3, then 4 and higher, the fractal dimension D increases from 1.75 to 1.96, then 2 in the Random Walk model or from 1.875 to 2.51, then 3 finally in the Percolation model.  From the analogy with the above mentioned, we reconsider the fact that the critical exponent γ of the permittivity determined in the temperature range above the transition point represents the degree of the phase transition diffuseness. The γ increase may be related with the fractal dimension D decrease.

Acoustic Emission (AE) Under Voltage
The electric field induced acoustic emission (AE) in ferroelectric ceramics can provide information on the domain dynamics and mechanical deformation, as well as cracking. This Section treats the normal AE, which is not related with the crack propagation. This is so-called "Kaiser" effect [4,7] .

Acoustic emission basics
Acoustic emission (AE) method is a non-destructive technique used to detect pulses of released elastic strain energy caused by deformation, crack growth and phase change in a solid. [19] The AE method has also been used to investigate fundamental domain motions in ferromagnetic and ferroelectric materials [20,21] , as well as mechanical and fatigue properties of materials. [22] In ferroelectric and piezoelectric materials, the AE method has been used to determine phase change [23][24][25] , to detect domain reorientation [26][27][28] and to monitor crack propagation. [29,30] However, our researches revealed a significant problem in most of the previous field-induced AE studies in the ferroelectrics. Acoustic waves were generated by a sample's mechanical vibration, which was caused by an electrical coupling between the power supply feedback and the piezoelectricity of sample. [31] The acoustic wave generated by the vibration was large enough to conceal the true (domain dynamics-related) AE signal and was constantly observed at high signal amplification of the AE measurement system. Improved AE measurement system: Figure 14: Improved Acoustic Emission (AE) measurement system developed in ICAT, Penn State University. [7] Figure 14 shows the measurement system developed in International Center for Actuators and Transducers (ICAT), The Pennsylvania State University. Refer to [4] and [7] for the details. The electric field induced AE and displacement can observed simultaneously. The AE signal is detected by an AE sensor with its resonant frequency of 450 kHz. The AE signal is amplified by 40 dB through a pre-amplifier and again up to 60 dB with a main amplifier. The pre-amplifier should be operated by batteries. Otherwise, the noise on the power cable conceals the real AE signal. The amplified signal is counted after passing through a high-pass filter (100 kHz) and a discriminator. Disc specimens of commercialized ferroelectric PZT ceramic (PZT-5A) as well as, electrostrictive ceramics 0.9PMN-0.1PT and soft piezoelectric PLZT (9/65/35), with 12.7 mm in diameter and 0.4 mm in thickness were employed for the measurements. Gold electrodes were formed on both surfaces by sputtering. The AE sensor was attached to the sample through brass-foil with 0.025 mm in thickness, and silicone grease was used as a couplant to make a close adhesion. An external resistor R (R = 20 MΩ) and a capacitor C (C = 0.01 μF) were installed to the power supply in order to eliminate the electrical coupling by increasing the time constant T (=R x C) of the voltage application system (i.e., to suppress the high frequency feedback components). It should be noted that the external resistor and capacitor did not cause any electrical interference to the AE measurement system and also the intrinsic AE generations in the ceramics as long as the induced displacement was maintained and the field was applied at much lower frequency than the time constant T. The threshold level and amplification of the AE signal were changed from 400 m V at 100 dB to 400 m V at 60 dB. The samples were driven under the bipolar field at various electric field levels. The driving frequency of 0.0015 Hz was employed to obtain the saturation ofthe field induced dis-placement and AE to avoid the effect of the supplied current limited by the external resistor and capacitor. In order to check complete elimination of the sample vibration and electrical noise, the fundamental "Kaiser Effect" in terms of the electric field should be monitored in ferroelectric ceramics prior to the detailed measurements [32] .Results are shown in Figure. 15 for a PZT disk sample. The field induced AE will not be present after poling as long as the applied field is maintained lower than that of the poling field (i.e., Kaiser Effect). The piezo-ceramic memorizes experienced electric field or previous history. Therefore, the previous reported ferroelectric domain contributions in the field induced AE generation should be re-examined through the new AE measurement technique. Figure 15: Kaiser effect observed in a piezoelectric ceramic. [32] The field induced AE will not be present after poling as long as the applied field is maintained lower than that of the poling field. Figure 16(a) and 16(b) show the AE event count rate and the induced displacement of the ferroelectric PZT ceramics as a function of applied field at 0.0015 Hz with field levels of Phenomena in Ferroelectrics ±25 kV/cm and ±35 kV/cm, respectively. The AE measurement conditions were set to be at 100dB with a threshold level of 400mV. A butterfly shape induced displacement was obtained. For both measurements, the AE was not observedat the coercive field at which the induced displacement took a minimum. The critical electric field, where the AE started to be generated, was slightly higher than the coercive field. This critical electric field corresponded to the point of inflection at which the second derivative of the displacement with respect to the field, d 2 (displacement)/dE 2 , changes the sign.  [7] In Figures 16(a) and 16(b), change of the first derivative of the displacement with respect to the field d(displacement)/dE is shown for the point of inflection, and the critical electric field corresponds to the point where |d(displacement)/dE| is the maximum. Since the displacement induction rate started to decrease in spite of the field simply increasing and the AE was present above this field, these results might indicate that a damping process for the sample deformation existed above this field and the AE was generated through that process. The AE event count rate increased with the applied field when a field of E = ±25 kV/cm was applied to the sample. For a field of E = ±35 kV/ cm, the maximum of AE event count rate was not observed at the maximum applied field. The maximum event rate occurred around 27.5 kV/cm. Considering that internal stress increases with the applied field in ferroelectrics, if this observed AE was caused only by the internal stress, the AE event count rate should have increased with the applied field. Thus, this decrease in the AE event rate might suggest that the internal stress could be the origin of the AE, but not only one source. Since the induced displacement in the ferroelectric ceramics consists of the domain reorientation related deformation and piezoelectric deformation without the domain reorientation, it was assumed that the field induced AE in the PZT ceramics was generated first through domain reorientation related deformation and, after domain reorientation was completed, the piezoelectric deformation unrelated to domain reorientation which was accompanied with induced stress was expected to be the origin of the AE. The decrease in the AE event might indicate the completion of the domain reorientation related deformation. Figure 17 shows the induced displacement and AE event count rate of the electrostrictive 0.9PMN-0.1PT ceramics. A quadratic induced displacement which showed much less hysteresis than that of ferroelectric PZT ceramics was observed. The AE count rate increased with the applied field. Because of the absence of the large domain in PMN-PT, the induced displacement might not be caused by the domain reorientation, but domain reorientation-unrelated electrostrictive deformation. It was assumed that the origin of the AE was the internal stress caused by the electrostrictive deformation unrelated to domain reorientation, resulting in the increased AE event rate with applied electric field.  Figure 18 shows the AE event count rate and the induced displacement of the soft piezoelectric PLZT(9/65/35) ceramics, which shows an intermediate displacement change between the ferroelectric PZT and electrostrictive PMN-PT ceramics. The AE generation started above the coercive field which was determined from the induced displacement. A peak of the AE event rate was observed around 10kV/cm, and the AE event rate increased with the applied field when a field higher than 10 kV/cm was applied. It was assumed that the peak of the AE event rate was due to the domain reorientation related deformation, and after the domain reorientation was completed, the AE was caused by the piezoelectric deformation unrelated to domain reorientation. These results suggests that both the domain reorientation related and unrelated deformations were the origins of AE in the PLZT (9/65/35) ceramics.

Fractal dimension of the electric field induced AE
A fractal analysis of the AE signal amplitude distribution, which is used in the AE method to determine the damage of materials [33] , is also applied for the ferroelectrics in our research. When the AE event rate f (x) is the minus D-th power of x, where x is the AE signal amplitude, The number of D is defined as the fractal dimension [32] .The integrated AE event F(x), is observed through changing the AE signal threshold level in the measurements. If D = 1, PZT: Figure 19 (a) Figure19(a) shows the AE event count per cycle as a function of the AE signal threshold level in the ferroelectric PZTceramics. The observed AE event count logarithmically decreases with the AE signal threshold level when a field of ±25kV/cm is applied Figure 19(a). Thus, a fractal dimension of D = 1 can be obtained for the logarithmic decrease in the ferroelectric PZT. Since repeatable and stable field induced AE is observed, it is supposed that the lowered fractal dimension is due to the existence of ferroelectric domain and the effect of the do-main reorientation related deformation, instead of the dam-age. When a field of ±35 kV/cm is applied, the higher AE signal threshold level parts does not follow to logarithmic change. The lower signal threshold level region still could be fitted to a logarithmic curve and the extended line also ends around the AE signal threshold level obtained for ±25 kV/cm. The fractal dimension of higher signal threshold level part can be found to be D = 2.8. Since the domain reorientation is supposed to be completed at higher electric field, the origin of the AE with the fractal dimension of 2.8 is assumed to be piezoelectric deformation without domain reorientation. It is shown that the fractal dimension of the field induced AE for the domain reorientation related deformation (D = 1) is lower than that for the piezoelectric deformation without domain reorientation (D = 2.8). Figure 19(b) shows the AE event count per cycle for the electrostrictive PMN-PT as a function of the AE signal threshold level. The same slopes for the AE event count are observed at various applied fields, and the critical fractal dimension of D = 2 can be obtained. It seems that the origin of AE in the electrostrictive PMN-PT ceramics does not change with the applied field level, and was assumed to be the internal stress caused by the electrostrictive deformation unrelated to macro-domain reorientation, but by micro/nano domain dynamics.  [7] www.ommegaonline.org 21 Phenomena in Ferroelectrics J Nanotech Mater Sci | volume 1: issue 1 Uchino, K., et al.

PLZT:
In the PLZT (9/65/35) ceramics, the AE event count fluctuated with the applied field Figure 19(c). This fluctuation in the AE is probably due to the field induced phase change and/or domain creation. It seems that the domain reorientation related deformation, the piezoelectric deformation unrelated to domain reorientation and the field induced domain creation as well as phase change related behavior cause the AE generations in the PLZT (9/65/35) ceramics. The fractal dimension of the PLZT (9/65/35) does not take a constant intermediate number between D = 1 and 2, but fluctuated. This is probably due to the change of AE generation source caused by the applied field. If we use 27.5 kV/cm AE values (close to linear), D = 1.7 can be obtained.

In partial conclusion
(1) The AE event rate in ferroelectric PZT and soft piezoelectric PLZT (9/65/35) ceramics show two origins: domain reorientation related deformation and piezoelectric deformation without domain reorient-action. On the contrary, the internal stress caused by the induced electrostrictive deformation (micro/nano domain dynamics) is supposed to be the dominant AE generation source in the electrostrictive0.9PMN-0.1PT ceramics.  An intuitive domain dynamics model is illustrated in Figure. 20 for explaining the fractal dimension. Figure 20(a) depicts macro-domain reorientation in PZT under low electric field. Clear plate-like 90º-domain walls move primarily in perpendicular to the plate (1-D motion). If the AE comes from this motion, the fractal dimensions D = 1 can be explained. However, when the electric field is larger than coercive field in PZT, mono-domain state will be realized. The deformation beyond this field level is caused by intrinsic piezoelectricity (elastic deformation); that is, 3-D deformation (1-D extension associated with 2-D shrinkage). This suggests the fractal dimension D = 2.8, close to the lattice dimension 3. On the contrary, in PMN-PT, spindle-like micro-domains breathe (i.e., fatter and slimmer like a rubber tube) by keeping the length, as shown in Figure 20(b) We can imagine D = 2 from this micro-domain dynamics.

Acoustic Emission (AE) Under Cracking
Though the damage or crack in ferroelectric ceramics is not our primary interest in this paper, brief description on "Felicity effect" may help with understanding the AE data practically measured.

Multilayer piezoelectric actuator
Multilayer ceramic actuators (MCA) are used in many applications, such as diesel injection valves for automobiles and probe positioners in scanning tunneling microscopy (STM), because of high force, electromechanical coupling and micro-displacement controllability. Thus, our primary research target in their reliability was to develop the monitoring system for clarifying the crack mechanisms. The acoustic emission (AE) method was adopted for this purpose [5,6] .  Figure 21 shows the MCA sample specifications. MCA of soft PZT is usually operated by a unipolar voltage from 0 to 60 V, which corresponds an electric field of 1.5 kV/mm, and exhibits an induced displacement of 1.6 μm. Since high local concentrated electric field around the edges of the internal electrodes cannot be determined, the applied voltage rather than field is used in this Section. Note that the internal electrode and surface buffer layers were designed intentionally for introducing the crack.  Figure 15, the first deformation related AE at poling was not observed. AE signals in the MCA might be easily attenuated by the internal electrodes, because of the acoustic impedance mismatch. On the contrary, when a poling voltage with a peak of 200 V was applied, the AE event count rate started to increase around 140 V Figure 22(b), and continuously observed even when the voltage was decreased. This is against the normal "Kaiser Effect", and considered to be higher internal stress/ strain and consequential cracking. We can conclude that the critical poling voltage for the actuator, which does not cause any major cracking, can be determined using the AE measurement.

Felicity effect
When elastic materials are damaged, the re-applied stress, where new emissions start, becomes lower than the previous maximum stress (The Felicity effect). [34] This decrease in the stress at the onset of AE count can be caused by a friction between free and damaged surface. The Felicity ratio is defined as follows: Felicity ratio = σ stress at onset of AE /σ previous maximum stress (36) 1st-cycle, (b) 2nd-cycle and (c) 10th-cycle. [5] Figure 23 shows the induced displacement and AE event count rate for a unipolar cyclic test with a peak voltage of 180 V. In the first cycle, a large number of AE event counts started from 160 V and could be caused by both stress and cracking Figure 23(a). The AE observed during the applied voltage decrease is proba-bly due to the friction at the damaged and un-damaged surfaces (i.e., delamination). In the second cycle, the total AE event count per cycle decreased. This is assumed to be the stress relaxation in the MCA due to cracking. A larger AE event rate than that of the first cycle was observed from 130 V to 160 V Figure 23(b). At the 10 th cycle, even the overall AE event decreased, the AE onset voltage stayed around 130 V Figure 23(c).
During the cyclic test, the total AE event count decreased with number of applied voltage cycle, but the AE onset voltage was found to be around 130 V. It is supposed that the damage in the MCA sample lowered the durability against applied voltage. This damage effect on the AE generation can indicate the Felicity effect in the MCA and the Felicity ratio is calculated to be 130 V/180 V = 0.72. The Felicity ratio in terms of the applied voltage can be used to evaluate the existing damages in the MCA and to determine the maximum safety driving voltage for the MCA which will not cause any major cracking.
After measuring the AE event count and observing the Felicity effect in the sample, the MCA was carefully sliced and polished for microscopic observation. Figure 24 shows a micrograph of the crack actually generated in the sample poled at 200 V. Since the crack was observed in the vicinity of the internal electrode edge at an active layer next to the inactive top layer, the stress concentration caused by the inactive layer should be the origin of the cracking. The failures which can raise the destruction probability of the MCA can be detected during the poling process using the AE monitoring.  [5] Fractal dimension for cracking During the above Felicity effect measurement, we have not paid attention to carefully calculate the fractal dimension. Thus, we will cite previous researches here for the references.
Hirata reported the compressive destruction test in Granite cylinder samples [35] . Rocks usually take three stages until total collapse: (1) First creep -fast speed deformation, (2) Second creep -constant and gradual speed deformation, then (3) Third creep -accelerated speed deformation, reaching to the fatal destruction. Hirata used AE measurement associated with these creep processes to generate the AE spatial distribution www.ommegaonline.org map, then calculated "accumulated 2-point correlation factor" to determine the fractal dimension. The accumulated 2-point correlation follows a power law for all three stages, which indicated the fractal phenomena. The exponent provides the fractal dimension of the AE spatial distribution: first creep -D=2.75, second creep -D=2.66, and third creep -D = 2.25. The fractal dimension decreases with the deformation progress, and the decrease is significant immediately prior to the fatal destruction. The original space structural fluctuation of rocks has a fractal dimension around 2.8, which indicates rather uniform random distribution. As the deformation progresses, the correlation between micro cracks gradually increases and the cracks cluster prior to the fatal collapse, as monitored in fractal dimension decrease [35] .
Dielectric breakdown phenomena may also be close to the piezoelectric actuator cracking. The dielectric breakdown is a successive growth phenomenon of the low-resistive part, which is initially generated by high voltage. Niemeyer et al. [36,37] simulated this growth process for a circular discharge phenomenon under a 2-D square lattice model: (1) Potential on the circle circumference (cathode) ϕ = 1 V, potential at the circle center (anode) and at the connected lattice points ϕ = 0. Potential at other (not connected) lattice point's ϕ can be calculated with the discrete Laplace equation. (2) The growth process (extending the connected lattice cluster to an adjacent lattice point) is decided by the random probability: P(r) = |Δ ϕ(r)| η /∑|Δϕ(r)| η . (37) For a case η = 1, the fractal dimensions D = 1.75, and D = 2.48 are reported, respectively for the Euclidean dimension d = 2 and 3.

Fractal-Related Light Scattering Studies
We review two representative papers on light scattering on PLZT and PMN in this chapter.

Light Scattering Intensity and Fractal Dimension
Schuh et al. analyzed the change on time in longitudinal elastic scattering of monochromatic light (λ = 0.63 μm) under unipolar rectangular electric field pulses (duration ranging from 1 ms to 10 s) [38] . It was shown that under application of pulse series the time dependence of light intensity I(t) scattered at a fixed angle ϕ is reproducible in details. This feature permitted to determine an instantaneous scattering indicatrixI(ϕ)at any time from a series of consecutive measurements of I(t) for light scattered inside various fixed angles {ϕ} selected by the set of ring diaphragms in angular range 0.5-7 degrees. The variation of the integrated scattered light intensity in the relaxor phase under application of a rectangular field pulse can be separated into two stages: at the switching stage (after field switch on) the intensity increases to I max (E); during back-switching stage after field removing the scattered light intensity spontaneously recovers the initial background level I min .
normal-transition ferroelectrics reveals an intermediate γvalue of 1.58. The critical exponent γ of the permittivity, determined in the temperature range from several to several hundreds °C above the transition point, may at least represent the degree of the phase transition diffuseness.

Fractal dimension in ferroelectrics
The fractal dimensions determined by acoustic emission (AE) in various normal and relaxor ferroelectrics are summarized in Table 4.The fractal dimension of the normal ferroelectric PZT ceramics consisted of the domain reorientation related D = 1 and the domain reorientation unrelated D = 2.8, when the electric field of ±35 kV /cm is applied. The fractal dimension of D = 2 is obtained in the electrostrictive PMN-PT (i.e., relaxor). The fluctuation of the fractal dimension is observed in PLZT (9/65/35) ceramics, which show an intermediate state between PZT and PMN-PT in terms of the induced displacement. The fractal dimensions in PLZT and PMN determined by light scattering experiments are also inserted for the reader's reference. Note that Koreeda's result on PMN was obtained from un-poled sample with micro-domains; our data from AE were from the poled PMN-PT or under a large electric field applied condition.  [38] ; # Reference [39] Though we do not have suitable theoretical tools/models at present to correlate the permittivity exponent γ and fractal dimension D, it is noteworthy that the systematic γ and D value change with phase transition diffuseness indicates a strong relationship. With promoting the domain size growth (from micro to macro) with temperature decrease and/or with external electric field increase, γ decreases from 2 to 1 and D increases from 2 to 3 (mono-domain state). However, we also point out that the critical exponent is the results for low voltage (domain reorientation unrelated), while the fractal dimension is the AE results measured under high electric field (domain reorientation related). Thus, the meaning of the fractal dimension seems to be different, as exemplified by D = 1 for large ferroelectric domain reorientation.