Piezoelectric Nonlinearity and Hysteresis Arising from Dynamics of Electrically Conducting Domain Walls

Macroscopic nonlinearity and hysteresis observed in the piezoelectric and dielectric responses of ferroelectric materials to external stimuli are commonly attributed to localized displacements of domain walls (DWs). The link between the macroscopic response and microscopic DW dynamics is provided by the well-known Rayleigh relations, extensively used to quantify the electrical and electromechanical behavior of ferroelectric ceramics and thin films under subswitching conditions. In this chapter, I will present an intriguing case where DWs exhibit enhanced electrical conductivity with respect to the bulk conductivity. By combin-ing experimental data and modeling, it will be shown that the local conductivity, related to accumulation of charged points defect at DWs, does not only affect DW dynamics through DW-defect pinning interactions, as we may expect, but goes beyond it by affecting the macroscopic nonlinearity and hysteresis in a more complex manner. The major characteristics and implications of the underlying nonlinear Maxwell-Wagner piezoelectric relaxation, triggered by the presence and dynamics of conducting DWs, will be presented, reviewed and discussed in the framework of systematic multiscale analyses on BiFeO 3 ceramics. The result may have implications in the development of promising BiFeO 3 -based compositions for high-temperature piezoelectric applications.


Introduction
Historically speaking, probably the most important examples of dielectric and piezoelectric nonlinearity and hysteresis in ferroelectrics can be found already in the initial lines of the introductory section of the paper by Joseph Valasek from 1921 [1], which we use today to celebrate the 100th anniversary of the discovery of ferroelectricity [2]. Despite struggling with moisture-sensitive Rochelle salt crystals, Valasek finally managed to publish the first ever case of a ferroelectric hysteresis loop. Today, such hysteresis of the polarization response to the external electric field, originating from ferroelectric domain switching, is essential in, e.g., ferroelectric random-access memories (FeRAMs) [3]. Hysteretic and nonlinear responses to external fields, however, are not only observed at switching conditions. Deviations from simple linear relationship between polarization (P) and electricfield (E) (dielectric response), P and mechanical stress (Π) (direct piezoelectricity) and mechanical strain (x) and E (converse piezoelectricity), are very common at subswitching driving conditions at which many important devices operate, such as capacitors, sensors and actuators. Such responses are measured in both bulk and thin-film ferroelectrics [4][5][6].
The origins of subswitching macroscopic nonlinearity and hysteresis alongside the frequency dependence of the property coefficients, which are directly relevant for device performance, are difficult to assess due to complex interrelated mechanisms often operating simultaneously over a wide length-scale range of the material. By risking of being too general, one could say that the major players of these mechanisms are domain walls (DWs) and similar interfaces [7], charged point defects [8][9][10][11], grain boundaries [12], spontaneous ferroelastic strains [13,14], oxygen octahedra tilts [15] and secondary phases [16], if only the widely discussed examples are listed. Most of these features are always present in polycrystalline ferroelectric materials and are capable of cross-interacting via electric and/or elastic fields (a classical example that every ferroelectrician should be aware of is the DW-defect electro-elastic interaction [17,18]). Complicating the picture is the fact that many of these features are strongly dependent on the synthesis conditions. Despite playing a key role, some of them, such as point defects, could not be directly visualized until recently [19][20][21][22], meaning that the interpretations of their role in the macroscopic response are supported by either indirect analysis or other reasoning (see next section). Without extensive knowledge in the area and with just a little intuition, it is pretty easy to understand why dielectric and piezoelectric nonlinearity as well as hysteresis in ferroelectric and similar materials are so difficult to understand and predict. Based on my personal opinion, the results shown in this chapter are the perfect example to illustrate the difficulties in elucidating the mechanisms responsible for the piezoelectric nonlinearity and hysteresis in a complex material, such as BiFeO 3 (BFO). The results that I am going to present have been acquired over eight long years of multidisciplinary research focusing on theoretical and experimental studies and involving a large number of researchers from different fields.
One of the most important and widely-recognized origins of nonlinearity and hysteresis in the subswitching electrical and electromechanical response of ferroelectrics is the displacement of domain walls (DWs) [5,[23][24][25][26]. These are interfaces separating domain regions in which the spontaneous polarization is directed in one of the symmetry-allowed directions. Displacements of DWs are mechanistically treated via the interaction with charged points defects, which act as pinning sites and, depending on their type, mobility and location, affect the movement of DWs under applied fields (this issue will be thoroughly discussed in the next section). Reversible and irreversible motion of DWs in a ferroelectric with randomly distributed pinning sites (defects) is described by the Rayleigh law (RL), which is vastly used to quantify and cross-compare the subswitching responses of ferroelectric ceramics and thin films [6,27]. For the converse piezoelectric effect (same is valid for the direct effect or dielectric response), the RL can be expressed by a set of two equations: For simplicity, the subscripts of the piezoelectric coefficient are omitted. Eq. (1) describes the linear dependence of the real component of the piezoelectric coefficient (d 0 ) on the electric-field amplitude (E 0 ) with d init and α representing the reversible and irreversible Rayleigh coefficients, respectively. Note that d init corresponds to the coefficient extrapolated at zero field or d init (2) is a multivalued function and describes the strain (x) versus electric-field (E) relationship where the positive and negative signs, separating the first and the second term of the equation, are related to the strain branches defined by descending and ascending E, respectively. The second quadratic term of Eq. (2) thus defines the hysteresis. Note that nonlinearity refers here to the nonlinear relationship between strain (x) and electric field (E) (second quadratic term of Eq. (2)), in which case the proportionality coefficient, d', is linearly dependent on the field amplitude (Eq. (1)).
By carefully inspecting the Rayleigh equations, it is not difficult to understand that the essence of this model is the intimate relationship between the nonlinearity and hysteresis. For example, in Eq. (1) the irreversible coefficient α defines the nonlinearity by the field dependence of the piezoelectric coefficient, however, α also appears in the second term of Eq. (2) describing the hysteresis (by the different up-field and down-field strain branches). Alternatively, one can easily realize the nonlinearity-hysteresis relation of RL in the second term of Eq. (2), which sets the nonlinear (quadratic) relation of strain to field while, at the same time, representing the hysteresis. In other words, hysteresis arises from nonlinearity and vice versa. A mathematical test is straightforward: in the absence of the irreversible contribution (α = 0), there are two important consequences: , meaning that the response is linear (i.e., the coefficient is constant), and (ii) the hysteresis vanishes (Eq. (2) reduces to the first anhysteretic term). Therefore, every irreversible DW displacement contributes to nonlinearity and hysteresis in a unique way set by the RL [27].
When mechanistically reasonable, the irreversible coefficient α, experimentally extracted from the measurements of, e.g., longitudinal converse piezoelectric coefficient d 33 as a function of field amplitude E 0 (se Eq. (1)), is commonly used to quantify (irreversible) DW contribution to piezoelectric properties (an analogous approach is used to quantify DW contribution to dielectric permittivity). A rigorous test to find out whether the measured data can be approximated by RL is to at least validate the nonlinearity-hysteresis relationship implied by RL. This can be done by verifying whether the experimental hysteresis can be fitted by Eq. (2) using the coefficient α that is determined from the d 33 versus E 0 slope [6,27]. Importantly, other parameters of the measured response can be checked against the RL predictions. For example, using the Fourier-series analysis of Rayleigh equations (derivations can be found in Ref. [27]) and the hysteresis-area analysis in complex mathematical formalism (see, e.g., references [24,28]), it can be shown that the ratio between the increment in the real piezoelectric coefficient ( 3π ) and does not depend on the external field amplitude. This can be easily verified if the experimental piezoelectric coefficients and piezoelectric phase angles (or dielectric permittivity and losses) are known for different driving fields [14]. Other response parameters that can be evaluated against the RL predictions are embedded in the third harmonic response; this is discussed in detail in Ref. [29].
It has been shown that the subswitching response of a number of ferroelectric materials well obey the RL. Examples include donor-doped soft Pb(Zr,Ti)O 3 (PZT) ceramics [5,27,28,30], coarse-grained undoped BaTiO 3 [31], high-Curie-temperature piezoceramics based on BiScO 3 -PbTiO 3 [24], textured (K 0.5 Na 0.5 )NbO 3 (KNN) ceramics [32], Aurivilius-type Nb-doped Bi 4 Ti 3 O 12 [31], some relaxor-ferroelectric Pb(Mg 1/3 Nb 2/3 )O 3 -PbTiO 3 compositions [33,34] and Pb-based thin films [35,36]. It has to be pointed out, however, that even in these cases, strictly speaking, a good match between the Rayleigh model and the measured subswitching response is typically observed in a limited driving field range, sometimes referred to as the "Rayleigh range" [5,27]. This is expected, considering that RL is an idealization and assumes DW motion in a perfectly random pinning potential. It is clear that real materials' responses may come close to this situation; however, in most cases deviations from RL predictions are naturally observed and should not be surprising. It is legitimate to think that these deviations have microscopic origins different from DW motion. However, this is not necessarily the case. A nice supporting example can be found in soft PZT, which often exhibit sublinear, quasi-saturating field dependency of the piezoelectric coefficient at relatively weak fields, clearly deviating from the perfect linear behavior predicted by Eq. (1) [4]. In this case, Preisach approach becomes useful as it shows that such response can be understood by considering a non-uniform pinning potential where the concentration of weakly pinned DWs is higher than those that are strongly pinned (note that the two should be equal in a perfect Rayleigh random pinning potential, resulting in a flat Preisach distribution function) [37]. Therefore, while the experimental data deviates from the RL, in some cases this deviation can still be explained by DW motion without necessarily invoking other mechanisms unrelated to DWs.
Apart from those cases where the response is close to RL predictions, the subswitching response can be clearly non-Rayleigh. Acceptor-doped hard PZT [30,38], Sm-doped PbTiO 3 [27] or monoclinic Pb(Mg 1/3 Nb 2/3 )O 3 -PbTiO 3 (PMN-PT) ceramic compositions with relaxor features [34] are just few of such examples. In this chapter, I am going to present another case that has attracted particular attention in recent times. The case is on BFO, a perovskite that is unique for at least two reasons: (i) because it contains electrically conducting DWs that are formed spontaneously and are attractive for applications in nanoelectronics [39] and (ii) due its high Curie temperature ($835°C) [40,41], which makes it the key perovskite in the development of novel compositions for high-temperature piezoelectric applications [42]. As it will be shown in the next section, enhanced conductivity at DWs completely alters the DW dynamics under applied subswitching fields. I will give evidence of a very large nonlinear contribution to the total piezoelectricity in BFO arising from irreversible DW dynamics, which is strongly dependent on the frequency of the external driving field. Along with a peculiar hysteretic behavior, characterized by a negative piezoelectric phase angle, and an unusual clock-wise rotational sense, I will strive to demonstrate how complex the dynamic response of electrically conducting DWs can be and how careful should we be in analyzing and treating materials that show unique behaviors. Somewhat surprisingly, I will show that despite the charged point defects, which are accumulated at DWs and are responsible for the local conductivity, act as pinning sites (as explained in Section 2), they can even increase the DW mobility and largely contribute to the enhanced nonlinear converse piezoelectric response. This may happen because the local conductive paths through the grains, set by the conductive DWs, lead to a redistribution of internal fields, effectively resulting in grains or grain families inside which the electric field is largely enhanced with respect to the field applied externally. The mechanism, called nonlinear Maxwell-Wagner piezoelectric effect (explained in Section 3), can be supported by simple analytical modeling (Section 4) and should be more seriously considered when engineering BFO-based materials for high-temperature piezoelectric applications.

Dynamic interactions between charged point defects and domain walls: a microscopic view
In this section, we will take a look at idealized microscopic pictures describing the interactions between DWs and charged point defects, which have profound implications in the DW dynamics and play a key role in the macroscopic dielectric and piezoelectric nonlinearity and hysteresis. Before discussing these interactions, however, it is timely to emphasize few important points. First, the reader should be aware that DW dynamics is not the only mechanism giving rise to electrical and electromechanical hysteresis. An example that should be familiar to every scientist and engineer working in the field of dielectrics or capacitors, is the DC electrical conductivity, which can directly contribute to the imaginary dielectric coefficient and thus dielectric losses and P-E hysteresis. In ionic compounds, electron and holes tend to be trapped at local sites in the lattice; they locally polarize and deform the lattice, creating the so-called polarons. Being self-trapped, these charges can move via an activated hopping process, similar to ionic conductivity [43]. Interestingly, in both physical and mathematical sense, hopping charge transitions are indistinguishable from pure dipole reorientations, meaning that hopping conductivity can contribute not only to losses (imaginary dielectric permittivity) but also to polarization (real dielectric permittivity). Mathematical descriptions of these problems are given in the book by Jonscher, which is highly recommended for learning conductivity phenomena in dielectrics [44].
Second, the effect of the electrical conductivity can go beyond a simple contribution to the complex dielectric permittivity and can strongly affect piezoelectricity, too. This case, which is less trivial and somewhat difficult to digest (even for experts in the field), is the Maxwell-Wagner (M-W) piezoelectric relaxation [45][46][47][48][49]. As an analogue of the well-known dielectric M-W relaxation [50], it originates, generally speaking, from electrically conductive paths present locally inside the material, which lead to internal electric-field redistribution as a function of the driving field frequency and, consequently, to electromechanical relaxation phenomena. The effect was predicted for the direct [4,48] and converse piezoelectric response [51], and even for purely elastic response [52]. It was experimentally observed in heterogeneous ceramic materials [48], two-phase polymer systems [53] and piezoelectricpolymer composites [47]. The piezoelectric M-W effect is described in detail in an older book on hysteresis [4], to which the reader should refer to as a support to the data presented in this chapter where the focus is on BFO ceramics.
To explain the dielectric M-W effect, either grain boundaries or sampleelectrode interfaces are commonly considered, as these are the regions where depletion layers are likely to be formed, giving rise to inhomogeneous field distribution inside the material and, consequently, to M-W relaxation [50]. In this chapter, I will present a different case where local conductive paths originate from charges concentrated at DWs. It is intuitively understood that this case must be unique because, unlike grain boundaries or sample-electrode interfaces, the conductive paths provided by the DWs are dynamic: they can move and locally displace under external fields. One may simplistically view this scenario as a local conductivity that is not spatially confined.
The third point to consider is rather general. While the discussion in this chapter is mostly about DWs, the reader should understand that also the dynamics of other types of interfaces may lead to the same macroscopic effects associated with subswitching nonlinearity and hysteresis. As an example, those other interfaces may be represented by boundaries separating regions consisting of phases with different crystallographic symmetries, which are present in morphotropic compositions [7]. The role of interface boundary motion during subcoercive loading of morphotropic BiScO 3 -PbTiO 3 ceramics was recently demonstrated by in situ X-ray diffraction (XRD) analysis [54].
We come back now to our case on DWs. In analogy to the well-known pinning of DWs in ferromagnetic materials [55], the same concept is used for ferroelectrics [56]. An important difference is that, on average, ferroelectric DWs are an order of magnitude thinner than ferromagnetic DWs [39,[57][58][59]. While it is clear that the subject is complex and cannot be reduced to few (over)simplified pictures, like those shown in Figure 1, it is surprising to realize that the behavior of many ferroelectric materials under applied fields can be qualitatively interpreted by considering relatively simple interactions between DWs and charged point defects. This is particularly true for high-field P-E (switching) hysteresis (see, e.g., reference [18]) and, to some degree, for subswitching responses (see, e.g., reference [4]). Precaution is, however, necessary because there are various parameters to consider, such as type, concentration, location and mobility of defects along with their binding state, as well as the presence of other pinning centers, including grain boundaries [12] or oxygen octahedra tilts [15]. Nevertheless, a simple exercise linking possible DW-defect interactions with the macroscopic response, which will be explained next, is very useful to assess the behavior of a given ferroelectric material.
Early studies identified three microscopic scenarios to describe the strong pinning and DW stabilization effects observed macroscopically by pinched and/or biased high-field P-E hysteresis loops of acceptor-doped PZT and BaTiO 3 [8][9][10][11]17]. The three scenarios are illustrated in Figure 1a and assume different locations of the pinning centers: (i) inside the domains (called "volume" or, sometimes, "bulk" effect), (ii) in DW regions ("domain wall" effect) and (iii) in or close grain boundary regions ("grain boundary" effect). The common feature of all these scenarios is that defects are arranged in a somewhat ordered manner and, by that, they stabilize the domain configuration via electric and elastic coupling [11,18], inhibiting domain wall motion. Ordered defects are a characteristic of the so-called "hard" ferroelectrics, exemplified by acceptor-doped PZT. In contrast, "soft" ferroelectrics, represented by donor-doped PZT, are usually described by assuming disordered defects as shown in Figure 1b. Hard-soft transitions, induced in hard-type compositions by disordering the otherwise ordered defects, can be achieved by either electric-field cycling [9,60] or quenching [30,61,62].
The "volume scenario" shown in Figure 1a (left schematic), which assumes binding of defects into complexes (see black arrows), is probably the only DWpinning mechanism that is highly accepted in our research community and is widely used to explain hardening and aging in acceptor-doped PZT, BaTiO 3 and similar perovskites (an excellent review on this topic can be found in the paper by Genenko et al. [63]). DW pinning effects mediated by orientation of defect complexes have been put in a theoretical framework more than 30 years ago [17]. More recent

Simplified microscopic pictures of interactions between charged point defects and domain walls in (a) hard and (b) soft ferroelectrics. P s and P d are spontaneous and defect polarization, respectively. Thick black lines denote
DWs and circles noted with + andsigns represent individual charged point defects (see also legend on the righthand side of the schematics). To be further noted is that P d represents the well-known acceptor-oxygen-vacancy defect complexes (see text for details). The pictures shown in panel (a) were adapted according to the pioneering studies by Jonker [10], Carl and Hardtl [9], Lambeck and Jonker [8] and Robels and Arlt [11]. The scenario in panel (b) assumes disordered defects in soft materials as discussed in, e.g., reference [4]. studies confirmed the formation of defect complexes and their orientation kinetics both theoretically (by first principles, [64]) and experimentally. In particular, electron paramagnetic resonance (EPR) spectroscopy was found to be a powerful tool due to its high sensitivity, not only to the binding state of certain defects (for details see [65]), but also to the orientation of the defect complexes in the material [65,66], allowing thus to directly monitor the alignment kinetics of the complexes under the applied field [67]. In most cases, the identified complexes are of the acceptoroxygen-vacancy type. The literature contains many examples of lead-based and lead-free ferroelectric compositions where such complexes have been identified by EPR; the list includes (but is probably not limited to) PbTiO 3 , PbZrO 3 , PZT, BaTiO 3 , KNN and (Bi,Na)TiO 3 (BNT). The dopants in those cases are typical B-site acceptors, such as Fe, Cu and Mn [65,66,[68][69][70][71].
The "DW scenario", which assumes defects accumulated in DW regions (Figure 1a, middle schematic), was first predicted by Postnikov more than 50 years ago [72]. Note that the original model does not assume defects bound into complexes. Interestingly, despite being seriously discussed in studies on fatigue mechanisms [73], the DW scenario received strong criticism (to be even labeled as "a wellknown speculation") from those that were strictly in favor of the volume scenario [74]. Nevertheless, very recent studies on BFO ceramics using atomic-resolution electron microscopy, which will be discussed later, support the idea of this scenario and indeed suggest that is active in this material [75].
The "grain-boundary scenario" (Figure 1a, right schematic), often referred to as the space-charge mechanism, was extensively discussed in a combined theoretical and experimental frame against the volume scenario originally proposed by Arlt et al. [17]. Key hardening characteristics of acceptor (Fe) doped PZT, such as the dopant-concentration dependent aging time, were shown to be well predicted by the space charge model [76,77]. A perhaps interesting observation is that phenomenological calculations revealed that the electrostatic pinning effects on DWs from the space charges can be two orders of magnitude stronger than the pinning effects provided by the aligned defect complexes (for the same charge-carrier concentration) [76]. Considering that grain boundaries are strained areas and sources of polarization discontinuities where point defects may be attracted due to electrostatic and elastic driving forces, the grain boundary scenario should probably be always considered when interpreting the macroscopic responses of ceramics to external driving fields, particularly if the ceramics are fine-grained.
It is clear that more than one of the three domain-stabilizing effects shown in Figure 1a may be active in a given material. Actually, the reason for considering only a single scenario is in the human nature seeking for the simplest explanation. A number of arguments supports the idea of multiple pinning mechanisms. I give two examples. In hard PZT, it was predicted by first principles that the acceptoroxygen-vacancy defect complexes should energetically prefer to be situated closer to DWs (as shown schematically with black arrows in Figure 1a, left schematic) [78]. This means that the pinning mechanism may consist of a combination of volume and DW effects. In another study based on conductive atomic-force microscopy (c-AFM) analyses, it has been shown that in BFO ceramics both DWs and grain boundaries exhibit enhanced electrical conductivity with respect to the bulk conductivity measured in the grain interiors [79]. This suggests the presence of mobile charges at both locations, which may act as pinning centers. In this same material, defect complexes based on oxygen vacancies have been also identified by EPR [80]. The overall data, therefore, point to an extremely complex situation where all three pinning scenarios shown in Figure 1a may be active. The challenge is to find whether and which mechanism is dominant. At present, conductive DWs appear to be the key features affecting the piezoelectric response of BFO, which will be discussed in detail in the next section.
In contrast to ordered defects in hard ferroelectrics, disordered defects in soft counterparts, as shown schematically in Figure 1b, are more difficult to be directly identified. Nevertheless, in the case of PZT, for example, the presence of disorder in donor-doped samples (in a pragmatic sense) is supported by the fact that the measured macroscopic subswitching nonlinearity and hysteresis are well consistent with the major predictions from the Rayleigh model [27]. Another argument is that EPR data on Gd-doped PbTiO 3 indicate no binding of the donor (Gd) dopant with the expected compensating Pb vacancies [81], supporting the results of firstprinciples studies on the same perovskite [78]. In addition, none of the two defects (donor dopants and Pb vacancies) are mobile below the Curie temperature of the material, where the defects can in principle be ordered because they are provided by electric and elastic driving forces originating from the spontaneous polarization and strain, respectively (see, e.g., the graphic explanations by Ren [18]). If the two defects are not mobile nor they show tendency of binding, then a "frozen-in" defect disorder state in donor-doped PT or PZT can be envisioned. It is not unusual, however, that the dopants segregate at the grain boundaries (a nice recent example is shown for Ti-doped BFO in Ref. [82]). The issue of defect segregation at or close to grain boundary regions should be more seriously considered, as recently pointed out by Slouka et al. [83].
Other discussions related to the nature of defects in soft PZT point to the likely possibility that the donor dopant reduces the concentration of oxygen vacancies in PZT (due to charge compensation), leading to a progressive transition from a state characterized by ordered defects (undoped PZT) to a state governed by disordered defects, as the donor dopant is added to PZT [29]. In this sense, experimental data even indicate that the dominant pinning centers affecting DW displacements and, consequently, nonlinearity, in both hard and soft PZT, could be oxygen vacancies. The lesson that can be learned from all these data is that the defect arrangements shown in Figure 1 should not be treated individually; as a matter of fact, all of them may be present, to some degree, in a given ferroelectric. Finally, it has to be emphasized that the situation in soft materials is more complex, and the true origins of softening are still not clear (if the reader is interested in these issues, it is recommended to consult the work of Dragan Damjanovic [29,78]).

Piezoelectric nonlinearity and hysteretic response of BFO ceramics: low-frequency nonlinearity and negative piezoelectric phase angle
One of the distinct characteristics of the piezoelectric response of BFO is the strong dependence of the nonlinearity on the frequency of the driving field. This is illustrated in Figure 2, which compares the converse piezoelectric response of BFO with that of hard and soft PZT. The PZT compositions correspond to the rhombohedral symmetry and were selected exactly to enable a comparison with the isostructural rhombohedral BFO, although symmetry is probably the only common aspect of these two perovskites. In the first place, it should be mentioned that the increasing tendency of d 33 with the electric-field amplitude E 0 observed in Figure 2a-c for the two PZTs (at all measured frequencies) and BFO (mostly at 1 and 0.2 Hz) is accompanied by an increase in the piezoelectric tanδ (not shown here). This suggests irreversibility in the response as implied by the intimate nonlinearity-hysteresis relationship predicted by RL (see Eqs. (1) and (2)). The nonlinear behavior in these samples can be thus related to irreversible domain wall contribution. This is supported by in situ XRD analysis on BFO [84] and similar PZT samples to those shown here [26].
As expected and described in the previous section, due to stronger DW pinning effects, hard PZT shows a lower absolute d 33 value, which is less dependent on the field magnitude (Figure 2a) than that of soft PZT (Figure 2b). Also, the d 33 field dependency in soft PZT is sublinear with respect to the rather linear dependency in hard PZT. This behavior is retained at all frequencies used in the measurements. The sublinear trend in soft PZT has been discussed to some degree in the introduction of this chapter in relation to non-uniform pinning potential.
The next information can be obtained by quantifying the data. The slope of the d 33 -vs-E 0 curves is for hard PZT $0.2 Á 10 À16 m 2 V À2 and is practically independent on the frequency. By contrast, in soft PZT, the slope (in this case calculated for 1 kV cm À1 , which corresponds to 10% of coercive field E C , i.e., E/E C = 0.1) is an order of magnitude higher, i.e., 8.1 Á 10 À16 m 2 V À2 at 90 Hz, and further increases to 9.9 Á 10 À16 m 2 V À2 at 0.2 Hz. The frequency dependence of the irreversible coefficient of soft PZT has been reported earlier and correlated with interface pinning in a disordered medium [85]. As easily assessed from Figure 2c, the nonlinearity and the associated irreversible contribution to piezoelectricity in BFO is completely different than in PZT. Here, the nonlinear behavior is strongly dependent on frequency: the d 33 -vs-E 0 slope is nearly zero at 100 Hz, precisely 0.003Á10 À16 m 2 V À2 , and increases to 0.27Á10 À16 m 2 V À2 when the frequency is reduced to 0.2 Hz (as in the previous case on PZT, the slope was calculated for the relative field E/E C = 0.1). Such a dramatic two-orders-of-magnitude increase in the nonlinearity coefficient with reduced frequency is obviously absent in the two PZT variants.
A direct comparison between the nonlinearity of BFO and PZT is shown in Figure 2d. As a quantitative measure, we use the fraction of the total piezoelectric coefficient that is due to irreversible DW displacements (X IR ). This parameter is essentially represented by the fraction of the field-dependent coefficient and can be derived from Eq. (1) as: Þis the coefficient at a given field amplitude E 0 and at zero field amplitude, respectively. As anticipated, the irreversible contribution in hard and soft PZT is practically independent on the frequency (Figure 2d, blue and red data). X IR reaches a maximum of 10% and 55% in hard and soft PZT, respectively. Interestingly, the irreversible contribution in BFO at 100 Hz is smaller than that of hard PZT and becomes higher at 0.2 Hz than that of soft PZT (for the same relative field E/E C ; see green data in Figure 2d). Note that the large contribution from the displacements of non-180°domain walls at sub-Hz driving frequencies in BFO was recently confirmed by in situ XRD stroboscopic analysis [84]. One could imagine BFO as very hard material at high frequencies and very soft at low frequencies. The data, therefore, present a distinct hard-to-soft transition in BFO induced by the driving frequency. As will be explained throughout the rest of the chapter, this transition originates from the presence and dynamics of conductive DWs.
Another distinct feature of the piezoelectric response of BFO is the negative piezoelectric phase angle (here and in subsequent discussion, the phase is represented as the tangent of the piezoelectric phase angle, tanδ p ). This rather unusual response is particularly strong in coarse-grained BFO and is presented in Figure 3. From these data it is clear that the piezoelectric response of BFO is rather complex and show strong frequency dependence (Figure 3a). When inspected from high to low frequencies (right to left), the response can be described as a sequence of increasing and decreasing d 33 in the frequency range 1-200 Hz and 0.1-1 Hz, respectively, giving rise to a d 33 maximum at 1.5 Hz (Figure 3a, black points). These two d 33 -vs-frequency behaviors are accompanied by a broad +tanδ p maximum at $10 Hz and a rather narrow -tanδ p minimum at 0.4 Hz, respectively (Figure 3a, red points). The overall behavior can be thus simplified to a sequence of retardation process (d 33 increasing with decreasing frequency accompanied by +tanδ p peak) and relaxation process (d 33 decreasing with decreasing frequency accompanied by -tanδ p peak) [4,46]. As it will be shown in the next section, this is consistent with the distinct features of M-W piezoelectric effect, which can be easily predicted by simple analytical modeling.
As it has been done in the case described here, the sign of the piezoelectric phase is very often determined by measurements performed using lock-in technique, which is capable of extracting the amplitude and phase of individual harmonic responses to external sinusoidal excitations [29]. For less experienced, this may be sometimes non-trivial as several instruments may reverse (by 180°) or somewhat affect the phase of the output signal. While the sign of tanδ p can be checked by, e.g., measurements of a sample with known response, such as, e.g., a donor-doped soft PZT where tanδ p should be positive (and, ideally, related to Eqs. (1) and (2)), it is useful to acquire the total signal containing all information by conventional oscilloscope imaging. Two examples of such acquired signals displayed either in time domain or as piezoelectric hysteresis loop (mechanical displacement ΔL versus electric field E plots) are shown in Figure 3b and c. The case in Figure 3b is shown for 10 Hz where a positive tanδ p was measured using the lock-in method. As expected, the output displacement ΔL signal lags behind the input driving field E signal (see arrow in Figure 3b), corresponding to counter-clockwise rotational sense of the piezoelectric hysteresis. In terms of hysteresis, this is a common situation observed during, e.g., conventional measurements of ferroelectric P-E hysteresis loop at switching fields (to give a popular example). The case that is less common is when tanδ p is negative as shown in Figure 3c. Here, instead of lagging, the output ΔL signal leads the input E signal, effectively resulting in a clock-wise rotation of the piezoelectric hysteresis.
Very often (and not to be blamed) the negative piezoelectric phase is misinterpreted because it gives a wrong impression that it violates the basic law of energy conservation. This is also the reason why is so often interpreted as a measurement artifact. An example is the clockwise hysteresis measured in the response  33 ) and tangent of the piezoelectric phase angle (tanδ p ) as a function of driving field frequency for coarse-grained BFO. The data were obtained at 11.9 kV/cm of driving field amplitude. Error bars represent measurement error. These errors alongside the details of the annealing procedure of the coarse-grained BFO sample (average grain size 16 μm), its microstructure, domain structure and local electrical conductivity are reported in Ref. [79].

The two insets show the time-domain driving field E signal (black sinusoidal curves) with overlaid mechanical displacement ΔL signals (blue or green sinusoidal curves) together with the corresponding ΔL-E piezoelectric hysteresis loops for the case of (b) positive (10 Hz; blue data) and (c) negative (1 Hz; green data) piezoelectric phase angle. The bigger arrows in panels (b) and (c) indicate
the lagging and leading ΔL signal with respect to the driving field signal for the case of positive and negative phase angle, respectively. The corresponding counter-clockwise and clock-wise rotational sense of the hysteresis, which is set by the lagging and leading output signals, respectively, is also noted on the respective loops. Numbered points on the time-domain plots and hysteresis are drawn to help identifying the rotational sense of the hysteresis. of ferroelectric field-effect transistor thin-film structures, which arise due to charge injection during measurements [86,87].
It is not within the scope of this contribution to provide a deep and detailed physical analysis of the negative piezoelectric phase, neither is such analysis within the main expertise of the author of this chapter. The reader is strongly advised to follow the discussion provided in the chapter on hysteresis by Damjanovic [4]. Nevertheless, one can envision a very simple reasoning based on the classical power dissipation principles in dielectrics (found in general textbooks, such as [88]) where the power loss is defined by the positive dielectric tanδ (considering that all other parameters in the equation for the dissipated power density, i.e., electric-field amplitude, frequency and real part of the dielectric permittivity, are positive by definition). Similarly, the area of the charge density (D)electric field (E) hysteresis represents the energy loss of a dielectric during an electric field cycle per unit volume of the material (units of D-E hysteresis is J m À3 ) [89]. In this analogy, one could understand a negative piezoelectric phase angle corresponding to energy gain that is represented by the piezoelectric hysteresis area. The simplest and perhaps strongest argument against this claim is that, unlike the dielectric D-E hysteresis area, the piezoelectric hysteresis area (either converse x-E or direct D-Π) does not have units of energy density. Therefore, the hysteresis shown in Figure 3c, where the phase is negative, does not directly represent an energy gain because the hysteresis area does not reflect energy density. This problem has been extensively elaborated by Holland [90]. The rigorous mathematical treatment therein shows that the necessary restriction of positive total power loss, which is represented by the sum of dielectric, elastic and mix piezoelectric components (proportional to respective imaginary coefficients), results into dielectric and elastic loss terms being always positive and bigger than the piezoelectric term. The latter was shown to not be restricted in its sign. Therefore, taking the longitudinal piezoelectric response of poled ferroelectric ceramics as an example, the imaginary piezoelectric coefficient d 33 " is permitted to be either positive or negative; for a positive real longitudinal coefficient d 33 ', as is the case of ferroelectric ceramics, the piezoelectric tanδ p (defined as tan δ p ¼ d 00

33
) can thus be either positive or negative. While not restricted in its sign, the piezoelectric term in the total power dissipation function is, however, restricted in magnitude. This essentially means that a negative piezoelectric phase angle measured in poled ceramics indicates a partial reduction of the total power loss. To avoid confusion, it should be noticed that in the power dissipation equation, reported in the paper by Holland, the term related to the piezoelectric coupling is proportional to EÁΠÁd", which must obviously possess units of energy density (subscripts are omitted for simplicity). In contrast, the piezoelectric hysteresis area is proportional to either the product of xÁE (converse effect) or DÁΠ (direct effect), none of which possess energy density units. A very nice experimental example of the reduction of the total power dissipation due to piezoelectric coupling was shown for Sm-doped PbTiO 3 [4].
Coming back to our case on BFO, it must be emphasized at this point that the negative converse piezoelectric phase angle has also been confirmed in the direct piezoelectric response of BFO [91]. In addition, a negative phase has also been detected in the lattice microstrain response of BFO ceramics to external electric field, which was characterized by in situ XRD stroboscopic analysis [84]. This is perhaps the most important evidence of a negative piezoelectric phase angle measured in any piezoelectric so far because, prior to that work, such response (to the best of the author's knowledge) was demonstrated only on the level of macroscopic measurements.
As anticipated in the initial discussion of the results shown in Figure 3, the negative piezoelectric phase angle is a strong indication of M-W piezoelectric process that is very common in, e.g., piezoelectric composites [46,47]. In analogy to the dielectric M-W relaxation, characterized by giant apparent dielectric permittivity [50], the modeling of Turik et al. showed the same effects in the piezoelectric M-W analogue [45]. As explained earlier, the M-W relaxation has origin in the internal electric field redistribution in a medium, where local regions exhibit different electrical conductivities; in the frame of modeling, such regions are often implemented in the form of M-W bilayer units, where the two layers are described by different conductivities [48,51]. This is the reason behind a large piezoelectric M-W effect in some Aurivillius phases, such as Bi 4 Ti 3 O 12 [48]. These ceramics, in fact, tend to show anisotropic microstructure with elongated plate-like grains that are characterized by different electrical conductivities in the direction parallel or perpendicular to the plane of the plates. The piezoelectric M-W effect in heterogeneous Bi 4 Ti 3 O 12 ceramics is thus easy to support using arguments of anisotropic conductivity.
Unlike in Bi 4 Ti 3 O 12 , the M-W features observed in the piezoelectric response of BFO (Figure 3) are more difficult to interpret because significant anisotropy in the conductivity in a homogeneous perovskite oxide, such as BFO, is not expected, at least not to a level as in layered Aurivillius-type structures. Also, BFO is characterized by a microstructure typically consisting of equiaxial grains (as it is illustrated in Figure 4 in the next section). It should be recalled, however, that the piezoelectric M-W effect in BFO is particular in that it shows, in addition to the negative piezoelectric phase angle (Figure 3), a very strong DW contribution observed only at low (sub-Hz) driving frequencies (Figure 2c). The overall data thus suggest a piezoelectric M-W effect provoked by DWs. The idea becomes reasonable when DWs are electrically conducting as the conductivity is what triggers the M-W effect. As it will be shown in the next section, it is exactly the conducting DWs that likely cause large anisotropy in the conductivity from grain to grain or across cluster of grains. Obviously, in BFO the situation is more complex than in other known piezoelectric M-W cases, simply because the features that are triggering the M-W effects, i.e., the DWs, can also move inside the grains under applied fields.

Analytical modeling of nonlinear Maxwell-Wagner piezoelectric relaxation arising from electrically conducting DWs
While, pragmatically speaking, the idea of conducting DWs giving rise to macroscopic piezoelectric M-W effects does fit into a reasonable explanation, the details of the mechanism are rather complicated. The reason is the complex situation in a polycrystalline ceramic matrix containing randomly oriented grains, each characterized by a domain structure with different DWs and corresponding DW planes oriented in various directions with respect to the reference external field axis. The DW planes are assumed to provide conductive paths through individual grains or in local regions inside the grains, modifying the internal electric fields. If further grain-to-grain elastic interactions are considered, which have a significant effect in polycrystalline piezoceramics [26], modeling may become extremely difficult. In the first approximation, therefore, the model can be simplified and reduced to a level that allows the major macroscopic parameters to be predicted and compared with experimental data. As demonstrated for the case of direct piezoelectric response of Aurivillius phases [48], analytical bilayer modeling is sometimes sufficient. A similar case study on BFO, but additionally elaborated to incorporate nonlinear effects, will be presented in this section.
In the first place, it is important to work out a physical picture onto which to construct the model. Following the original discovery of conducting DWs in epitaxial BFO thin films [92], conduction at DWs in BFO polycrystalline samples have been determined using the same c-AFM measurements, which make it possible to probe the local conductivity of the sample. In BFO films, the epitaxial growth determines the crystallographic orientation of the film with respect to the substrate plane, so that the determination of the type of DWs in rhombohedral (R) BFO (71°, 109°or 180°) is facilitated. In ceramics, this is more difficult as the orientation of the grains is random and it is thus not known when imaging the surface with microscopy techniques. A simple way to accomplish the task is to determine the orientation of the polar [111] axis of the BFO R symmetry in adjacent domains in a selected grain using, e.g., electron back-scattered diffraction (EBSD) analysis. By indexing the Kikuchi patterns obtained in individual domains, as those highlighted in blue circles in Figure 4a, with the available R BFO space group, the angle between the [111] vectors in adjacent domains can be determined. Once these angles are known and without going into detailed geometrical analyses, which may be non-trivial (see the example in the supplementary material 1 of the paper [93]), a separation between 180°and non-180°DW-type is possible. In the former case, the angle should be zero ([111] vectors in adjacent domains are parallel), while in the second is non-zero. Note that further analysis is complicated by the fact that EBDS cannot determine the orientation of the ferroelectric spontaneous polarization, but only the ferroelastic distortion.
EBSD analysis was used to determine the 180°and non-180°DWs in the example shown in Figure 4a (see also reference [79] for further details). In the next step, these same regions were analyzed by c-AFM (Figure 4b) to finally confirm that both 180°and non-180°DWs in R BFO ceramics exhibit enhanced electrical conductivity with respect to that measured in the domains (see also the electric-current profiles in Figure 4b measured across DWs as indicated with the dashed white lines in respective c-AFM maps). The results are fully consistent with the atomic-scale microscopy analysis, which confirmed the presence of p-type charge carriers, identified as Fe 4+ states, concentrated inside all three DW variants (71°, 109°and 180°) of BFO [19], resulting in the domains mostly likely depleted from the charges. The p-type conductive nature of the DWs is supported by annealing studies in controlled atmospheres, reported in the same paper, while the dynamic pinning effects of the p-type carriers is reported in Ref. [75].
The analysis of the effect of conductive DWs on the piezoelectric response can be reduced, in the first approximation, to a two-grain problem as shown in Figure 5a (left schematic). Since ceramics are composed of randomly oriented grains containing different type of DWs with various orientations, it is legitimate to analyze a hypothetical case of two grains oriented along [100] and [111] directions with respect to the externally applied electric field vector. The choice of these orientations will become evident in the subsequent discussion. If the two grains contain only 71°DWs percolating the grains, then the situation should closely resemble the left schematic shown in Figure 5a (detailed geometrical analysis of the DW angles in the two grains is reported elsewhere [84]). The key element of the model is the different orientation of conductive DWs in the two grains with respect to the external field axis, i.e., in the top grain 1 (red) the DWs are parallel to E, while in the bottom grain 2 (blue) they form a different angle. The net result is that the conductivity, measured along E, of grain 1 should be higher than that of grain 2 because in the former the charges may migrate along the vertical conductive DWs (see notations σ 1 and σ 2 in Figure 5a indicating the conductivity measured vertically in grain 1 and 2, respectively). Using the same reasoning, the conductivity of grain 2 should be higher when measured along a direction away from the E axis (not parallel). It is this anisotropy in the electric conductivity that it is assumed here to lead to M-W-like internal field redistribution (see down-arrowed E 1 and up-arrowed E 2 notations in Figure 5a, denoting the internal fields). Note also that the bottom grain 2 is oriented in a way to give rise to a stronger DW contribution than the upper grain 1; in the ideal case, due to orientational constraints, grain 1 should exhibit only lattice strain as a response to the external E (in Figure 5a the microstrains arising from the two grains due to field application are noted as ΔL 1 and ΔL 2 ). Further details regarding the model and mathematical derivations can be found in the supplementary Section 5 of the paper by Makarovic et al. [94].
The assembly of the two grains represents a basic M-W bilayer unit. As in classical dielectric M-W modeling [50], the picture can be rationalized in terms of an equivalent circuit consisting of two leaky capacitors connected in series, where each leaky capacitor may be represented by an ideal capacitor and an ideal resistor connected in parallel (Figure 5a, right-hand schematic). Obviously, the piezoelectric effect is added to the two layers (grains). Due to the expected DW contribution in the bottom grain 2, the RL relations (Eqs. (1) and (2)) were used to calculate the piezoelectric response of this grain. In contrast, grain 1 is assumed to respond via intrinsic lattice piezoelectric effect, so its response can be modeled by the linear constitutive piezoelectric equation. For a set of dielectric, piezoelectric and conductivity parameters, which are within the margins reasonable for BFO (see reference [94]), the results of the modeling are shown in Figure 5b Driven by the anisotropy in the electrical conductivity, defined as the conductivity ratio of the two grains ( σ 1 σ 2 ), the internal fields in the grains (E 1 , E 2 ) will be redistributed as a function of driving field frequency (Figure 5b). In accordance to the conductive behavior related to the different DW orientation in the two grains, the sinusoidal electric field applied externally to the bilayer serial structure will be redistributed at low driving frequencies (<10 Hz) such that the electric field in the top grain 1 (red color coding) will be reduced with respect to the nominal (external)  [100] pc and [111] pc oriented with respect to field axis E and are denoted as grain 1 (red) and grain 2 (blue), respectively. Notation "pc" indicates "pseudocubic". The dashed lines inside the grains represent DWs. The DW orientation is assumed such that the DW planes in grain 1 are parallel to E, while those in grain 2 are not, forming an angle of $35°with respect to E axis. These angles should correspond to those of 71°DWs present in the two grains (for detail geometrical analysis, see reference [84]). Due to the different orientation of the electrically conductive DW in the two grains, the conductivity measured along E axis will be higher in grain 1 than in grain 2 (see σ 1 and σ 2 notation with upward and downward arrow in grain 1 and grain 2, respectively). This will result in the drop of the electric field in grain 1 (noted by down-arrowed E 1 ) and rise of the electric field in grain 2 (up-arrowed E 2 ). The same difference in the orientation of the DW with respect to E will result in a lattice strain with no contribution from DW motion in grain 1 (see ΔL 1 ), while displacements of DWs is expected in grain2, in addition to lattice strain (see larger ΔL 2 ). To incorporate DW contributions, the calculations were performed assuming RL relations in grain 2 as described in detail in Ref. [94]. The two grains may be viewed as a M-W bilayer unit where each layer is assumed as a leaky piezoelectric element, characterized by its own piezoelectric coefficient d i , dielectric permittivity ε i , specific DC electrical conductivity σ i and volume fraction ν i (i = 100 or 111 indicating the two grain orientations field, while that of the bottom grain 2 (blue) will be increased instead (Figure 5b). In other words, due to leakage in grain 1 caused by the vertical orientation of conductive DWs, the field inside this grain will drop and will be thus transferred from the leaky grain 1 to the less-leaky grain 2 (see Figure 5a). This will happen at low driving frequencies as such driving conditions provide to the charges sufficient time to migrate along conductive DWs. Being proportional to the internal field, the microstrains of the two grains (Figure 5c) will show the same frequency behavior as that of respective internal fields, leading to either retardation and a peak in the positive piezoelectric phase angle (see tanδ 2 in Figure 5d), or relaxation and a peak in the negative piezoelectric phase angle (see tanδ 1 in Figure 5d). This M-W mechanism will ultimately lead to nonlinear effects: the increased internal electric field in grain 2 (E 2 in Figure 5b) will boost the DW motion in this grain, resulting in effective nonlinearity enhanced at low driving frequencies as shown in Figure 5e (compare also with experimental data, Figure 2c). The piezoelectric M-W effect has thus a nonlinear character or, in other words, the piezoelectric nonlinearity (i.e., DW motion) is boosted due to the electrical conductivity.
At this stage, I have to point out that the presented model used to understand the experimental data on BFO (as those shown in Figures 2c and 3) is largely simplified: (i) it considers only two isolated and unconstrained grains connected in series, neglecting the true elastic and electric boundary conditions of the analyzed grains set by the presence of other surrounding grains, (ii) it does not consider elastic coupling of the two isolated grains and transverse piezoelectric effects (as was done in, e.g., reference [51]), (iii) it assumes very simple domain structure in the two grains and only one type of DW (i.e., 71°), and (iv) the conductivity of the grains set by the conductive paths along DWs is assumed to be fixed, although a dynamic conductivity is not unreasonable considering that the conductive DWs may switch locally (through a RL-like irreversible jump, for example), resulting in a modified conductive path through the switched DW. Despite all these limitations, however, it is surprising to realize that the simple two-grain model predicts all the key experimental observations. First, the model predicts both retardation and relaxation processes (Figure 5c and d), which are likely convoluted in the experimental response (as discussed for Figure 3a). Second, the model shows that the lattice strain response should exhibit a negative phase with respect to sinusoidal external field, thus leading the field signal (Figure 5d). Not only is this consistent with macroscopically measured piezoelectric strain leading the field signal (Figure 3c), but this is also exactly the behavior of the lattice strain deconvoluted from the total converse piezoelectric response of BFO ceramics using synchrotron XRD measurements [84]. Third, the important outcome of the model related to the low-frequency nonlinearity is a natural consequence of introducing RL relations, lining up with the macroscopic experimental data (Figure 2c) and the low-frequency DW contribution determined by synchrotron XRD measurements [84].
The low-frequency nonlinearity has recently been shown to be a response parameter that can be controlled by designing the fraction of conductive DWs in BFO ceramics via doping [94]. The idea was triggered by the outcomes of the model itself, further reinforcing the value and importance of simple modeling. It is probably needless to say that controlling nonlinearity and hysteresis is very important for the development of high-temperature piezoceramics based on BFO. This is indeed supported by a study on BFO showing that the temperature dependent piezoelectric response of these ceramics is strictly controlled by the same M-W processes described in this chapter [95]. Importantly, it was shown that the strong temperature dependence of the piezoelectric nonlinearity and hysteresis, which is of a direct relevance for the device operation, has origin in the thermally activated nature of the local electrical conductivity in BFO ceramics.
Simple analytical modeling was found to be a promising first step toward understanding the complex piezoelectric behavior of BFO arising from conductive DWs. It could be interesting, however, to model a more complex situation to account for the many different parameters that are necessarily neglected in the simple analytical approach. As shown in Figure 5b-d, the effect of the anisotropic parameter is crucial as it determines the strength of the M-W effect. For more advanced modeling, this parameter could be viewed as varying from one M-W unit to another in a complex ceramic matrix composed of interacting grains forming units with different time constants (τ is proportional to the ratio of the weighted sum of permittivity and conductivity of the individual layers in the bilayer unit [48,84]). This could eventually account for the different regions inside the ceramics exhibiting different levels of electric-field redistribution (as shown by the different curves in Figure 5b). 3-dimensional (3D) finite element modeling based on a phenomenological approach has been recently demonstrated to be a powerful tool in predicting local electric fields in 3D ceramic matrices with defined porosity [96]. For solving complex problems, such as those encountered in elastically and electrically coupled grains in ferroelectric ceramics, multiscale modeling approaches show great promises [97][98][99].
Another interesting point that could be considered in the future is to push the relaxation to higher frequencies (equal to decreasing the τ value). This would make it possible to use the large response in a higher frequency range that is more relevant for piezoelectric applications (in the example shown in Figure 5, τ = 0.13 s, corresponding to a relaxation frequency of f relax = 1.2 Hz). As discussed in the preceding paragraph, the anisotropy in the conductivity is also an important factor in tailoring the usable frequency range. It is thus tempting to consider designing a matrix with charged DWs, which can possess metallic-like conductivity, as recently demonstrated for BaTiO 3 [100]. In this case, the DW conductivity may exceed the bulk conductivity for impressive 8-10 orders of magnitude, perhaps providing an opportunity in designing piezoelectric properties.

Summary and conclusions
In this chapter, I have reviewed and discussed a case study on BFO ceramics explaining how the presence and dynamics of DWs showing enhanced electrical conductivity with respect to bulk conductivity can have a crucial effect on the macroscopic piezoelectric response of BFO. The mechanism goes far beyond the expected DW-defects pinning interactions, described in Section 2 of this chapter, and reflect itself in the nonlinear and hysteretic piezoelectric response of BFO. The unusual features of the response, i.e., a hard-to-soft transition induced by lowering the driving electric-field frequency and negative piezoelectric phase angle, can be explained by nonlinear piezoelectric Maxwell-Wagner effects. Simple analytical modeling confirms these macroscopic-response features and show that the key entities leading to such effects are conductive DWs. I could envision few points that can be drawn from the presented results. First, the data clearly show that local conductive paths (such as those along DWs) should be considered more seriously in addition to conventional bulk conductivity, which is mostly discussed in the literature on BFO. This is particularly important for the development of next-generation BFO-based piezoceramics for high-temperature applications as the local conductivity is what makes the response of BFO unstable in terms of driving field parameters (amplitude and frequency) and, most importantly, temperature. Second, the key problem related to M-W effects is that the response is boosted only at quasi-static driving conditions, as shown earlier by modeling of ceramic-polymer composites [45]. While certainly limited in the frequency range, it could be interesting to test anisotropic effects and piezoelectric enhancements in engineered matrices containing highly conducting charged DWs. It seems reasonable, though, that one should first validate the idea by modeling.