Ferroelectric Domain Structure and Local ...

9 downloads 18161 Views 12MB Size Report
Jan 7, 2017 - guideline for the understanding of domain processes in lead-free ..... The bi-domain structure was created in a single grain using ±30 V.
materials Review

Ferroelectric Domain Structure and Local Piezoelectric Properties of Lead-Free (Ka0.5Na0.5)NbO3 and BiFeO3-Based Piezoelectric Ceramics Denis Alikin 1 , Anton Turygin 1 , Andrei Kholkin 1,2 and Vladimir Shur 1, * 1 2

*

School of Natural Sciences and Mathematics, Ural Federal University, Ekaterinburg 620000, Russia; [email protected] (D.A.); [email protected] (A.T.); [email protected] (A.K.) Department of Physics, CICECO-Aveiro Institute of Materials, University of Aveiro, Aveiro 3810-193, Portugal Correspondence: [email protected]; Tel.: +7-343-261-7436

Academic Editor: Ingo Dierking Received: 25 November 2016; Accepted: 30 December 2016; Published: 7 January 2017

Abstract: Recent advances in the development of novel methods for the local characterization of ferroelectric domains open up new opportunities not only to image, but also to control and to create desired domain configurations (domain engineering). The morphotropic and polymorphic phase boundaries that are frequently used to increase the electromechanical and dielectric performance of ferroelectric ceramics have a tremendous effect on the domain structure, which can serve as a signature of complex polarization states and link local and macroscopic piezoelectric and dielectric responses. This is especially important for the study of lead-free ferroelectric ceramics, which is currently replacing traditional lead-containing materials, and great efforts are devoted to increasing their performance to match that of lead zirconate titanate (PZT). In this work, we provide a short overview of the recent progress in the imaging of domain structure in two major families of ceramic lead-free systems based on BiFeO3 (BFO) and (Ka0.5 Na0.5 )NbO3 (KNN). This can be used as a guideline for the understanding of domain processes in lead-free piezoelectric ceramics and provide further insight into the mechanisms of structure–property relationship in these technologically important material families. Keywords: KNN; BFO; MPB; PPB; domain structure; piezoelectric properties; local switching

1. Introduction Piezoelectric materials exhibit a unique ability to expand under external electric field or to develop a charge under applied mechanical stress, combining high coupling coefficients with exceptional stability and low cost. They are used in a variety of devices, including resonators, ultrasound generators, and actuators realizing precise nano-motion [1–5]. The most widely used materials are currently lead zirconate titanate (PZT)-based ceramics, which are well-known for their excellent piezoelectric properties. However, considering the toxicity of lead and its derivatives, there is a general trend for the development of environmentally-friendly lead-free materials as regulated by several incentives of the European Union. The international efforts in removing toxic substances from everyday applications have been pursued in the last decade. The EU passed the “Waste Electrical and Electronic Equipment” (WEEE) and “Restriction of the use of certain Hazardous Substances in electrical and electronic equipment” (RoHS) initiatives in 2003 [6]. While the WEEE regulates the disposal, reuse, and recycling of the aforementioned devices, RoHS requires that this can be accomplished safely without endangering the environment or people’s health. As a result, two classes of piezoelectric materials are now being considered as potentially attractive alternatives to PZTs for specific applications: (i) perovskites—i.e., (Bi0.5 Na0.5 )TiO3 Materials 2017, 10, 47; doi:10.3390/ma10010047

www.mdpi.com/journal/materials

Materials 2017, 10, 47

2 of 23

(BNT), BaTiO3 (BT), BiFeO3 (BFO), KNbO3 , NaNbO3 , and solid-solutions on their base; (ii) non-perovskites—i.e., bismuth layer structured ferroelectrics and tungsten-bronze-type compounds. However, their piezoelectric properties were found to be inferior compared to those of PZT. It was shown that—similarly to PZT—the properties could be improved by forming solid solutions of KNbO3 with NaNbO3 (KNN), so as to sinter ceramics with a composition close to the region of structural instability known as the polymorphic phase boundary (PPB) regulated by the temperature and the morphotropic phase boundary (MPB) controlled by doping. Similar to PZT, the properties of ceramics are also governed by the domain structure and the grain morphology, which are in turn determined by the defect transport and controlled by the sintering method. Such complex interplay of the physical and chemical properties of lead-free piezoelectric ceramics has been recently rationalized in terms of domain wall conductivity and the diffusion of charged defects [7]. Taking recent progress in domain visualization [8,9] and structure refining into account, we believe that the review on the domain-related properties of lead-free piezoelectric ceramics such as BFO and KNN solid solutions is quite timely and can be useful for the specialists working in this field. We will start with a short description of the sintering methods and phase diagrams of these ceramics, and continue with the review of the modern domain visualization techniques. In the following, we will overview the current status of domain studies in BFO and KNN ceramics, and will end with an analysis of the effect of domain structure on the physical properties of ceramics, such as switching behavior and dielectric constant. 2. Sintering and Crystal Structure 2.1. BiFeO3 System Bismuth ferrite oxide (BFO) is a well-known room temperature multiferroic material belonging to the space group R3c [10]. It is attractive due to its very high value of spontaneous polarization of about 100 µC/cm2 [10]. However, BFO prepared by a solid state synthesis is typically characterized by a high leakage current due to the high concentration of Fe2+ ions and oxygen vacancies [11]. Therefore, it is difficult to polarize the ceramics in order to obtain sufficiently high piezoelectric and dielectric properties. The phase diagram Bi2 O3 –Fe2 O3 shows the existence of two equilibrium phases (Fe and Bi saturated): orthorhombic Bi2 Fe4 O9 (Pbam), rhombohedral perovskite BiFeO3 (R3c), which decomposes to Bi2 Fe4 O9 and liquid phase at 935 ◦ C, and cubic Bi25 FeO39 , which decomposes to Bi2 O3 and liquid phase at 790 ◦ C [12] (Figure 1). A number of experimental and theoretical works have shown the temperature instability of BFO at T > 700 ◦ C with decomposition to Bi2 O3 , Fe2 O3 , and Bi2 Fe4 O9 [13]. It has been shown that in the temperature range from 447–767 ◦ C, Bi25 FeO39 and Bi2 Fe4 O9 phases are more thermodynamically stable than BiFeO3 , and three phases can coexist [11]. During ceramics synthesis in the temperature range 447–767 ◦ C, BiFeO3 partially converts into other phases [14]. The reaction is reversible, and additional annealing returns a BiFeO3 phase. A multi-phase state is detected experimentally in a much wider temperature range that can be explained by the small change of Gibbs potential being the driving force for the reaction [15]. Moreover, it has been shown by Valant et al. [16] that the multi-phase ceramic state may be a result of interaction with different types of impurity oxides (AOx ). Different doping strategies (alkaline earth ions, rare earth ions, and others) have been used for the formation of MPBs resulting in the reduction of leakage currents, and therefore, improvement of the ferroelectric and piezoelectric properties [17–22].

Materials 2017, 10, 47 Materials 2017, 10, 47

3 of 23 3 of 23

Figure 1. 1. Phase diagram of of the BiBi 2O23O –Fe 2O32 O system. Open circles show thethe data points obtained byby Figure Phase diagram the Open circles show data points obtained 3 –Fe 3 system. differential thermal analysis (DTA). The dotted line above thethe liquidus represents thethe approximate differential thermal analysis (DTA). The dotted line above liquidus represents approximate temperature limit notnot to to be be surpassed in order to avoid decomposition, otherwise correct equilibrium temperature limit surpassed in order to avoid decomposition, otherwise correct equilibrium DTA peaks areare nono longer observed upon a second heating. The α, α, β, and γ phases areare rhombohedral, DTA peaks longer observed upon a second heating. The β, and γ phases rhombohedral, orthorhombic, and cubic, respectively. Adapted from [23], with permission from © 2008 American orthorhombic, and cubic, respectively. Adapted from [23], with permission from © 2008 American Physical Society. Physical Society.

2.2. KNbO 3–NaNbO 3 System 2.2. KNbO 3 –NaNbO 3 System (K,(K, Na)NbO 3 is a solution of of ferroelectric KNbO 3 and antiferroelectric NaNbO 3, exhibiting Na)NbO a solid solution ferroelectric KNbO antiferroelectric NaNbO 3 issolid 3 and 3 , exhibiting 2), and 2 ◦ sufficiently high Curie temperature (T c = 420 °C), good ferroelectric properties (P r = 33 μC/cm sufficiently high Curie temperature (Tc = 420 C), good ferroelectric properties (Pr = 33 µC/cm ), large coupling factorsfactors [24]. The diagram for theforKNN system is shown in andelectromechanical large electromechanical coupling [24].phase The phase diagram the KNN system is shown Figure 2. At 2. room MPBs lie at 17.5%, 47.5% contents. It is notable in Figure At temperature, room temperature, MPBs lie at 32.5%, 17.5%, and 32.5%, andNaNbO 47.5% 3NaNbO It is 3 contents. ◦ that phase transition temperatures between ferroelectric phases at ~200 °C and between ferroelectric notable that phase transition temperatures between ferroelectric phases at ~200 C and between ◦ Cindependent and paraelectric at ~400 °C are almost of the composition (in contrast the ferroelectric andphases paraelectric phases at ~400 are almost independent of the composition (in to contrast composition-dependent transitiontransition temperature of PZT).ofOnly of sodium for to the composition-dependent temperature PZT).small Only substitutions small substitutions of sodium potassium in NaNbO 3 cause a transition to ferroelectric from pure antiferroelectric sodium niobate for potassium in NaNbO3 cause a transition to ferroelectric from pure antiferroelectric sodium [25,26]. The piezoelectric data for the air-fired around 33 = 80 pC/N, and density the niobate [25,26]. The piezoelectric data for the samples air-fired are samples are daround d33 = 80 pC/N, andof density 3 [24,27]. One 3 sample is around 4.25 g/cm of the main obstacles for the development of potassium of the sample is around 4.25 g/cm [24,27]. One of the main obstacles for the development of sodium niobate solid solution (KNN) as a commercial ceramic material by conventional potassium sodium niobate solid solution (KNN) aspiezoelectric a commercial piezoelectric ceramic material by method is the method difficulty in difficulty processing and densification. Egerton and co-workers reported the conventional is the in processing and densification. Egerton and co-workers reported electrical properties of KNN, in which they indicated relatively low dielectric constants over a wide the electrical properties of KNN, in which they indicated relatively low dielectric constants over a compositional range [24]. Hence, achieve densification, hot-pressed KNN ceramics (~99% wide compositional range [24].to Hence, to sufficient achieve sufficient densification, hot-pressed KNN ceramics of (~99% the theoretical density)density) have been reported to possess a higha Curie temperature (Tc =(T 420 of the theoretical have been reported to possess high Curie temperature = 420 ◦aC), c °C), large piezoelectric longitudinal response (d 33 = 160 pC/N), and a high planar coupling coefficient a large piezoelectric longitudinal response (d33 = 160 pC/N), and a high planar coupling coefficient (kp(k=p 45%). KNN samples have been prepared byby conventional airair sintering in in order to to reach high = 45%). KNN samples have been prepared conventional sintering order reach high densities (over 95%), which yielded superior piezoelectric properties (d 33 = 100 pC/N) as compared densities (over 95%), which yielded superior piezoelectric properties (d33 = 100 pC/N) as compared to to those obtained by by thethe same method [28].[28]. It is Itimportant to note that KNN material prepared by those obtained same method is important to note that KNN material prepared spark plasma sintering showed significantly higher higher dielectric and piezoelectric properties than those by spark plasma sintering showed significantly dielectric and piezoelectric properties than prepared by conventional method (ε ~ 700 and d 33 ~ 148 pC/N) [29,30]. Saito et al. fabricated texturedthose prepared by conventional method (ε ~700 and d33 ~148 pC/N) [29,30]. Saito et al. fabricated based KNN ceramics byceramics the reactive grain-growth method, which resulted d33 value as as textured-based KNN by the reactive grain-growth method, whichinresulted in as d33high value ~416 pC/N [31].pC/N [31]. high as ~416

Materials 2017, 10, 47 Materials 2017, 10, 47

4 of 23 4 of 23

Figure 2. 2. Phase Phase diagram diagramfor forthe thesystem systemKNbO KNbO3–NaNbO –NaNbO3.. Adapted Adapted from from [32], [32], with with permission permission from from Figure 3 3 © 2009 AIP Publishing LLC. © 2009 AIP Publishing LLC.

Furthermore, the volatility of potassium oxide makes it difficult to maintain stoichiometry and Furthermore, the volatility of potassium oxide makes it difficult to maintain stoichiometry and high density of ceramics [33]. In order to optimize the processing conditions and to obtain high density of ceramics [33]. In order to optimize the processing conditions and+to obtain reproducible reproducible properties, KNN ceramics are doped with suitable elements: Li and alkali-earth ions, properties, KNN ceramics are doped with suitable elements: Li+ and alkali-earth ions, especially 2+. The doping changes the cell parameters, promotes densification, decreases the phase especially Sr Sr2+ . The doping changes the cell parameters, promotes densification, decreases the phase transition transition temperatures, and improves the electrical properties [34–37]. The addition of LiSbO3 or temperatures, and improves the electrical properties [34–37]. The addition of LiSbO3 or LiTaO3 to LiTaO3 to KNN leads to sufficient enhancement of the dielectric, piezoelectric, and ferroelectric KNN leads to sufficient enhancement of the dielectric, piezoelectric, and ferroelectric properties of properties of ceramics [38,39]. For more details, we refer the reader to the following exhaustive ceramics [38,39]. For more details, we refer the reader to the following exhaustive reviews [33,40,41]. reviews [33,40,41]. 3. Methods of Domain Structure Visualization 3. Methods of Domain Structure Visualization Rapid development of microscopy techniques brought a lot of possibilities for domain observation Rapid development of The microscopy techniques brought a lot ofmethods, possibilities for domain in ferroelectric materials [8]. domains can be visualized by several including optical observation in ferroelectric materials [8]. The domains can be visualized by several methods, microscopy, scanning electron microscopy (SEM) with electron backscatter diffraction (EBSD) [42], including optical microscopy, scanning microscopy (SEM) with electron Transmission electron microscopy (TEM) electron [43], X-ray diffraction [44,45], and variousbackscatter scanning diffraction (EBSD) [42], Transmission electron microscopy (TEM) [43], X-ray diffraction and probe microscopy (SPM) techniques [46–49], such as piezoresponse force microscopy [44,45], (PFM) [50], various Raman scanning probe microscopy techniques [46–49], such as [48]. piezoresponse force confocal microscopy (CRM) [51],(SPM) and electric force microscopy (EFM) The most useful microscopy (PFM) [50], confocal Raman microscopy (CRM) [51], and electric force microscopy (EFM) methods for domain visualization in ceramics with high spatial resolution are SEM after selective [48]. The etching, most useful methods forCRM. domain visualization in ceramics with high spatial resolution are chemical TEM, PFM, and SEM after selective chemical etching, TEM, PFM, and CRM. 3.1. Scanning Electron Microscopy after Selective Chemical Etching 3.1. Scanning Electron Microscopy after Selective Chemical Etching Surface chemical etching was the earliest method for the visualization of the static domain Surface chemical etching was the and earliest method for the of theisstatic structure in ferroelectric single crystals ceramics [52,53]. Thisvisualization surface treatment baseddomain on the structure in ferroelectric single crystals and ceramics [52,53]. This surface treatment is based on the different etching rates of the opposite polarities of a polarization dipole [52,54]. Typically, several acids different the opposite of aused polarization dipole [52,54].depending Typically, on several (HF, HCl,etching HNO3 )rates and ofalkalis (NaOH,polarities KOH) are as etchant solutions, the acids (HF, HCl, HNO 3 ) and alkalis (NaOH, KOH) are used as etchant solutions, depending on the material [55–57]. Domains were visualized in BFO after etching with HNO3 at room temperature for material [55–57]. Domains were visualized in BFO after etching with HNO 3 at room temperature for 2–4 h [55] or with 0.5% HF at room temperature for 45 s [58]. KNN ceramics were chemically etched 2–4 h [55] with 0.5% HF attemperature room temperature 45 s [58]. ceramics were chemically in 48% HF or solution at room for five for minutes [59].KNN The resulting nm-scale change etched of the in 48% HF solution at room temperature for five minutes [59]. The resulting nm-scale change of the the surface relief can be visualized with different microscopy methods [60–62]. The spatial resolution of surface technique relief can depends be visualized different method microscopy [60–62]. of etching on thewith visualization and methods is typically below 2The nmspatial in SEMresolution registration the etching technique depends on the visualization method and is typically below 2 nm in SEM scattering of electrons from the surface relief (Figure 3a–f). The disadvantages of this method are its registrationinfluence scatteringonofthe electrons theand surface relief (Figure 3a–f). The disadvantages of this destructive sample from surface possible partial back-switching of domain structure method are its destructive on theused sample surface and possible back-switching of upon etching [63]. However,influence it is frequently together with EBSD to getpartial information on domain domain structure upon etching [63]. However, it is frequently used together with EBSD to get orientations in randomly-oriented ceramic grains [64–66]. information on domain orientations in randomly-oriented ceramic grains [64–66].

Materials 2017, 10, 47 Materials 2017, 10, 47

5 of 23 5 of 23

Figure 3. SEM images of domain patterns for the poled NaNbO3 (KNN) ceramics. The magnifications Figure 3. SEM images of domain patterns for the poled NaNbO3 (KNN) ceramics. The magnifications are (a,b) ×3000; (c,e) ×8000; (d) ×10,000 respectively; (f) The partially enlarged view of SEM image. are (a,b) ×3000; (c,e) ×8000; (d) ×10,000 respectively; (f) The partially enlarged view of SEM image. Models of domain configurations with the front planes as the observation planes: (g) Model for Models of domain configurations with the front planes as the observation planes: (g) Model for domain configuration type I; and (h) for domain configuration type II. The grains with domain domain configuration type I; and (h) for domain configuration type II. The grains with domain configuration type I are in (c,d), grains with domain configuration type II are in (e,f). Adapted from configuration type I are in (c,d), grains with domain configuration type II are in (e,f). Adapted from [67], [67], with permission from © 2013 AIP Publishing LLC. with permission from © 2013 AIP Publishing LLC.

3.2. Transmission Electron Microscopy 3.2. Transmission Electron Microscopy TEM has a very high spatial resolution (below 1 nm), but requires very thin samples; thus, it is TEM has a very spatial resolution (below 1 nm), requires very thin samples; thus, it is very sensitive to thehigh sample preparation [68]. Domains can but be irreversibly modified upon polishing andsensitive focused ion beam etching. TEM can be used for the imaging of domains, modified domain walls, local very to the sample preparation [68]. Domains can be irreversibly uponand polishing phase distribution in different piezoelectric ceramics, including BFO and KNN [69,70]. Contrast in and focused ion beam etching. TEM can be used for the imaging of domains, domain walls, and local this method is provided by different mechanisms, such as the scattering of high energy electrons phase distribution in different piezoelectric ceramics, including BFO and KNN [69,70]. Contrast ininthis the local electric and stress fields [71]. Scanning transmission electron microscopy (STEM)ininthe method is provided by different mechanisms, such as the scattering of high energy electrons aberration corrected mode coupled with electron-energy loss spectroscopy can be used for direct local electric and stress fields [71]. Scanning transmission electron microscopy (STEM) inthe aberration measurement of the atomic displacement in domain wall regions and investigation of defect structure corrected mode coupled with electron-energy loss spectroscopy can be used for the direct measurement local strains [7,72,73]. Itinisdomain of particular importance ceramics, because the structure relation between of and the atomic displacement wall regions andfor investigation of defect and local defects (such as dislocations) and domain configurations can be easily identified. Defects and strains [7,72,73]. It is of particular importance for ceramics, because the relation between defects disorder at the grain boundaries can be also seen by TEM. Modern TEM microscopes provide the (such as dislocations) and domain configurations can be easily identified. Defects and disorder at possibility of directly applying an electric field in the microscope camera, and, therefore, to study the grain boundaries can be also seen by TEM. Modern TEM microscopes provide the possibility of domain kinetics in situ [74]. directly applying an electric field in the microscope camera, and, therefore, to study domain kinetics in 3.3. situPiezoresponse [74]. Force Microscopy PFM is one of Microscopy the most useful methods for the visualization of the domain structure in 3.3. Piezoresponse Force ferroelectric materials, due to its simple sample preparation, generally no need of vacuum or other PFM is one of the most useful methods for the visualization of the domain structure in special conditions, high signal-to-noise ratio, nanometer spatial resolution, and variety of different ferroelectric materials, to itsthe simple sample preparation, generally no needand of vacuum spectroscopic modes due allowing measurement of local material ferroelectric dielectric or other special[50,75]. conditions, signal-to-noise ratio, nanometer and variety of properties PFM ishigh a strain-based scanning probe microscopyspatial [76,77]resolution, mode, where application different spectroscopic modes allowing the measurement of local material ferroelectric and dielectric of a modulated electric field to the conductive SPM tip results in the appearance of in-phase surface properties [50,75]. is a strain-based scanning probe microscopy [76,77]variations mode, where application displacement in aPFM pm–nm range. A lock-in detection technique in different is used for the of measurements a modulated electric field toand thephase conductive SPM tip results in the appearance in-phase of amplitude of piezoresponse [78]. The phase signal canofbe linked tosurface the displacement a pm–nm range. A lock-in detection variations is used for the spontaneousinpolarization orientation, while PFM technique amplitudeinisdifferent a function of local effective piezoelectric coefficient [79].and PFMphase allows only the visualization domain structure, also to thethe measurements of amplitude ofnot piezoresponse [78]. Theof phase signal can bebut linked quantitativepolarization determination of the polarization analysis of out-of-plane spontaneous orientation, while PFMorientation amplitudeby is asimultaneous function of local effective piezoelectric

Materials 2017, 10, 47

6 of 23

coefficient [79]. PFM allows not only the visualization of domain structure, but also the quantitative determination of47the polarization orientation by simultaneous analysis of out-of-plane and in-plane Materials 2017, 10, 6 of 23 piezoresponse signals (Vector PFM [80], 3D-PFM [80], and angle-resolved PFM [81]). Application of and enough in-plane DC piezoresponse (Vector PFM [80], [80], andfor angle-resolved PFMreversal) [81]). a high bias to the signals SPM tip (higher than the3D-PFM threshold field the polarization a high enough DC bias direction to the SPM tip (higher thanunder the threshold field forallows the canApplication reverse theof spontaneous polarization locally in the area the tip [82]. This polarization reversal) can reverse the spontaneous polarization direction locally in the area under thethe the measurement of local hysteresis loops of material and the study of domain wall motion in tip [82]. This allows the measurement of local hysteresis loops of material and the study of domain electric field of the probe. These measurements were quite rarely done on ceramics [77] because of the wall motion in the electric field of the probe. These measurements were quite rarely done on ceramics unknown orientation of a particular grain and interception with the grain boundaries. However, it was [77] because of the unknown orientation of a particular grain and interception with the grain possible to determine the intragrain domain wall velocity and other parameters, such as dimensionality boundaries. However, it was possible to determine the intragrain domain wall velocity and other of domain walls [83]. parameters, such as dimensionality of domain walls [83]. 3.4. Confocal Raman Microscopy 3.4. Confocal Raman Microscopy CRM is based on the study of the Raman spectra variations across the material surface, with spatial CRM is based on the study of the Raman spectra variations across the material surface, with resolution of about 300 nm provided by confocal microscopy [84]. The significant change of the Raman spatial resolution of about 300 nm provided by confocal microscopy [84]. The significant change of spectra (shifts and change of the intensity of the Raman bands) in the vicinity of the domain walls was the Raman spectra (shifts and change of the intensity of the Raman bands) in the vicinity of the shown in crystals, thin ceramics [9,51].and Moreover, method allowsthe notmethod only the domainsingle walls was shown in films, single and crystals, thin films, ceramicsthe [9,51]. Moreover, imaging 4) structure [51], but also the4)extraction of information about mechanical allowsof notdomain only thestructure imaging (Figure of domain (Figure [51], but also the extraction of information stresses [85,86] and defect concentration [51]. Polarized Raman scattering can yield knowledge about about mechanical stresses [85,86] and defect concentration [51]. Polarized Raman scattering can yield theknowledge orientationabout of spontaneous polarization in distinct grains of ceramics [9,87]. the orientation of spontaneous polarization in distinct grains of ceramics [9,87].

Figure 4. Characterization of KNN ceramics by confocal Raman spectroscopy: (a,b) Average Raman Figure 4. Characterization of KNN ceramics by confocal Raman spectroscopy: (a,b) Average Raman spectra of adjacent striped domains separated by ◦a 90° domain wall. The insets show magnified spectra of adjacent striped domains separated by a 90 domain wall. The insets show magnified Raman Raman spectra and Lorentzian fits of domain structure in the frequency range between 500 and 700−1 spectra and Lorentzian fits of domain structure in the frequency range between 500 and 700 cm . cm−1. These spectra are fitted to the sum of two Lorentzian peaks, ascribed to the Eg (υ2) and A1g (υ1) These spectra are fitted to the sum of two Lorentzian peaks, ascribed to the Eg (υ2 ) and A1g (υ1 ) Raman Raman modes, respectively; (c) Raman map of domain structure of the KNN exhibiting clear modes, respectively; (c) Raman map of domain structure of the KNN exhibiting clear differences differences between average spectra of adjacent striped domains separated by a 90° domain wall. The between average spectra of adjacent striped domains separated by a 90◦ domain wall. The Raman map Raman map was derived by summing up the total spectral pixel intensity from 100 to 1000 cm−1. was derived by summing up the total spectral pixel intensity from 100 to 1000 cm−1 . Adapted from [59], Adapted from [59], with permission from © 2012 Royal Society of Chemistry. with permission from © 2012 Royal Society of Chemistry.

4. Domain Structure in BiFeO3 4. Domain Structure in BiFeO3 4.1. Undoped BiFeO3 4.1. Undoped BiFeO3 The spontaneous polarization in BFO is oriented along the equivalent crystallographic direction Theand spontaneous polarization in BFO Based is oriented along thetypes equivalent crystallographic direction {111} has eight possible orientations. on this, three of domain walls are possible in {111} and has eight possible orientations. Based on this, three types of domain walls are possible BFO: 180°, 109°, and 71° (Figure 5). These angles are the rotation angles between neighboring in ◦ , and 71◦ (Figure 5). These angles are the rotation angles between neighboring domains. BFO: 180◦ , 109 domains. The permissible domain wall orientations are, therefore, {110} for 109°, {001} for 71°, and ◦ , {001} for The permissible domain orientations are,for therefore, {110} for 71◦ , and any plane any plane parallel to thewall polarization vector 180° domains. In a109 non-perfect crystal, however, it ◦ can be that thevector actualfor wall deviateInslightly from the crystallographically-predicted parallel to expected the polarization 180may domains. a non-perfect crystal, however, it can be expected planes. The smaller thedeviate wall area, the larger deviation. Four domain boundary types or that the actual wall may slightly fromthe thepossible crystallographically-predicted planes. The smaller possible, where “head-to-head” or domain “tail-to-tail” domaintypes wallsor are charged, whileare theconfigurations wall area, theare larger the possible deviation. Four boundary configurations “head-to-tail” “tail-to-head” walls are non-charged (canwalls be considered neutral). possible, where or “head-to-head” or “tail-to-tail” domain are charged, while “head-to-tail” or “tail-to-head” walls are non-charged (can be considered neutral).

Materials 2017, 10, 47 Materials 2017, 10, 47 Materials 2017, 10, 47

7 of 23 7 of 23 7 of 23

5. The schematic diagram of (001)-oriented BiFeO3 crystal structure and direction of the FigureFigure 5. The schematic diagram of (001)-oriented BiFeO3 crystal structure and direction of spontaneous polarization corresponding to the (a) 180°; (b) 109° and (c) 71° domains. Adapted from Figure 5. The schematic diagram of (001)-oriented BiFeO the the spontaneous polarization corresponding to the (a)3 crystal 180◦ ; structure (b) 109◦ and anddirection (c) 71◦ofdomains. [88], with permission from © 2006 Nature Publishing Group. spontaneous polarization corresponding to 2006 the (a) 180°; (b) 109° and (c) 71° domains. Adapted from Adapted from [88], with permission from © Nature Publishing Group. [88], with permission from © 2006 Nature Publishing Group.

Domain structure in BFO ceramics typically represents a mixture of domains with irregular

shapes separated by walls (so called, watermarks) and regular lamellar domains separated by Domain structure in 180° BFO ceramics typically represents a mixture of domains with irregular shapes Domain structure in BFO ceramics typically represents a mixture of domains with irregular ◦ ◦ non-180° walls (Figure 6) [89]. Being up to 10 μm wide, regular domains are an order of magnitude separated 180 walls (so called, watermarks) and regular lamellar domains separated by non-180 shapesby separated by 180° walls (so called, watermarks) and regular lamellar domains separated by larger than those reported for BFO thin films [89]. PFM contrast corresponding to the domain wallsnon-180° (Figure walls 6) [89]. Being up to 10 µm domains are an of magnitude (Figure 6) [89]. Being up towide, 10 μmregular wide, regular domains areorder an order of magnitudelarger structure varies among the grains, as[89]. expected non-oriented ceramics. The than thosesignificantly reported for BFO thin films PFM for contrast corresponding thenanoscale domain than larger those reported for BFO thin films [89]. PFM contrast corresponding to thetodomain structure domain structure at the boundary intersections observed in BFO by TEM can lead to high mechanical structure varies significantly among the grains, as expected for non-oriented ceramics. The nanoscale varies significantly among the grains, as expected for non-oriented ceramics. The nanoscale domain stress [90]. Recently, the existenceintersections of meta-stable polarization and vortex structures has been domain structure at the boundary BFOstates by high mechanical structure at the boundary intersections observed observed in BFO byinTEM canTEM leadcan to lead hightomechanical stress [90]. shown in BFO ceramics activation [91]. The domain sizes arehas typically stress [90]. Recently, theproduced existence by of mechanochemical meta-stable polarization states and vortex structures been Recently, existence of meta-stable polarization statesofand vortexceramics structures has the been shown in muchthe smaller than the grain sizes, domain contrast as-grown reflects intricate shown in BFO ceramics produced by and mechanochemical activation [91]. The domain sizes are typically BFO ceramics produced by mechanochemical activation [91]. The domain sizes are typically much interplay of mechanical stresses, charges, and defects accumulated at the grain much smaller than the grain sizes, uncompensated and domain contrast of as-grown ceramics reflects the intricate smaller than the grain sizes, and domain contrast of as-grown ceramics reflects the intricate interplay boundaries [92]. interplay of mechanical stresses, uncompensated charges, and defects accumulated at the grain of mechanical uncompensated charges, and defects accumulated at the grain boundaries [92]. boundariesstresses, [92].

Figure 6. (a) Out-of-plane and (b) in-plane PFM images of lamellar (red squares) and blotch (blue squares) domain structure in BiFeO3. Adapted from [89], with permission from © 2007 AIP Publishing Figure 6. (a) Out-of-plane and (b) in-plane PFM images of lamellar (red squares) and blotch (blue FigureLLC. 6. (a) Out-of-plane and (b) in-plane PFM images of lamellar (red squares) and blotch squares) domain structure in BiFeO3. Adapted from [89], with permission from © 2007 AIP Publishing (blue squares) domain structure in BiFeO3 . Adapted from [89], with permission from © 2007 AIP LLC.

Publishing 4.2. DopingLLC. by Rare Earth Ions

4.2. Doping by Rarethe Earth Ions Commonly, influence of rare earth element doping on the domain and phase structure

4.2. Doping by Earth Ions of the periodical domain structure to speckle-like domains with smaller results in:Rare (1) transformation Commonly, the influence of rare earth element doping on the domain and phase structure

sizes; in: (2) (1) appearance of non-polar/antipolar phase clusters (PFM signaland isdomains close noise with results transformation the periodical domain structure to speckle-like smaller Commonly, the influence ofofrare earth element doping on the domain phasetowith structure results increasing doping level. sizes; (2) appearance of non-polar/antipolar phase clusters (PFM signal is close to noise with in: (1) transformation of the periodical domain structure to speckle-like domains with smaller sizes; TEM doping and PFM studies of BixRe1-xFeO3 (Re = Sm, Gd, Dy) demonstrated concerted change of both increasing level. (2) appearance of non-polar/antipolar phase clusters (PFM signal is close to noise with increasing crystal structure and piezoresponse [21]. The change of the ferroelectric domain TEM and PFM studies of BixRe1-xcontrast FeO3 (Redistribution = Sm, Gd, Dy) demonstrated concerted change of both doping level. structure and phase composition as a function of the doping level was demonstrated for Sm-doped crystal structure and piezoresponse contrast distribution [21]. The change of the ferroelectric domain TEM and PFM [92]. studies Bix%Re Sm, Gd, is Dy) demonstrated concerted changeand of both 1-x FeO 3 (Re = structure BFO (8%–18%) At 8of mol Sm, domain mainly comprised of regular structure and phase composition as athe function of the doping level was demonstrated for lamella Sm-doped crystalwedges structure and piezoresponse contrast distribution [21]. The change of the ferroelectric domain of about 100–500 nm%inSm, size, to structure those observed in pristine BFO. the Sm content BFO (8%–18%) [92]. At 8 mol thesimilar domain is mainly comprised of As regular lamella andis increased to 12 and 14 mol %, the domains become progressively smaller (~50–200 nm) and more structure and phase composition as a function of the doping level was demonstrated for Sm-doped wedges of about 100–500 nm in size, similar to those observed in pristine BFO. As the Sm content is BFO irregularly shaped. The number ofdomains visible domains is also reduced as the composition (8%–18%) [92].toAt mol 14 % Sm, the domain structure is mainly comprised of (~50–200 regular lamella and wedges increased 128 and mol %, the become progressively smaller nm)approaches and more the MPB (ranging from 8 to 14 mol % Sm), to the extent that almost no domains are visible at 15.5 mol of about 100–500 nm in size, similar to those observed in pristine BFO. As the Sm content is increased to irregularly shaped. The number of visible domains is also reduced as the composition approaches % Sm (Figure 7). However, the regions with piezoelectric contrast keep a switchable behavior typical the14MPB from 8 tobecome 14 mol %progressively Sm), to the extent that almost no domains visible at 15.5 mol 12 and mol(ranging %, the domains smaller (~50–200 nm) andare more irregularly shaped. % Sm (Figure 7). However, regions with piezoelectric contrast keep a switchable behavior typical The number of visible domainsthe is also reduced as the composition approaches the MPB (ranging from 8 to 14 mol % Sm), to the extent that almost no domains are visible at 15.5 mol % Sm (Figure 7). However, the regions with piezoelectric contrast keep a switchable behavior typical of ferroelectric materials.

Materials 2017, 10, 47 Materials 2017, 10, 47

8 of 23 8 of 23

ferroelectric At the same significantly time, the phase concentration significantly dependsmethod on the [93]. At theofsame time, thematerials. phase concentration depends on the ceramic preparation method [93]. PFM contrast was attributed to consecutive R3c–Pbam–Pnma phase PFM ceramic contrastpreparation was attributed to consecutive R3c–Pbam–Pnma phase transformations. A similar effect of transformations. A similar effect of doping was revealed by TEM in Ndand Sm-doped BFO doping was revealed by TEM in Nd- and Sm-doped BFO compositions [94,95]. The compositions with compositions [94,95]. The compositions with x ≤ and 10% exhibited were rhombohedral with R3c and and x ≤ 10% were rhombohedral with R3c symmetry, superstructure andsymmetry, orientational exhibited superstructure and orientational and translational domains characteristic of an antiphasetranslational domains characteristic of an antiphase-tilted ferroelectric perovskite [95]. At the phase tilted ferroelectric perovskite [95]. At the phase boundary between the orthoferrite and rhombohedral boundary between the orthoferrite and rhombohedral cells in Nd- and Sm-doped systems, a new cells in Nd- and Sm-doped systems, a new structure is stabilized with a quadrupled unit cell similar structure stabilized with a quadrupled unit cell similar to the case of PbZrO3 [95]. to theiscase of PbZrO 3 [95].

Figure 7. Dependence of the domain structure on Sm2+ doping fraction. (a–d) PFM out-of-plane

Figure 7. Dependence of the domain structure on Sm2+ doping fraction. (a–d) PFM out-of-plane piezoresponse images obtained for 8 mol % Sm, 12 mol % Sm, 14 mol % Sm, and 15.5 mol % Sm, piezoresponse images obtained 8 mol % Sm, 12 % Sm, selected 14 mol % Sm, and 15.5 mol % Sm, respectively; (e,f) Bright field for (BF)-TEM images andmol associated area electron diffraction respectively; (e,f) Bright (BF)-TEM images and associated selected area electron diffraction (SAED) patterns fromfield 12 mol % Sm. (e) A region identified as R3c by corresponding SAED in [001] (SAED) pc patterns mol % Sm. (e) A region identified as R3c byregion corresponding in phase [001] pc zonefrom axis 12 (inset), where regular domains are seen. (f) A identified SAED as Pbam byzone axis (inset), where regular domains are seen. (f) Anano-sized region identified as observed. Pbam phase corresponding corresponding SAED (inset), where complicated features are Theby green-to-blue transition represents the change wt % ratioThe of R3c and Pbam color phases. PC SAEDcolor (inset), whereapproximately complicated nano-sized features areinobserved. green-to-blue transition denotes pseudocubic [92]. in wt % ratio of R3c and Pbam phases. PC denotes pseudocubic approximately representsnotation the change notation [92]. Local piezoresponse in Dy-substituted BFO is approximately three times weaker than in undoped ceramics, thus pointing to a smaller value of the spontaneous polarization and effective Local piezoresponse incan Dy-substituted is approximately times weaker than in piezocoefficient [96]. This be attributed to BFO the large difference in ionicthree radii of Bi3+ and Dy3+ ions, undoped ceramics, thus pointing to a smaller value of the spontaneous polarization and effective hampering the formation of homogeneous solid solutions. The coexistence of the regions 3+ and Dy3+ piezocoefficient [96]. This can attributed large difference ionic radii was of Biobserved demonstrating a distinct PFMbecontrast with to thethe areas showing a zero in piezoresponse the x = 0.15the compound. ions, for hampering formation of homogeneous solid solutions. The coexistence of the regions The piezoresponse in contrast Bi0.9Gd0.1FeO wasareas approximately weaker as was compared to for demonstrating a distinct PFM with3 the showing atwo zerotimes piezoresponse observed undoped BFO ceramics, due to a smaller value of the spontaneous polarization [97]. the x = 0.15 compound.

The piezoresponse in Bi0.9 Gd0.1 FeO3 was approximately two times weaker as compared to 4.3. Doping by Alkaline Earth and Heavy Metal Ions undoped BFO ceramics, due to a smaller value of the spontaneous polarization [97]. Domain structure of BFO ceramics with heterovalent substitution by Ca, Sr, Pb, and Ba ions was studiedbyinAlkaline Reference [17]. domain represented mostly non-oriented speckle-like 4.3. Doping Earth andThe Heavy Metalstructure Ions domains. At the same time, the average domain size strongly depended on the sintering

Domain structure of BFO ceramics with heterovalent substitution by Ca, Sr, Pb, and Ba ions was conditions [89]. studied inMn-doped ReferenceBFO [17].samples The domain structure mostly non-oriented speckle-like domains. exhibit a higherrepresented volume density of the domain walls than those of At theundoped same time, averagethat domain size depended onthe the sintering ones,the suggesting the Mn ionstrongly can effectively reduce domain size conditions in BFO [98]. [89]. Negative self-polarization (or polarization offset) was foundofinthe thedomain Pr and Sc co-substituted Mn-doped BFO samples exhibit a higher volume density walls than those of Bi 0.9 Pr 0.1 Fe 1−x Sc x O 3 (0.01 ≤ x ≤ 0.07), which had a maximum for about 3% Sc and 2% Pr [99]. undoped ones, suggesting that the Mn ion can effectively reduce the domain size in BFONegative [98]. Negative self-polarization (or polarization offset) was found in the Pr and Sc co-substituted Bi0.9 Pr0.1 Fe1−x Scx O3 (0.01 ≤ x ≤ 0.07), which had a maximum for about 3% Sc and 2% Pr [99]. Negative self-polarization in these samples can be a result of the built-in internal bias field Eint

Materials 2017, 10, 47

9 of 23

generated by excess electrons and charge defects like oxygen vacancies. These electrons trapped near the interface with the bottom electrode can form defect dipoles aligned during the formation of domain structure, and result in an internal bias field oriented towards the electrode. Domain structure in Co-, Ni-, Zn-, Nb-, and W-modified multiferroic BiFeO3 represented stripe domains with different orientations [100]. The domains studied by TEM were stable under the action of the electron beam, and their size did not vary much with the composition [100]. Obtained diffraction pattern in periodical domains revealed by high-resolution TEM showed apparent splitting of the electron diffraction spots perpendicular to the (110) planes, indicating the formation of the (110) domain walls [100]. 4.4. BiFeO3-x LaFeO3 -0.05La2/3 TiO3 Domain structure in BiFeO3-x LaFeO3 -0.05La2/3 TiO3 ceramics has been separated into domains with different length scales [101]: fine scale ferroelectric/ferroelastic twin domains (10–20 nm) and larger regions (100–200 nm), which define the domain structure associated with antiphase tilting. This fact suggests that the local direction of polarization and strain are inconsistent with the rhombohedral distortion of the macroscopic tilt system and symmetry. According to the proposed model, each tilt domain consisted of many tens of finer scale ferroelastic/ferroelectric twins (10–20 nm), which average polarization vector and spontaneous strain are consistent with the symmetry of the macroscopic tilt system, but with lower local symmetry. 4.5. Temperature Dependence of Local Piezoelectric Response As expected, the values of piezoelectric coefficients in RE-doped BFO ceramics decrease significantly with decreasing temperature [69,102]. Stabilization of the ceramics in a wide temperature range and improvement of the properties at close to room temperature is very important for various applications [103]. Pure BFO ceramics demonstrate a significant increase of the piezoelectric and dielectric properties at elevated temperatures [104]. The value of piezoelectric charge coefficient d33 measured at 1 Hz ranged from 33 pm/V at 24 ◦ C to 118 pm/V at 262 ◦ C. Dielectric permittivity demonstrates a similar trend and reaches a value above 3000 at elevated temperatures and low frequencies. This effect was attributed to Maxwell–Wagner relaxation from the grain boundaries and reversible motion of the conductive non-180◦ domain walls. The addition 9% of Dy stabilized piezoelectric properties up to 350 ◦ C, while BFO doped with Sm and Gd demonstrated strong decrease of d33 at temperatures above 200 ◦ C [69]. 4.6. Grain Size–Domain Size Relation Generally, the properties of ferroelectric and piezoelectric ceramics strongly depend on the grain size [105]. The effect of grain size on domain structure, dielectric, and piezoelectric properties in different materials has been extensively studied by Arlt and coworkers [60,106–108]. Castillo et al. [109] revealed that a decrease of the grain size in pure BFO ceramics led to an increase in the elementary cell volume and a decrease of the local piezoresponse. Grains with size below a few hundred nanometers are predominantly single domain [109]. At the same time, complex domain patterns including both 180◦ and ferroelastic domain walls are typical for micron-sized grains, similar to BaTiO3 -based ceramics [60]. 4.7. Local Switching by PFM The main disadvantage of BFO ceramics for piezoelectric applications is the combination of high electrical conductivity and high coercive field [11]. Due to the high leakage current and low breakdown field, it is difficult to reach saturation during the poling process. Several authors reported non-saturated polarization-electric-field (P-E) loops, which are often misinterpreted because of the high contribution of the conductive current [11]. Thus, the local measurement of piezoresponse hysteresis and local switching by the SPM tip look like an attractive method for studying polarization reversal properties in BFO ceramics [21,97].

Materials 2017, 10, 47

10 of 23

Materials 2017, 10, 47local switching by the SPM tip look like an attractive method for studying polarization 10 of 23 hysteresis and

reversal properties in BFO ceramics [21,97]. The original approach was proposed in BFO for the study of polarization reversal during The original approach was proposed in BFO for the study of polarization reversal during successive poling of a square area (10 × 10 μm2) under stepwise-increasing DC voltage (Figure 8) [89]. successive poling of a square area (10 × 10 µm22) under stepwise-increasing DC voltage (Figure 8) [89]. After each poling procedure, the 20 × 20 μm area was scanned again without DC bias. Thus, an After each poling procedure, the 20 × 20 µm2 area was scanned again without DC bias. Thus, an analog analog of the macroscopic polarization reversal was realized, and the main characteristic of this of the macroscopic polarization reversal was realized, and the main characteristic of this poling process poling process could be extracted from the average piezoelectric signal inside the poling square could be from the average piezoelectric signal inside the poling square versus applied DC bias. versusextracted applied DC bias.

Figure 8. Local switching by progressive step-by-step poling of the area. Out-of-plane PFM images of Figure 8. Local switching by progressive step-by-step poling of the area. Out-of-plane PFM images BiFeO3 ceramic after poling an inner area of 10 × 10 μm2 subsequently at different bias voltages: of BiFeO3 ceramic after poling an inner area of 10 × 10 µm2 subsequently at different bias voltages: (a) Virgin state Udc = 0; (b) +20; (c) +60; (d) 0; (e) −20; (f) −60; (g) 0; and (h) +20 V. The average PFM (a) Virgin state Udc = 0; (b) +20; (c) +60; (d) 0; (e) −20; (f) −60; (g) 0; and (h) +20 V. The average PFM signal of this area after these and more poling procedures in zero bias yields a hysteresis curve vs. Udc signal of this area after these and more poling procedures in zero bias yields a hysteresis curve vs. Udc (central inset). Adapted from [89], with permission from © 2007 AIP Publishing LLC. (central inset). Adapted from [89], with permission from © 2007 AIP Publishing LLC.

Another approach for the evaluation of the switching properties of doped BFO-based ceramics was demonstrated in Reference [93]. Theof bi-domain structure was created in a BFO-based single grainceramics using Another approach for the evaluation the switching properties of doped ±30 V DC voltagein bias voltage and by PFM.structure After that,was thecreated piezoresponse histograms of poled was demonstrated Reference [93].scanned The bi-domain in a single grain using ±30 V and bias co-substituted samples were created to evaluate average piezoresponse value. The DCBFO voltage voltage and scanned by PFM. After that, the the piezoresponse histograms of poled BFO dependence of the difference between negative and positive piezoelectric signals on the concentration and co-substituted samples were created to evaluate the average piezoresponse value. The dependence of Scdifference dopant demonstrated a maximum for 1% piezoelectric Sc and minimum for the 7% Sc co-substituted of the between negative and positive signals on the concentration of sample. Sc dopant

demonstrated a maximum for 1% Sc and minimum for the 7% Sc co-substituted sample. 5. Domain Structure and Local Piezoelectric Properties of KNN-Based Ceramics 5. Domain Structure and Local Piezoelectric Properties of KNN-Based Ceramics 5.1. Domain Structure before Poling 5.1. Domain Structure before Poling of KNN ceramics is represented by cubic grains with faceted grain The typical grain structure boundaries (Figure 9a)structure [110]. Such possess orthorhombic symmetry at room The typical grain ofceramics KNN ceramics is represented by cubic grainstemperature, with faceted which suggests the existence of 60°, 90°, 120°, and 180° domain walls in the Bmm2 structure [67,111]. grain boundaries (Figure 9a) [110]. Such ceramics possess orthorhombic symmetry at room The spontaneous polarization Ps of the orthorhombic phase is parallel to the {100} C direction, and 12 ◦ ◦ ◦ ◦ temperature, which suggests the existence of 60 , 90 , 120 , and 180 domain walls in the Bmm2 polarization orientations are permissible in the orthorhombic phase (Figure 10). Using pseudocubic structure [67,111]. The spontaneous polarization Ps of the orthorhombic phase is parallel to the {100}C coordinates, charged and uncharged 90° domain walls are confined to {100}C planes, while the direction, and 12 polarization orientations are permissible in the orthorhombic phase (Figure 10). charged 60° domain walls and uncharged 120° domain walls are limited to {110}C planes [67]. The Using pseudocubic coordinates, charged and uncharged 90◦ domain walls are confined to {100}C 180° domain walls are oriented parallel to the spontaneous polarization Ps. The indices of uncharged planes, while the charged 60◦ domain walls and uncharged 120◦ domain walls are limited to {110}C 60° domain walls and charged 120° domain walls depend on the piezoelectric or electrostrictive planes [67]. The 180◦ domain walls are oriented parallel to the spontaneous polarization Ps. The indices of uncharged 60◦ domain walls and charged 120◦ domain walls depend on the piezoelectric or electrostrictive coefficients [112,113]. Using pseudocubic coordinates, charged and uncharged 90◦

Materials 2017, 10, 47

11 of 23

Materials 2017, 10, 47

11 of 23

Materials 2017, 10, 47

11 of 23

coefficients [112,113]. Using pseudocubic coordinates, charged and uncharged 90° domain walls are domain walls are confined to {100}C planes, while the charged 60◦ domain walls and uncharged 120◦ confined to [112,113]. {100}C planes, while the charged 60° domain walls and coefficients Using pseudocubic coordinates, charged anduncharged uncharged120° 90° domain walls are domain walls are limited to {110}C planes [67]. The 180◦ domain walls are oriented parallel to the limited toto {110} planes while [67]. The 180° domain walls walls are oriented parallel to domain the spontaneous confined {100}CC planes, the charged 60° domain and uncharged 120° walls are spontaneous polarization Ps . The indices of uncharged 60◦ domain walls and charged 120◦ domain polarization Ps. The indices of uncharged 60° domain walls charged 120° domain depend limited to {110} C planes [67]. The 180° domain walls areand oriented parallel to the walls spontaneous walls onPsthe piezoelectric or electrostrictive coefficients ondepend the piezoelectric or electrostrictive coefficients [112,113]. polarization . The indices of uncharged 60° domain walls and [112,113]. charged 120° domain walls depend on the piezoelectric or electrostrictive coefficients [112,113].

Figure 9. Scanning electron microscopy of the (a) grain structure and (b) domain structure revealed

Figure 9. Scanning electron microscopy of the (a) grain structure and (b) domain structure revealed by by chemical etchingelectron in KNNmicroscopy ceramics. Adapted from [59],with permission from © 2012 Royalrevealed Society Figure 9. Scanning of the (a) grain structure and (b) domain structure chemical etching in KNN ceramics. Adapted from [59], with permission from © 2012 Royal Society of Chemistry. by chemical etching in KNN ceramics. Adapted from [59],with permission from © 2012 Royal Society of Chemistry. of Chemistry.

Figure 10. The schematic diagram of KNN crystal structure and possible directions of the spontaneous polarization. Figure 10. The schematic diagram of KNN crystal structure and possible directions of the spontaneous

Figure 10. The schematic diagram of KNN crystal structure and possible directions of the polarization. spontaneous The mostpolarization. usual domain structure types in KNN ceramics are so called “watermarks” and “herringbone” (Figure 9b). However, existence of “zigzag” and “square net” are also The most [60,67,70] usual domain structure types inthe KNN ceramics are so called “watermarks” and The most domain structure types inthe KNN ceramics are called “watermarks” reported in ausual few publications [59,114]. Such structures result freesoenergy minimization. Theand “herringbone” [60,67,70] (Figure 9b). However, existence offrom “zigzag” and “square net” are also “herringbone” [60,67,70] (Figure 9b). However, the existence of “zigzag” and “square net” herringbone contains parallel strips subdivided narrow V-shaped domains, which form reported in apattern few publications [59,114]. Such structuresby result from free energy minimization. Thea are alsoherringbone reported inpattern a fewby publications [59,114]. Such result from free energy minimization. 120° angle bisected the long domain wall. The structures wallsbyfor both short lamellar domains and long contains parallel strips subdivided narrow V-shaped domains, which form a domains arebisected oriented {211}parallel directions [70]. Domain was found to bedomains, dependent onlong the The120° herringbone pattern contains strips subdivided byboth narrow which form angle byalong the long domain wall. The walls structure for shortV-shaped lamellar domains and ◦ angle average grain size [115]. According to thewall. thermodynamic theory [116], thelamellar equilibrium domain size domains are oriented {211}domain directions [70].The Domain structure was found to be dependent on the a 120 bisected byalong the long walls for both short domains and long is proportional to the square root of domain wall energy, since domain size is determined by a balance average grain size [115]. According to the thermodynamic theory [116], the equilibrium domain size domains are oriented along {211} directions [70]. Domain structure was found to be dependent on between the energy of domain and energies ofsince depolarization and elastic fields by proportional to the[115]. square rootwall of domain wall energy, domain size is determined bycaused a balance the is average grain size According to the the thermodynamic theory [116], the equilibrium domain the spontaneous polarization and strain [117]. It wall wasofshown insince KNNdomain that elastic thesize fine exhibit the energy domain wall and energies depolarization and fields caused by by sizebetween is proportional to of the square root of the domain energy, isgrains determined predominant lamellar twinning, while in [117]. coarseItgrains, twinning with athat banded structure released the spontaneous polarization and strain was shown in KNN the fine grains exhibit a balance between the energy of domain wall and the energies of depolarization and elastic fields [118]. PFM study of KNN-based revealed a non-zero the unpoled state predominant lamellar twinning,ceramics while inalso coarse grains, twinningpiezoresponse with a bandedinstructure released caused by the spontaneous polarization and strain [117]. It was shown in KNN that the fine grains (self-polarization effect) [119]. ceramics also revealed a non-zero piezoresponse in the unpoled state [118]. PFM study of KNN-based exhibit predominant lamellar twinning, while in coarse grains, twinning with a banded structure (self-polarization effect) [119]. released [118]. Structure PFM study of KNN-based ceramics also revealed a non-zero piezoresponse in the 5.2. Domain after Poling unpoled state (self-polarization effect) [119]. 5.2. Domain Structure after Poling The domain structure in most grains becomes periodical after poling (Figure 3a–f). Some grains only a simple parallel of grains stripesbecomes extending over the whole grain,(Figure while other The domain structure in set most periodical after poling 3a–f).grains Someexhibit grains 5.2. exhibit Domain Structure after Poling two or more of parallel stripes [67,120]. over This the domain structure was other associated the exhibit only asets simple paralleldomain set of stripes extending whole grain, while grainswith exhibit Theordomain structure in most grains becomes periodical after poling (Figure 3a–f). grains formation of sets ferroelastic orthorhombic domains, termed 90°, or 120° domains [67]. TheSome average two more of parallel domain stripes [67,120]. This60°, domain structure was associated with the exhibit only a parallel set of stripes extending over90°, the grain, [67]. while other grains formation of simple ferroelastic orthorhombic domains, termed 60°, or whole 120° domains The average

exhibit two or more sets of parallel domain stripes [67,120]. This domain structure was associated

Materials 2017, 10, 47

12 of 23

with the formation of ferroelastic orthorhombic domains, termed 60◦ , 90◦ , or 120◦ domains [67]. The average domain Materials 2017, 10, 47 width is strongly dependent on the grain size, ranging from 100 nm to 1–3 12 ofµm 23 [67]. When several parallel domain stripes exist in one polycrystalline grain, the intersection angles formed ◦ or domain is strongly on the grain size,are ranging from 100 nm 1–3135 μm◦ .[67]. When to between the width adjacent sets ofdependent parallel domain stripes either around 45to According several parallel domain stripes exist in one polycrystalline grain, the intersection angles formed the intersection angles, these domain patterns were classified as domain configuration of type I and between the adjacent sets of parallel domain stripes are either around 45° or 135°. According to the type II, respectively (Figure 3g,h) [67]. Two models were proposed to explain such domain patterns. intersection angles,◦these◦ domain patterns were classified as domain configuration of type I and type One is composed of 90 , 60 , and 120◦ domain walls; the other one is composed of 180◦ , 90◦ , and 120◦ II, respectively (Figure 3g,h) [67]. Two models were proposed to explain such domain patterns. One domain walls. However, diverse crystallographic can one be observed at the polishing plane due to is composed of 90°, 60°, and 120° domain walls;planes the other is composed of 180°, 90°, and 120° the random orientation alignment of the polycrystalline grains, and values of the intersection domain walls. However, diverse crystallographic planes can be observed at the polishing plane dueangle between sets of parallel domain stripesofmay also show a large variation [67]. to the random orientation alignment the polycrystalline grains, and values of the intersection angle between of parallel domain stripesproperties may also show a large variation [67].Na0.50 )0.935 Li0.065 NbO3 and The poorsets stability of piezoelectric (aging) found in (K0.50 The poor stability piezoelectric (aging) effect found attributed in (K0.50Na0.50 0.065NbO3 and (K0.50 Na0.50 )0.92 Li0.08 NbO3 of ceramics is dueproperties to the depoling to)0.935 theLiformation of 180◦ (K0.50Na 0.50)0.92Li0.08NbO3 ceramics is due to the depoling effect to the formation 180°of the domains (Figure 11) [70]. The reorientation of 90◦ domains is attributed much more difficult thanofthat domains (Figure 11) [70]. The reorientation of 90° domains is much more difficult than that of the ◦ ◦ 180 domains in orthorhombic phase, being practically forbidden for 90 single domains due to 180° domains in orthorhombic phase, being practically forbidden for 90° single domains due to clamping [70]. The formation of 180◦ domains is commonly attributed to the reduction of electrostatic clamping [70]. The formation of 180° domains is commonly attributed to the reduction of electrostatic energy, and their presence significantly simplifies the depoling process. energy, and their presence significantly simplifies the depoling process.

Figure 11. In situ TEM observation during electrical poling: (a) 8 kV/cm; (b) 10 kV/cm; (c) 14 kV/cm;

Figure 11. In situ TEM observation during electrical poling: (a) 8 kV/cm; (b) 10 kV/cm; (c) 14 kV/cm; (d) 18 kV/cm. The direction of poling fields is indicated by the bright arrow in (b). Representative (d) 18 kV/cm. The direction of poling fields is indicated by the bright arrow in (b). Representative selected area diffraction pattern at each poling field are shown in the insets. Adapted from [70], with selected area diffraction pattern at each LLC. poling field are shown in the insets. Adapted from [70], permission from © 2013 AIP Publishing with permission from © 2013 AIP Publishing LLC. 5.3. Coexistence of Tetragonal and Orthorombic Phases

5.3. Coexistence of Tetragonal and Orthorombic Phases phases in KNN ceramics was found to depend The coexistence of orthorhombic and tetragonal on the doping conditions and phase content [59,70,121].phases Domain in poled samples containing The coexistence of orthorhombic and tetragonal inpatterns KNN ceramics was found to depend these two phases consist of lamellar domains that are slightly bent, and therefore, domain walls on the doping conditions and phase content [59,70,121]. Domain patterns in poled samples containing deviate from straight lines (Figure 12a). Such domains are usually discussed in terms of lattice these two phases consist of lamellar domains that are slightly bent, and therefore, domain walls deviate distortion between the orthorhombic and tetragonal crystal phases [121]. from straight lines (Figure 12a). Such domains are usually discussed in terms of lattice distortion between the orthorhombic and tetragonal crystal phases [121].

Materials 2017, 10, 47

13 of 23

Materials 2017, 10, 47

13 of 23

Materials 2017, 10, 47

13 of 23

Figure 12. (a) Bright-field image of NKNS-0.0375LT ceramic sample; (b) Scheme of domain Figure 12. (a) Bright-field image of NKNS-0.0375LT ceramic sample; (b) Scheme of domain morphology morphology evolution from single T phase to coexisted O and T phases. Adapted from [121], with evolution from single T phase to coexisted O and T phases. Adapted from [121], with permission from permission from © 2011 AIP Publishing LLC. Figure (a) Bright-field © 2011 AIP12. Publishing LLC. image of NKNS-0.0375LT ceramic sample; (b) Scheme of domain morphology evolution from single T phase to coexisted O and T phases. Adapted from [121], with The polymorphic phase boundary typically results in a decrease of polarization rather than in permission from © 2011 AIP Publishing LLC.

polymorphic phase boundary typically results inasafor decrease ofpiezoceramics polarization rather than anThe enhancement of dielectric and piezoelectric properties, MPBs in [59]. One of in anthe enhancement of dielectric and piezoelectric properties, as for MPBs in piezoceramics [59]. One of possible current understandings of the origin of the polymorphic phase boundary in KNN-type The polymorphic phase boundary typically results in a decrease of polarization rather than in thematerials possibleiscurrent understandings of distribution the origin of[59]. the A polymorphic boundary informed KNN-type basedofon specificand domain schematic ofphase structure[59]. at an enhancement dielectric piezoelectric properties, as for MPBs indomain piezoceramics One of materials is based on specific domain distribution [59]. A schematic of domain structure formed the PPB is represented in Figure 13 [59]. This shows two striped regions separated by a 90° domain the possible current understandings of the origin of the polymorphic phase boundary in KNN-type The is striped region on left-hand side contains 180° domains with the polarization direction at wall. the PPB represented inthe Figure 13distribution [59]. This [59]. shows two striped regions separated by aat90◦ materials is based on specific domain A schematic of domain structure formed perpendicular the surface, on right-hand side,180 an◦ alternation ofby tetragonal and domain wall. Thetostriped on side contains domains with the polarization the PPB is represented inregion Figurewhereas 13 the [59].left-hand Thisthe shows two striped regions separated a 90° domain orthorhombic phases istopresent, giving rise to contains 60° or 120° domains. The tetragonal direction perpendicular thethe surface, whereas on the right-hand side, anhigh-temperature alternation of tetragonal and wall. The striped region on left-hand side 180° domains with the polarization direction ◦ or 120the ◦ domains. domain structure serves as a template that constrains formation of the orthorhombic domains orthorhombic phases is present, giving rise to 60 The high-temperature tetragonal perpendicular to the surface, whereas on the right-hand side, an alternation of tetragonal and and eventually preserves tetragonal structure compensating the energy of the ceramic domain structure serves asthe a template constrains the formation of the orthorhombic domains orthorhombic phases is present, giving that rise to 60° by or 120° domains. Thesurface high-temperature tetragonal grains. In addition to the stability of the tetragonal symmetry at temperatures below that of theceramic TO–T and eventually preserves by compensating the energy of the domain structure servesthe as tetragonal a template structure that constrains the formation ofsurface the orthorhombic domains transition, the unusual relaxor behavior of this transition was also explained. The proposed domain and eventually preserves tetragonal structure by compensating surface energy of thethat ceramic grains. In addition to the the stability of the tetragonal symmetry atthe temperatures below of the structure clearly explains a reduction of the piezoelectric properties in lead-free ceramics due to grains. In addition to the stability of behavior the tetragonal symmetry at temperatures below that of the Tthe O–T TO–T transition, the unusual relaxor of this transition was also explained. The proposed limitation of the domain dynamics. On the other hand, excellent properties of the composition transition, the unusual of this was alsoproperties explained.inThe proposed domain domain structure clearly relaxor explainsbehavior a reduction of transition the piezoelectric lead-free ceramics due proposed by Saito et al. [31] were attributed to the following facts: (1) Tlead-free O–T transition is shifted close structure clearly explains a reduction of the piezoelectric properties in ceramics due to the to the limitation of the domain dynamics. On the other hand, excellent properties of the composition to room temperature, thusdynamics. the ceramic have highexcellent tetragonality ratio; (2) texturing of the limitation of theetdomain Onsamples the other hand, properties of the composition proposed by Saito al. [31] were attributed to the following facts: (1) T transition is shifted close to O–T of ceramics results in an excellent path for stress relief by the alignment polar orientations and proposed by Saitothus et al. [31] were attributed to the following facts: (1) TO–T(2) transition is of shifted close room temperature, the ceramic samples have high tetragonality ratio; texturing the ceramics probably by increasingthus the grain size thatsamples contributes tohigh better accommodation stresses. to room temperature, thestress ceramic ratio;of(2) texturing of theby results in an excellent path for relief by thehave alignmenttetragonality of polar orientations and probably ceramics results in an excellent path for stress relief by the alignment of polar orientations and increasing the grain size that contributes to better accommodation of stresses. probably by increasing the grain size that contributes to better accommodation of stresses.

Figure 13. Schematic representation of the three-dimensional domain structure. Adapted from [59], with permission from © 2012 Royal Society of Chemistry.

Figure 13.13. Schematic domainstructure. structure.Adapted Adapted from [59], Figure Schematicrepresentation representationof ofthe the three-dimensional three-dimensional domain from [59], with permission from ©©2012 with permission from 2012Royal RoyalSociety Societyof ofChemistry. Chemistry.

5.4. Temperaute Dependences of Local Piezoelectric Response and Grain Size–Domain Size Relation

Domain structure transformation as a function of temperature has been studied by PFM [122]. 5.4. Temperaute Dependences LocalPiezoelectric Piezoelectric Response Relation 5.4.The Temperaute ofofLocal and Grain Grain Size–Domain Size–Domain Size Relation transitionDependences into the tetragonal phase occursResponse at approximately 65 °C, and isSize followed up by the reconstruction of the domain structureas in the temperature range between 50 and 130 by °C.PFM At higher Domain structure transformation function of [122]. Domain structure transformation as a function of temperature temperaturehas hasbeen beenstudied studied by PFM [122]. temperatures, domains start to disappear first in larger grains, being more stable in smaller grains. ◦ C,and The transition intothe thetetragonal tetragonalphase phase occurs occurs at at approximately approximately 65 upup byby thethe The transition into 65°C, andisisfollowed followed Finally, at 210 °C, onlydomain some residual nanodomains are observed at the grain boundaries 14). ◦(Figure reconstruction the structure in the the temperature temperature range AtAt higher reconstruction ofofthe domain structure in rangebetween between5050and and130 130°C. C. higher temperatures, domains start to disappear first in larger grains, being more stable in smaller grains. Finally, at 210 °C, only some residual nanodomains are observed at the grain boundaries (Figure 14).

Materials 2017, 10, 47

14 of 23

temperatures, domains start to disappear first in larger grains, being more stable in smaller grains. ◦ C, only some residual nanodomains are observed at the grain boundaries (Figure 14). Finally, at 210 Materials 2017, 10, 47 14 of 23

Figure 14. Vertical and lateral piezoresponse force microscopy (VPFM and LPFM, respectively) Figure 14. Vertical and lateral piezoresponse force microscopy (VPFM and LPFM, respectively) images of the polished surface of 0.95(Na0.49 K0.49Li0.02)(Nb0.8Ta0.2)O3-0.05CaZrO3 ceramics taken upon images of the polished surface of 0.95(Na0.49 K0.49 Li0.02 )(Nb0.8 Ta0.2 )O3 -0.05CaZrO3 ceramics taken heating at 30 °C, 145◦°C, and◦210 °C. Adapted from [122], with permission from © 2014 AIP Publishing upon heating at 30 C, 145 C, and 210 ◦ C. Adapted from [122], with permission from © 2014 AIP LLC. Publishing LLC.

Domain structure transformation with temperature was also studied in Domain structure transformation with temperature was also studied in 0.96(K 0.4Na0.6)(Nb0.96Sb0.04)O3-0.04Bi0.5K0.5Zr0.85Sn0.15O3 ceramics [123]. The process leads to the 0.96(K )O3 -0.04Bibetween O3 50ceramics [123]. The process leads 0.4 Na0.6 )(Nb 0.96 Sb0.04structures 0.5 K0.5 Zr0.85 reorganization of domain 28 Sn °C0.15 and °C, resulting in different temperature◦ C and 50 ◦ C, resulting in different to the reorganization of domain structures between 28 dependent behaviors of regular and irregular domains [123]. At high temperature (above 100 °C), the temperature-dependent behaviors of regular and irregular domains [123]. At high temperature amplitude of piezoresponse significantly decreased, accompanied with obscured ferroelectric ◦ (above 100 C), the amplitude of piezoresponse significantly decreased, accompanied with obscured domain boundaries. At lower temperatures (50–100 °C), indiscriminate behaviors of the multi-scale ◦ C), indiscriminate behaviors of the ferroelectric domain boundaries. At lower temperatures (50–100 domains may be related to the single-phase state (T state). The nano-scale domains possess improved multi-scale domains be related to thewall single-phase statecould (T state). Therespond nano-scale domains stimuli, possess flexibility due to themay reduced domain energy, and easily to external improved flexibility due to the reduced domain wall energy, and could easily respond to external contributing to the piezoelectric performance. stimuli, contributing to the piezoelectric performance. Li0.02(K0.45Na0.55)0.98NbO3 (LKNN) ceramics demonstrated a clear dependence of domain structure Li NbO3 (LKNN) ceramics a clear dependence of size domain 0.02 (K0.45 0.55 )0.98 on the grain sizeNa [124]. A decrease of the period of 90°demonstrated lamellar domains on decreasing grain was ◦ lamellar domains on decreasing structure on the grain size [124]. A decrease of the period of 90 observed down to 3 μm. Below 3 μm, an opposite dependence of domain period on grain size was grain observed down 3 µm. Below 3 µm, an opposite dependence of domain period on found.size Thewas grain size effect wastostudied in 0.95(K 0.5Na0.5)NbO3-0.05BaTiO3 (KNN-BT) ceramics with grain sizesecondary was found.milling The grain size effect was studiedgrain in 0.95(K Na0.5 )NbOof 0.5Reduction 3 -0.05BaTiO 3 (KNN-BT) different times, leading to different sizes. the grain size from 14 ceramics with different secondary milling times, leading to different grain sizes. Reduction of the grain to 1 μm increased the maximum dielectric permittivity from 2600 to 5000. This fact was attributed to size from 14 to 1 µm increased the maximum dielectric permittivity from 2600 to 5000. This fact was high internal non-uniformly distributed stress in grains resulting in an easier domain wall motion attributed to high internal non-uniformly distributed stress in grains resulting in an easier domain [115]. wall motion [115]. 5.5. Local Switching by PFM 5.5. Local Switching by PFM Dependence of the local domain switching on the domain structure period was observed in Dependence of the local domain switching on the domain structure period was observed in Reference [82]. It was demonstrated that the decrease of the average domain size was accompanied Reference [82]. It was demonstrated that the decrease of the average domain size was accompanied by by the decrease of coercive voltage [115]. The temperature increase significantly modified hysteresis the decrease of coercive voltage [115]. The temperature increase significantly modified hysteresis loop loop shape [122]. The remnant and maximal piezoresponse reached their peak values at shape [122]. The remnant and maximal piezoresponse reached their peak values at approximately approximately 55 °C, which was attributed to the orthorhombic–tetragonal transition that occurred 55 ◦ C, which was attributed to the orthorhombic–tetragonal transition that occurred in the studied in the studied composition at approximately 65 °C. composition at approximately 65 ◦ C. Rapid decrease of the piezoresponse amplitude and remarkable broadening of the local loops Rapid decrease of the piezoresponse amplitude and remarkable broadening of the local loops were obtained in the tetragonal phase. It was assumed that nucleation of new domains contributed were obtained in the tetragonal phase. It was assumed that nucleation of new domains contributed to to polarization reversal, meaning higher remnant polarization and maximal switchable polarization reversal, meaning higher remnant polarization and maximal switchable piezoresponse. piezoresponse. The transition into the tetragonal state was accompanied by a change in local The transition into the tetragonal state was accompanied by a change in local polarization switching polarization switching kinetics: the domain nucleation under the PFM tip occurred at a larger bias kinetics: the domain nucleation under the PFM tip occurred at a larger bias voltage, while propagation voltage, while propagation of the new domain was faster than in the orthorhombic phase. The less saturated local piezoelectric hysteresis loops and increased difference between the maximal and remnant piezoresponse indicated that the switching was incomplete and unstable. Thus, higher temperature promoted the effect of backswitching of the reversed domains.

Materials 2017, 10, 47

15 of 23

of the new domain was faster than in the orthorhombic phase. The less saturated local piezoelectric hysteresis loops and increased difference between the maximal and remnant piezoresponse indicated that the switching was incomplete and unstable. Thus, higher temperature promoted the effect of backswitching of the reversed domains. 6. Domain Structure Input to the Dielectric Permittivity and Piezoelectricity Relation between domain structure and dielectric permittivity in wide frequency range is well known in piezoelectric ceramics [107,125–128]. A number of models were suggested to describe the mechanism of the effect: 1. 2.

3. 4.

Domain wall vibrations under the action of applied electric field (resonance domain wall oscillations) [127–130] described by empirical Rayleigh law [126–128]. Emission of elastic shear waves from ferroelastic domain walls [129,131–135]. Ferroelastic domain walls are displaced in an applied electric field, which causes a shift of matter on both sides of the wall in opposite directions parallel to the domain wall. Thus, the domain walls behave like shear wave transducers; the shear waves emitted into the adjacent grains at high frequencies cause considerable dielectric losses, resulting in dielectric relaxation. Piezoelectric sound generation by laminar stacks of 180◦ domains (piezoelectric resonance of domains) [133–135]. Maxwell–Wagner effect of the conductive domain walls [136].

However, some experimental data exists related to direct measurements of the influence of domain structure on the resulting properties of ferroelectric ceramics. Strong dielectric dispersion in the MHz to GHz range was found in pristine KNN ceramics, and was attributed to the dynamics of ferroelastic–ferroelectric domain walls [137]. The oscillation of sets of laminar ferroelastic domain walls was proposed to emit transverse acoustic waves [128]. The nonlinear dielectric response was found in KNN-0.05LT piezoelectric ceramics that underwent an abrupt fall at the polymorphic phase transition due to the change of domain configuration [137]. The dielectric Rayleigh coefficient in the tetragonal phase zone was smaller than that in the orthorhombic one. Additionally, during the polymorphic phase transition (PPT), both the intrinsic and the extrinsic contributions to the dielectric constant significantly fluctuated with temperature variation. The nonlinear dielectric response at low frequencies and irreversible extrinsic contribution to dielectric permittivity found in KNN-LTS ceramics were strongly dependent on the domain configuration [130], determined by the crystallographic structure. The spontaneous distortion in tetragonal phase is larger than that in the orthorhombic one, which constrains the domain wall motion and results in smaller irreversible extrinsic contribution for the tetragonal compositions. The behavior in the PPT region is controlled by two effects: on one hand, the high domain wall concentration decreases the wall mobility due to self-clamping, while on the other hand, the additional spontaneous polarization emerging from the phase coexistence enhances the domain wall mobility. Due to the combination of both effects, the irreversible extrinsic contribution in the PPT is similar to that in the tetragonal region, and lower than that in the orthorhombic one. Dielectric permittivity demonstrated non-monotonic dependence on the domain wall concentration: first increasing with density, and then falling down [130]. The intrinsic dielectric permittivity increased, while irreversible extrinsic contribution decreased after poling [130]. In the Sb5+ -modified KNN-based ceramics, domain-wall vibrations were found to change from the resonance to relaxation mode due to a substantial increase of damping constant [129]. Doping by Mn influences the mobility of ferroelastic domains [138]. Sufficiently strong domain motion was proposed due to clamping along the vertical direction to grain bars. Internal friction occurred near the grain boundaries and ferroelectric domain walls during polarization reversal in the external electric field. Therefore, the friction force significantly affects the dielectric response [139].

Materials 2017, 10, 47

16 of 23

The role of domain structure and nanodomains appeared close to the PPT region in KNN ceramics remains poorly studied. Significant variation in the electrical and mechanical properties was observed in KNN before and after poling [140]. The highest coercive field and the lowest field-induced strain were obtained in the pristine samples of the PPB composition, while the opposite trend was detected after poling [140]. Unusual behavior of KNN ceramics close to PPT is usually attributed to the appearance of nanodomains due to two coexisting phases. High domain wall concentration leads to a larger domain wall energy barrier in the PPT region, making the domain wall motion difficult. Materials 2017, 10, 47 16 of 23 Therefore, the coercive field of the pristine samples reaches a maximum in the vicinity of the PPT region During nanodomains transformtransform irreversibly into micron-sized domains. the PPT[140]. region [140].poling, Duringthe poling, the nanodomains irreversibly into micron-sized The domains in the PPT region move more easily due to the existence of additional polarization states. domains. The domains in the PPT region move more easily due to the existence of additional Therefore, thestates. coercive field of the ceramics has a minimum value has at the PPT composition [140]. polarization Therefore, the poled coercive field of the poled ceramics a minimum value at the At the same time, it was proposed that nanodomains could be responsible for enhanced electric PPT composition [140]. At the same time, it was proposed that nanodomains could be responsible for field-induced strain [141]. enhanced electric field-induced strain [141]. Hayati et al. found that the addition of of nanosized nanosized ZnO ZnO additive additive increased increased the the grain grain size size of of the Hayati et al. found that the addition the KNN ceramics, with a corresponding increase of the domain size. This increase was accompanied KNN ceramics, with a corresponding increase of the domain size. This increase was accompanied 2+ ions entered in B-positions, with with aa reduction reduction of of the the coercive coercive field field [142]. [142]. For For higher higher doping doping levels, levels, Zn Zn2+ ions entered in B-positions, resulting in the formation of oxygen vacancies that pinned domain walls and concomitant increase of resulting in the formation of oxygen vacancies that pinned domain walls and concomitant increase the coercive field [142]. of the coercive field [142]. An exhaustive discussion about about indirect indirect measurements measurements of of domain domain wall wall input input to to dielectric An exhaustive discussion dielectric and and piezoelectric properties of BFO was reviewed by Rojac and co-authors [11]. The poling of rhombohedral piezoelectric properties of BFO was reviewed by Rojac and co-authors [11]. The poling of BiFeO3 ceramics through 71◦ and 109◦ domain wall reorientation was shown and quantified using rhombohedral BiFeO 3 ceramics through 71° and 109° domain wall reorientation was shown and analysis ofusing X-rayanalysis diffraction band diffraction intensities band [11]. intensities The non-linear effects related to the motion quantified of X-ray [11]. The non-linear effects related of to conductive domain walls in BFO bulk ceramics (non-linear Maxwell–Wagner effect) were the motion of conductive domain walls in BFO bulk ceramics (non-linear Maxwell–Wagnerrecently effect) demonstrated [105]. The domain grainwall boundary conductivity been measured by were recently demonstrated [105].wall The and domain and grain boundary have conductivity have been conductive atomic force microscopy (Figure 15), and mobility of domain walls under sub-switching measured by conductive atomic force microscopy (Figure 15), and mobility of domain walls under electric fields has beenfields determined BFO ceramics [7,136]. The high resolution imaging of sub-switching electric has beenindetermined in BFO ceramics [7,136]. The highSTEM resolution STEM the individual wallsdomain was done together with inspection of their conductivity depends on the imaging of thedomain individual walls was done together with inspection of their conductivity sintering conditions [7]. The conductivity of domain walls in bulk BFO ceramics was suggested be depends on the sintering conditions [7]. The conductivity of domain walls in bulk BFO ceramicsto was 4+ 3+ p-type, related the polaron mechanism between Fe and Fe sites suggested to betop-type, relatedhopping to the polaron hopping mechanism between Fe4+[7]. and Fe3+ sites [7].

Figure 15. (a) Out-of-plane PFM amplitude and (b) conductive atomic force microscopy (c-AFM) Figure 15. (a) Out-of-plane PFM amplitude and (b) conductive atomic force microscopy (c-AFM) images of selected regions with domains in BiFeO3 annealed in O2 atmosphere. Adapted from [7], images of selected regions with domains in BiFeO3 annealed in O2 atmosphere. Adapted from [7], with permission from © 2016 Nature Publishing Group. with permission from © 2016 Nature Publishing Group.

7. Conclusions 7. Conclusions The domain dynamics plays an important role in the physical properties of lead-free The domain dynamics plays an important role in the physical properties of lead-free piezoelectric piezoelectric ceramics properties, and, as such, should be rigorously studied. The rapid development ceramics properties, and, as such, should be rigorously studied. The rapid development of modern of modern microscopic techniques gives the opportunity not only to visualize domain structure, but microscopic techniques gives the opportunity not only to visualize domain structure, but also to also to quantitatively determine the orientation and values of spontaneous polarization and quantitatively determine the orientation and values of spontaneous polarization and piezoelectric piezoelectric coefficient locally. The morphotropic and polymorphic phase boundaries existing both coefficient locally. The morphotropic and polymorphic phase boundaries existing both in BFO and in BFO and KNN are principally coined to the transformation of the domain assemblages and phase structure. Ferroelectric nanodomains and inhomogeneously distributed mixed phases appear due to doping and peculiarities of the sintering process, and are extremely important for future applications. This review is one of the first attempts to describe the domain structure and local piezoelectric properties and their impact on the macroscopic performance of lead-free piezoelectric ceramics based

Materials 2017, 10, 47

17 of 23

KNN are principally coined to the transformation of the domain assemblages and phase structure. Ferroelectric nanodomains and inhomogeneously distributed mixed phases appear due to doping and peculiarities of the sintering process, and are extremely important for future applications. This review is one of the first attempts to describe the domain structure and local piezoelectric properties and their impact on the macroscopic performance of lead-free piezoelectric ceramics based on BFO and KNN. Acknowledgments: The equipment of the Ural Center for Shared Use “Modern Nanotechnology” UrFU has been used. The research was made possible by the Ministry of Education and Science of Russian Federation (UID RFMEFI58715X0022). Conflicts of Interest: The authors declare no conflict of interest.

References 1. 2. 3. 4. 5. 6.

7.

8. 9. 10. 11. 12. 13. 14. 15. 16. 17.

18. 19.

Hong, C.-H.; Kim, H.-P.; Choi, B.-Y.; Han, H.-S.; Son, J.S.; Ahn, C.W.; Jo, W. Lead-free piezoceramics—Where to move on? J. Materiomics 2016, 2, 1–24. [CrossRef] Irschik, H.; Krommer, M.; Vetyukov, Y. On the use of piezoelectric sensors in structural mechanics: Some novel strategies. Sensors 2010, 10, 5626–5641. [CrossRef] [PubMed] Safari, A. Novel piezoelectric ceramics and composites for sensor and actuator applications. Mater. Res. Innov. 1999, 52, 263–269. [CrossRef] Aksel, E.; Jones, J.L. Advances in lead-free piezoelectric materials for sensors and actuators. Sensors 2010, 10, 1935–1954. [CrossRef] [PubMed] Spanner, K.; Koc, B. Piezoelectric motors, an overview. Actuators 2016, 5, 6. [CrossRef] Directive 2002/95/EC of the European Parliament and of the Council of 27 January 2003 on the Restriction of the Use of Certain Hazardous Substances in Electrical and Electronic Equipment; Directive 2002/95/EC; European Union: Brussels, Belgium, 2003. Rojac, T.; Bencan, A.; Drazic, G.; Sakamoto, N.; Ursic, H.; Jancar, B.; Tavcar, G.; Makarovic, M.; Walker, J.; Malic, B.; et al. Domain-wall conduction in ferroelectric BiFeO3 controlled by accumulation of charged defects. Nat. Mater. 2016. [CrossRef] [PubMed] Potnis, P.R.; Tsou, N.T.; Huber, J.E. A review of domain modelling and domain imaging techniques in ferroelectric crystals. Materials 2011, 4, 417–447. [CrossRef] Shur, V.Y.; Zelenovskiy, P.S. Micro- and nanodomain imaging in uniaxial ferroelectrics: Joint application of optical, confocal Raman, and piezoelectric force microscopy. J. Appl. Phys. 2014, 116, 066802. [CrossRef] Liu, H.; Wang, X. Large electric polarization in BiFeO3 film prepared via a simple sol-gel process. J. Sol-Gel Sci. Technol. 2008, 47, 154–157. [CrossRef] Rojac, T.; Bencan, A.; Malic, B.; Tutuncu, G.; Jones, J.L.; Daniels, J.E.; Damjanovic, D. BiFeO3 ceramics: Processing, electrical, and electromechanical properties. J. Am. Ceram. Soc. 2014, 97, 1993–2011. [CrossRef] Maître, A.; François, M.; Gachon, J.C. Experimental study of the Bi2 O3 -Fe2 O3 pseudo-binary system. J. Phase Equilibria Diffus. 2004, 25, 59–67. [CrossRef] Achenbach, G.D.; James, W.J. Preparation of single-phase polycrystalline BiFeO3 . J. Am. Ceram. Soc. 1967, 50, 437. [CrossRef] Selbach, S.M.; Einarsrud, M.A.; Grande, T. On the thermodynamic stability of BiFeO3 . Chem. Mater. 2009, 21, 169–173. [CrossRef] Selbach, S.M.; Tybell, T.; Einarsrud, M.; Grande, T. Size-dependent properties of multiferroic BiFeO3 nanoparticles. Chem. Mater. 2007, 19, 6478–6484. [CrossRef] Valant, M.; Axelsson, A.; Alford, N.; Bifeo, P. Peculiarities of a solid-state synthesis of multiferroic polycrystalline BiFeO3 . Chem. Mater. 2007, 19, 5431–5436. [CrossRef] Khomchenko, V.A.; Kiselev, D.A.; Vieira, J.M.; Jian, L.; Kholkin, A.L.; Lopes, A.M.L.; Pogorelov, Y.G.; Araujo, J.P.; Maglione, M. Effect of diamagnetic Ca, Sr, Pb, and Ba substitution on the crystal structure and multiferroic properties of the BiFeO3 perovskite. J. Appl. Phys. 2008, 103, 024105. [CrossRef] Nalwa, K.S.; Garg, A.; Upadhyaya, A. Effect of samarium doping on the properties of solid-state synthesized multiferroic bismuth ferrite. Mater. Lett. 2008, 62, 878–881. [CrossRef] Leist, T.; Granzow, T.; Jo, W.; Rödel, J. Effect of tetragonal distortion on ferroelectric domain switching: A case study on La-doped BiFeO3 -PbTiO3 ceramics. J. Appl. Phys. 2010, 108, 014103. [CrossRef]

Materials 2017, 10, 47

20. 21. 22.

23. 24. 25. 26. 27. 28. 29.

30.

31. 32. 33. 34.

35. 36. 37. 38. 39. 40. 41.

42. 43.

18 of 23

Ozaki, T.; Matsuo, H.; Noguchi, Y.; Miyayama, M.; Mori, S. Microstructures related to ferroelectric properties in (Bi0.5 K0.5 )TiO3 -BiFeO3 . Jpn. J. Appl. Phys. 2010, 49, 09MC05. [CrossRef] Khomchenko, V.A.; Paixão, J.A.; Kiselev, D.A.; Kholkin, A.L. Intermediate structural phases in rare-earth substituted BiFeO3 . Mater. Res. Bull. 2010, 45, 416–419. [CrossRef] Fujino, S.; Murakami, M.; Anbusathaiah, V.; Lim, S.H.; Nagarajan, V.; Fennie, C.J.; Wuttig, M.; Salamanca-Riba, L.; Takeuchi, I. Combinatorial discovery of a lead-free morphotropic phase boundary in a thin-film piezoelectric perovskite. Appl. Phys. Lett. 2008, 92, 202904. [CrossRef] Palai, R.; Katiyar, R.S.; Schmid, H.; Tissot, P.; Clark, S.J.; Robertson, J.; Redfern, S.A.T.; Catalan, G.; Scott, J.F. β phase and γ-β metal-insulator transition in multiferroic BiFeO3 . Phys. Rev. B 2008, 77, 014110. [CrossRef] Jaeger, R.E.; Egerton, L. Hot pressing of potassium-sodium niobates. J. Am. Ceram. Soc. 1962, 45, 209–213. [CrossRef] Ahtee, M.; Hewat, A.W. Structural phase transitions in sodium-potassium niobate solid solutions by neutron powder diffraction. Acta Crystallogr. A 1978, 34, 309–317. [CrossRef] Ahtee, M.; Glazer, A.M. Lattice parameters and tilted octahedra in sodium-potassium niobate solid solutions. Acta Crystallogr. A 1976, 32, 434–446. [CrossRef] Jaffe, B.; Cook, W.R.; Jaffe, H.L. Piezoelectric Ceramics; Academic Press: Cambridge, MA, USA, 1971. Birol, H.; Damjanovic, D.; Setter, N. Preparation and characterization of (K0.5 Na0.5 )NbO3 ceramics. J. Eur. Ceram. Soc. 2006, 26, 861–866. [CrossRef] Zhang, B.P.; Li, J.F.; Wang, K.; Zhang, H. Compositional dependence of piezoelectric properties in Nax K1-x NbO3 lead-free ceramics prepared by spark plasma sintering. J. Am. Ceram. Soc. 2006, 89, 1605–1609. [CrossRef] Li, J.F.; Wang, K.; Zhang, B.P.; Zhang, L.M. Ferroelectric and piezoelectric properties of fine-grained Na0.5 K0.5 NbO3 lead-free piezoelectric ceramics prepared by spark plasma sintering. J. Am. Ceram. Soc. 2006, 89, 706–709. [CrossRef] Saito, Y.; Takao, H.; Tani, T.; Nonoyama, T.; Takatori, K.; Homma, T.; Nagaya, T.; Nakamura, M. Lead-free piezoceramics. Nature 2004, 432, 84–87. [CrossRef] [PubMed] Baker, D.W.; Thomas, P.A.; Zhang, N.; Glazer, A.M. A comprehensive study of the phase diagram of Kx Na1-x NbO3 . Appl. Phys. Lett. 2009, 95, 91903. [CrossRef] Maliˇc, B.; Koruza, J.; Hrešˇcak, J.; Bernard, J.; Wang, K.; Fisher, J.; Benˇcan, A. Sintering of lead-free piezoelectric sodium potassium niobate ceramics. Materials 2015, 8, 8117–8146. [CrossRef] Fuentes, J.; Portelles, J.; Pérez, A.; Durruthy-Rodríguez, M.D.; Ostos, C.; Raymond, O.; Heiras, J.; Cruz, M.P.; Siqueiros, J.M. Structural and dielectric properties of La- and Ti-modified K0.5 Na0.5 NbO3 ceramics. Appl. Phys. A 2012, 107, 733–738. [CrossRef] Zhou, J.J.; Li, J.F.; Wan, K.; Zhang, X.W. Phase structure and electrical properties of (Li,Ta)-doped (K, Na)NbO3 lead-free piezoceramics in the vicinity of Na/K = 50/50. J. Mater. Sci. 2011, 46, 5111–5116. [CrossRef] Malic, B.; Bernard, J.; Holc, J.; Kosec, M. Strontium doped K0.5 Na0.5 NbO3 based piezoceramics. Ferroelectrics 2005, 314, 149–156. [CrossRef] Malic, B.; Bernard, J.; Holc, J.; Jenko, D.; Kosec, M. Alkaline-earth doping in (K, Na)NbO3 based piezoceramics. J. Eur. Ceram. Soc. 2005, 25, 2707–2711. [CrossRef] Guo, Y.; Kakimoto, K.; Ohsato, H. (Na0.5 K0.5 )NbO3 -LiTaO3 lead-free piezoelectric ceramics. Mater. Lett. 2005, 59, 241–244. [CrossRef] Fang, B.; Du, Q.; Jiang, N.; Wu, J. Effects of the co-addition of LiSbO3 -LiTaO3 on the densification of (Na1/2 K1/2 )NbO3 lead free ceramics by atmosphere sintering. J. Alloys Compd. 2011, 509, 2420–2424. Wu, J.; Xiao, D.; Zhu, J. Potassium-sodium niobate lead-free piezoelectric materials: Past, present, and future of phase boundaries. Chem. Rev. 2015, 115, 2559–2595. [CrossRef] [PubMed] Li, J.F.; Wang, K.; Zhu, F.Y.; Cheng, L.Q.; Yao, F.Z. (K, Na) NbO3-based lead-free piezoceramics: Fundamental aspects, processing technologies, and remaining challenges. J. Am. Ceram. Soc. 2013, 96, 3677–3696. [CrossRef] Reichmann, A.; Mitsche, S.; Zankel, A.; Poelt, P.; Reichmann, K. In situ mechanical compression of polycrystalline BaTiO3 in the ESEM. J. Eur. Ceram. Soc. 2014, 34, 2211–2215. [CrossRef] Bradt, R.C.; Ansell, G.S. Transmission electron microscopy of ferroelectric domain boundaries in barium titanate. IEEE Trans. Electron Devices 1969, 16, 594–597. [CrossRef]

Materials 2017, 10, 47

44.

45.

46. 47. 48. 49.

50. 51.

52. 53. 54. 55. 56. 57. 58. 59. 60. 61. 62.

63.

64.

65. 66.

67.

19 of 23

Roshchupkin, D.V.; Irzhak, D.V.; Antipov, V.V. Study of LiNbO3 and LiTaO3 ferroelectric domain structures using high-resolution x-ray diffraction under application of external electric field. J. Appl. Phys. 2009, 105, 024112. [CrossRef] Antipov, V.V.; Blistanov, A.A.; Roshchupkina, E.D.; Tucoulou, R.; Ortega, L.; Roshchupkin, D.V. High-resolution X-ray topography and diffraction study of bulk regular domain structures in LiNbO3 crystals. Appl. Phys. Lett. 2004, 85, 5325. [CrossRef] Takashige, M.; Hamazaki, S.-I.; Shimizu, F.; Kojima, S. Observation of 90◦ domains in BaTiO3 by atomic force microscopy. Ferroelectrics 1997, 196, 211–214. [CrossRef] Son, J.Y.; Kim, B.Q.; Kim, C.H.; Cho, J.H. Writing polarization bits on the multiferroic BiMnO3 thin film using Kelvin probe force microscope. Appl. Phys. Lett. 2004, 84, 4971–4973. [CrossRef] Hong, J.W.; Park, S.; Khim, Z.G. Measurement of hardness, surface potential, and charge distribution with dynamic contact mode electrostatic force microscope. Rev. Sci. Instrum. 1999, 70, 1735–1739. [CrossRef] Zeng, H.; Shimamura, K.; Kannan, C.V.; Villora, E.G.; Takekawa, S.; Kitamura, K.; Yin, Q. Acoustic imaging of ferroelectric domains in BaTiO3 single crystals using atomic force microscope. Jpn. J. Appl. Phys. 2007, 46, 286–288. [CrossRef] Soergel, E. Piezoresponse force microscopy (PFM). J. Phys. D Appl. Phys. 2011, 44, 464003. [CrossRef] Shur, V.Y.; Zelenovskiy, P.S.; Nebogatikov, M.S.; Alikin, D.O.; Sarmanova, M.F.; Ievlev, A.V.; Mingaliev, E.A.; Kuznetsov, D.K. Investigation of the nanodomain structure formation by piezoelectric force microscopy and Raman confocal microscopy in LiNbO3 and LiTaO3 crystals. J. Appl. Phys. 2011, 110, 052013. [CrossRef] Hooton, J.A.; Merz, W.J. Etch patterns and ferroelectric domains in BaTiO3 single crystals. Phys. Rev. 1955, 98, 409–413. [CrossRef] Cook, W.R. Domain twinning in barium titanate ceramics. J. Am. Ceram. Soc. 1956, 1201, 17–19. [CrossRef] Nassau, K.; Levinstein, H.J.; Loiacono, G.M. The domain structure and etching of ferroelectric lithium niobate. Appl. Phys. Lett. 1965, 6, 228–229. [CrossRef] Kubel, F.; Schmid, H. Growth, twinning and etch figures of ferroelectric/ferroelastic dendritic BiFeO3 single domain crystals. J. Cryst. Growth 1993, 129, 515–524. [CrossRef] Nassau, K.; Levinstein, H.J.; Loiacono, G.M. Ferroelectric lithium niobate. 1. Growth, domain structure, dislocations and etching. J. Phys. Chem. Solids 1966, 27, 983–988. [CrossRef] Belitz, R.K. Differential etching of BaTiO3 by H3 PO4 . J. Am. Ceram. Soc. 1962, 45, 617–618. [CrossRef] Johann, F.; Morelli, A.; Vrejoiu, I. Epitaxial BiFeO3 nanostructures fabricated by differential etching of BiFeO3 films. Appl. Phys. Lett. 2011, 99, 082904. [CrossRef] Rubio-Marcos, F.; del Campo, A.; López-Juárez, R.; Romero, J.J.; Fernández, J.F. High spatial resolution structure of (K, Na)NbO3 lead-free ferroelectric domains. J. Mater. Chem. 2012, 22, 9714–9720. [CrossRef] Arlt, G.; Sasko, P. Domain configuration and equilibrium size of domains in BaTiO3 ceramics. J. Appl. Phys. 1980, 51, 4956–4960. [CrossRef] Ohnishi, N.; Iizuka, T. Etching study of microdomains in LiNbO3 single crystals. J. Appl. Phys. 1975, 46, 1063–1067. [CrossRef] Iwata, M.; Katsuraya, K.; Suzuki, I.; Maeda, M.; Yasuda, N.; Ishibashi, Y. Domain observation of Pb(Zn1/3 Nb2/3 )O3 -PbTiO3 mixed crystals by scanning probe microscopy. Jpn. J. Appl. Phys. 2003, 42, 6201–6204. [CrossRef] Shur, V.Y.; Lobov, A.I.; Shur, A.G.; Kurimura, S.; Nomura, Y.; Terabe, K.; Liu, X.Y.; Kitamura, K. Rearrangement of ferroelectric domain structure induced by chemical etching. Appl. Phys. Lett. 2005, 87, 022905. [CrossRef] Fu, Z.; Yang, J.; Yu, H.; Wang, Y.; Lu, P.; Yang, Q.; Xu, F.; Li, Y. Periodic configuration of δ-boundaries and ferroelectric domains in Li-modified (K, Na)NbO3 lead-free piezoelectric single crystals by solid state crystal growth. J. Eur. Ceram. Soc. 2017, 37, 593–598. [CrossRef] Faryna, M.; Sztwiertnia, K.; Sikorski, K. Simultaneous WDXS and EBSD investigations of dense PLZT ceramics. J. Eur. Ceram. Soc. 2006, 26, 2967–2971. [CrossRef] Uetsuji, Y.; Kuramae, H.; Tsuchiya, K.; Kamlah, M. Electron backscatter diffraction crystal morphology analysis and multiscale simulation of piezoelectric materials. Int. J. Comput. Methods Exp. Meas. 2013, 1, 199–211. [CrossRef] Qin, Y.; Zhang, J.; Gao, Y.; Tan, Y.; Wang, C. Study of domain structure of poled (K,Na)NbO3 ceramics. J. Appl. Phys. 2013, 113, 204107. [CrossRef]

Materials 2017, 10, 47

68. 69.

70. 71. 72. 73.

74.

75.

76. 77. 78. 79. 80. 81. 82. 83.

84.

85. 86.

87. 88.

89.

20 of 23

Ayache, J.; Beaunier, L.; Boumendil, J.; Ehret, G.; Laub, D. Sample Preparation Handbook for Transmission Electron Microscopy; Springer: Berlin, Germany, 2010. Walker, J.; Ursic, H.; Bencan, A.; Malic, B.; Simons, H.; Reaney, I.; Viola, G.; Nagarajan, V.; Rojac, T. Temperature dependent piezoelectric response and strain–electric-field hysteresis of rare-earth modified bismuth ferrite ceramics. J. Mater. Chem. C 2016, 4, 7859–7868. [CrossRef] Guo, H.; Zhang, S.; Beckman, S.P.; Tan, X. Microstructural origin for the piezoelectricity evolution in (K0.5 Na0.5 )NbO3 -based lead-free ceramics. J. Appl. Phys. 2013, 114, 154102. [CrossRef] Fultz, B.; Howe, J. Transmission Electron Microscopy and Diffractometry of Materials, 4th ed.; Springer: New York, NY, USA, 2013. Tsai, F.; Khiznichenko, V.; Cowley, J.M. High-resolution electron microscopy of 90◦ ferroelectric domain boundaries in BaTiO3 and Pb(ZrO0.52 Ti0.48 )O3 . Ultramicroscopy 1992, 45, 55–63. [CrossRef] Gao, P.; Britson, J.; Nelson, C.T.; Jokisaari, J.R.; Duan, C.; Trassin, M.; Baek, S.-H.; Guo, H.; Li, L.; Wang, Y.; et al. Ferroelastic domain switching dynamics under electrical and mechanical excitations. Nat. Commun. 2014, 5, 3801. [CrossRef] [PubMed] Winkler, C.R.; Damodaran, A.R.; Karthik, J.; Martin, L.W.; Taheri, M.L. Direct observation of ferroelectric domain switching in varying electric field regimes using in situ TEM. Micron 2012, 43, 1121–1126. [CrossRef] [PubMed] Kholkin, A.; Kalinin, S.; Roelofs, A.; Gruverman, A. Review of ferroelectric domain imaging by piezoresponse force microscopy. In Scanning Probe Microscopy: Electrical and Electromechanical Phenomena at the Nanoscale; Kalinin, S., Gruverman, A., Eds.; Springer Science & Business Media: New York, NY, USA, 2007; pp. 173–214. Li, J.; Li, J.-F.; Yu, Q.; Chen, Q.N.; Xie, S. Strain-based scanning probe microscopies for functional materials, biological structures, and electrochemical systems. J. Mater. 2015, 1, 3–21. [CrossRef] Kholkin, A.L.; Bdikin, I.K.; Kiselev, D.A.; Shvartsman, V.V.; Kim, S.H. Nanoscale characterization of polycrystalline ferroelectric materials for piezoelectric applications. J. Electroceram. 2007, 19, 81–94. [CrossRef] Jungk, T.; Hoffmann, Á.; Soergel, E. Quantitative analysis of ferroelectric domain imaging with piezoresponse force microscopy. Appl. Phys. Lett. 2006, 89, 163507. [CrossRef] Kalinin, S.; Bonnell, D. Imaging mechanism of piezoresponse force microscopy of ferroelectric surfaces. Phys. Rev. B 2002, 65, 125408. [CrossRef] Kalinin, S.V.; Rodriguez, B.J.; Jesse, S.; Shin, J.; Baddorf, A.P.; Gupta, P.; Jain, H.; Williams, D.B.; Gruverman, A. Vector piezoresponse force microscopy. Microsc. Microanal. 2006, 12, 206–220. [CrossRef] [PubMed] Kim, K.L.; Huber, J.E. Mapping of ferroelectric domain structure using angle-resolved piezoresponse force microscopy. Rev. Sci. Instrum. 2015, 86, 013705. [CrossRef] [PubMed] Cho, Y.; Hashimoto, S.; Odagawa, N.; Tanaka, K.; Hiranaga, Y. Nanodomain manipulation for ultrahigh density ferroelectric data storage. Nanotechnology 2006, 17, 137–141. Pertsev, N.A.; Kiselev, D.A.; Bdikin, I.K.; Kosec, M.; Kholkin, A.L. Quasi-one-dimensional domain walls in ferroelectric ceramics: Evidence from domain dynamics and wall roughness measurements. J. Appl. Phys. 2011, 110, 052001. [CrossRef] Shur, V.Y.; Shishkin, E.I.; Nikolaeva, E.V.; Nebogatikov, M.S.; Alikin, D.O.; Zelenovskiy, P.S.; Sarmanova, M.F.; Dolbilov, M.A. Study of nanoscale domain structure formation using Raman confocal microscopy. Ferroelectrics 2010, 398, 91–97. [CrossRef] DeWolf, I. Micro-Raman spectroscopy to study local mechanical stress in silicon integrated circuits. Semicond. Sci. Technol. 1996, 11, 139–154. Myers, G.A.; Hazra, S.S.; de Boer, M.P.; Michaels, C.A.; Stranick, S.J.; Koseski, R.P.; Cook, R.F.; Delrio, F.W. Stress mapping of micromachined polycrystalline silicon devices via confocal Raman microscopy. Appl. Phys. Lett. 2014, 104, 191908. [CrossRef] Borodavka, F.; Pokorny, J.; Hlinka, J. Combined piezoresponse force microscopy and Raman scattering investigation of domain boundaries in BiFeO3 ceramics. Phase Transit. 2016, 89, 746–751. [CrossRef] Zhao, T.; Scholl, A.; Zavaliche, F.; Lee, K.; Barry, M.; Doran, A.; Cruz, M.P.; Chu, Y.H.; Ederer, C.; Spaldin, N.A.; et al. Electrical control of antiferromagnetic domains in multiferroic BiFeO3 films at room temperature. Nat. Mater. 2006, 5, 823–829. [CrossRef] [PubMed] Shvartsman, V.V.; Kleemann, W.; Haumont, R.; Kreisel, J. Large bulk polarization and regular domain structure in ceramic BiFeO3 . Appl. Phys. Lett. 2007, 90, 172115. [CrossRef]

Materials 2017, 10, 47

90.

91.

92.

93.

94. 95. 96.

97.

98. 99.

100. 101. 102. 103. 104. 105. 106. 107. 108. 109.

110. 111.

21 of 23

Tu, C.S.; Chen, C.S.; Chen, P.Y.; Wei, H.H.; Schmidt, V.H.; Lin, C.Y.; Anthoniappen, J.; Lee, J.M. Enhanced photovoltaic effects in A-site samarium doped BiFeO3 ceramics: The roles of domain structure and electronic state. J. Eur. Ceram. Soc. 2016, 36, 1149–1157. [CrossRef] Alikin, D.O.; Turygin, A.P.; Walker, J.; Bencan, A.; Malic, B.; Rojac, T.; Shur, V.Y.; Kholkin, A.L. The effect of phase assemblages, grain boundaries and domain structure on the local switching behavior of rare-earth modified bismuth ferrite ceramics. Acta Mater. 2017, 125, 265–273. [CrossRef] Walker, J.; Simons, H.; Alikin, D.O.; Turygin, A.P.; Shur, V.Y.; Kholkin, A.L.; Ursic, H.; Bencan, A.; Malic, B.; Nagarajan, V.; et al. Dual strain mechanisms in a lead-free morphotropic phase boundary ferroelectric. Sci. Rep. 2016, 6, 19630. [CrossRef] [PubMed] Alikin, D.O.; Turygin, A.P.; Walker, J.; Rojac, T.; Shvartsman, V.V.; Shur, V.Y.; Kholkin, A.L. Quantitative phase separation in multiferroic Bi0.88 Sm0.12 FeO3 ceramics via piezoresponse force microscopy. J. Appl. Phys. 2015, 118. [CrossRef] Karimi, S.; Reaney, I.M.; Levin, I.; Sterianou, I. Nd-doped BiFeO3 ceramics with antipolar order. Appl. Phys. Lett. 2009, 94, 112903. [CrossRef] Karimi, S.; Reaney, I.M.; Han, Y.; Pokorny, J.; Sterianou, I. Crystal chemistry and domain structure of rare-earth doped BiFeO3 ceramics. J. Mater. Sci. 2009, 44, 5102–5112. [CrossRef] Khomchenko, V.A.; Karpinsky, D.V.; Kholkin, A.L.; Sobolev, N.A.; Kakazei, G.N.; Araujo, J.P.; Troyanchuk, I.O.; Costa, B.F.O.; Paixão, J.A. Rhombohedral-to-orthorhombic transition and multiferroic properties of Dy-substituted BiFeO3 . J. Appl. Phys. 2010, 108, 074109. [CrossRef] Khomchenko, V.A.; Kiselev, D.A.; Bdikin, I.K.; Shvartsman, V.V.; Borisov, P.; Kleemann, W.; Vieira, J.M.; Kholkin, A.L. Crystal structure and multiferroic properties of Gd-substituted BiFeO3 . Appl. Phys. Lett. 2008, 93, 262905. [CrossRef] Yan, F.; Xing, G.; Wang, R.; Li, L. Tailoring surface phase transition and magnetic behaviors in BiFeO3 via doping engineering. Sci. Rep. 2015, 5, 9128. [CrossRef] [PubMed] Coondoo, I.; Panwar, N.; Bdikin, I.; Puli, V.S.; Katiyar, R.S.; Kholkim, A.L. Structural, morphological and piezoresponse studies of Pr and Sc co-substituted BiFeO3 ceramics. J. Phys. D Appl. Phys. 2012, 45, 055302. [CrossRef] Azough, F.; Freer, R.; Thrall, M.; Cernik, R.; Tuna, F.; Collison, D. Microstructure and properties of Co-, Ni-, Zn-, Nb- and W-modified multiferroic BiFeO3 ceramics. J. Eur. Ceram. Soc. 2010, 30, 727–736. [CrossRef] Khesro, A.; Boston, R.; Sterianou, I.; Sinclair, D.C.; Reaney, I.M. Phase transitions, domain structure, and pseudosymmetry in La- and Ti-doped BiFeO3 . J. Appl. Phys. 2016, 119, 054101. [CrossRef] Chen, J.; Cheng, J.; Dong, S. Review on high temperature piezoelectric ceramics and actuators based on BiScO3 –PbTiO3 solid solutions. J. Adv. Dielectr. 2014, 4, 1430002. [CrossRef] De, U.; Sahu, K.R.; De, A. Ferroelectric materials for high temperature piezoelectric applications. Solid State Phenom. 2015, 232, 235–278. [CrossRef] Rojac, T.; Makarovic, M.; Walker, J.; Ursic, H.; Damjanovic, D.; Kos, T. Piezoelectric response of BiFeO3 ceramics at elevated temperatures. Appl. Phys. Lett. 2016, 109, 042904. [CrossRef] Martirena, H.T.; Burfoot, J.C. Grain-size effects on properties of some ferroelectric ceramics. J. Phys. C Solid State Phys. 1974, 7, 3182–3192. [CrossRef] Pertsev, N.A.; Arlt, G. Theory of the banded domain structure in coarse-grained ferroelectric ceramics. Ferroelectrics 1992, 132, 27–40. [CrossRef] Arlt, G.; Hennings, D.; de With, G. Dielectric properties of fine-grained barium titanate ceramics. J. Appl. Phys. 1985, 58, 1619–1625. [CrossRef] Arlt, G.; Dederichs, H. Complex elastic, dielectric and piezoelectric constants by domain wall damping in ferroelectric ceramics. Ferroelectrics 1980, 29, 47–50. [CrossRef] Escobar Castillo, M.; Shvartsman, V.V.; Gobeljic, D.; Gao, Y.; Landers, J.; Wende, H.; Lupascu, D.C. Effect of particle size on ferroelectric and magnetic properties of BiFeO3 nanopowders. Nanotechnology 2013, 24, 355701. [CrossRef] [PubMed] López, R.; González, F.; Cruz, M.P.; Villafuerte-Castrejon, M.E. Piezoelectric and ferroelectric properties of K0.5 Na0.5 NbO3 ceramics synthesized by spray drying method. Mater. Res. Bull. 2011, 46, 70–74. [CrossRef] Qin, Y.; Zhang, J.; Yao, W.; Wang, C.; Zhang, S. Domain structure of potassium-sodium niobate ceramics before and after poling. J. Am. Ceram. Soc. 2015, 98, 1027–1033. [CrossRef]

Materials 2017, 10, 47

22 of 23

112. Fousek, J.; Janovec, V. The orientation of domain walls in twinned ferroelectric crystals. J. Appl. Phys. 1969, 40, 135–142. [CrossRef] 113. Hirohashi, J.; Yamada, K.; Kamio, H.; Uchida, M.; Shichijyo, S. Control of specific domain structure in KNbO3 single crystals by differential vector poling method. J. Appl. Phys. 2005, 98, 034107. [CrossRef] 114. Zhang, J.; Tian, X.; Gao, Y.; Yao, W.; Qin, Y.; Su, W. Domain structure of poled (K0.50 Na0.50 )1-x Lix NbO3 ceramics with different stabilities. J. Am. Ceram. Soc. 2015, 98, 990–995. [CrossRef] 115. Cho, J.; Lee, Y.; Han, K.; Chun, M.; Nam, J.; Kim, B. Effect of domain size on the coercive field of orthorhombic (Li, K, Na)NbO3 ceramics. J. Korean Phys. Soc. 2010, 57, 971–974. 116. Arlt, G. Twinning in ferroelectric and ferroelastic ceramics: Stress relief. J. Mater. Sci. 1990, 25, 2655–2666. [CrossRef] 117. Rossetti, G.A.; Khachaturyan, A.G.; Akcay, G.; Ni, Y. Ferroelectric solid solutions with morphotropic boundaries: Vanishing polarization anisotropy, adaptive, polar glass, and two-phase states. J. Appl. Phys. 2008, 103, 114113. [CrossRef] 118. Cho, J.; Lee, Y.; Kim, B. Domain structure of orthorhombic (Li,K,Na)NbO3 ceramics. J. Ceram. Process. Res. 2010, 11, 237–240. 119. Rai, R.; Coondoo, I.; Rani, R.; Bdikin, I.; Sharma, S.; Kholkin, A.L. Impedance spectroscopy and piezoresponse force microscopy analysis of lead-free (1 − x) K0.5 Na0.5 NbO3 − xLiNbO3 ceramics. Curr. Appl. Phys. 2013, 13, 430–440. [CrossRef] 120. Qin, Y.; Zhang, J.; Tan, Y.; Yao, W.; Wang, C.; Zhang, S. Domain configuration and piezoelectric properties of (K0.50 Na0.50 )1−x Lix (Nb0.80 Ta0.20 )O3 ceramics. J. Eur. Ceram. Soc. 2014, 34, 4177–4184. [CrossRef] 121. Fu, J.; Zuo, R.; Xu, Z.; Fu, J.; Zuo, R.; Xu, Z. High piezoelectric activity in (Na,K)NbO3 based lead-free piezoelectric ceramics: Contribution of nanodomains. Appl. Phys. Lett. 2011, 99, 21–24. [CrossRef] 122. Gobeljic, D.; Shvartsman, V.V.; Wang, K.; Yao, F.; Li, J.F.; Jo, W.; Rödel, J.; Lupascu, D.C. Temperature dependence of the local piezoresponse in (K,Na)NbO3 -based ceramics with large electromechanical strain. J. Appl. Phys. 2014, 116, 066811. [CrossRef] 123. Zhou, J.; Wang, K.; Yao, F.; Zheng, T.; Wu, J. Lead-free piezoceramics with large piezoelectricity. J. Mater. Chem. C 2015, 3, 8780–8787. [CrossRef] 124. Hunpratub, S.; Yamwong, T.; Srilomsak, S.; Maensiri, S.; Chindaprasirt, P. Effect of particle size on the dielectric and piezoelectric properties of 0–3BCTZO/cement composites. Ceram. Int. 2015, 40, 1209–1213. [CrossRef] 125. Taylor, D.V.; Damjanovic, D. Evidence of domain wall contribution to the dielectric permittivity in PZT thin films at sub-switching fields. J. Appl. Phys. 1997, 82, 1973–1975. [CrossRef] 126. Rayleigh, L. XXV. Notes on electricity and magnetism.—III. On the behaviour of iron and steel under the operation of feeble magnetic forces. Philos. Mag. Ser. 5 1887, 23, 225–245. [CrossRef] 127. Pertsev, N.A.; Arlt, G. Forced translational vibrations of 90-degrees domain-walls and the dielectric-dispersion in ferroelectric ceramics. J. Appl. Phys. 1993, 74, 4105–4112. [CrossRef] 128. Buixaderas, E.; Bovtun, V.; Kempa, M.; Savinov, M.; Nuzhnyy, D.; Kadlec, F.; Vanˇek, P.; Petzelt, J.; Eriksson, M.; Shen, Z. Broadband dielectric response and grain-size effect in K0.5 Na0.5 NbO3 ceramics. J. Appl. Phys. 2010, 107, 014111. [CrossRef] 129. Zhang, J.; Hao, W.; Gao, Y.; Qin, Y.; Tan, Y.; Wang, C. Large decrease of characteristic frequency of dielectric relaxation associated with domain-wall motion in Sb5+ -modified (K, Na)NbO3 -based ceramics. Appl. Phys. Lett. 2012, 101, 252905. [CrossRef] 130. Huan, Y.; Wang, X.; Li, L.; Koruza, J. Strong domain configuration dependence of the nonlinear dielectric response in (K, Na)NbO3 -based ceramics. Appl. Phys. Lett. 2015, 107, 202903. [CrossRef] 131. Arlt, G.; Böttger, U.; Witte, S. Dielectric dispersion of ferroelectric ceramics and single crystals at microwave frequencies. Ann. Phys. 1994, 506, 578–588. [CrossRef] 132. Arlt, G.; Böttger, U.; Witte, S. Emission of GHz shear waves by ferroelastic domain walls in ferroelectrics. Appl. Phys. Lett. 1993, 63, 602–604. [CrossRef] 133. Arlt, G.; Böttger, U.; Witte, S. Dielectric dispersion of ferroelectric ceramics and single crystals by sound generation in piezoelectric domains. J. Am. Ceram. Soc. 1995, 78, 1097–1100. [CrossRef] 134. Arlt, G. Strong ultrasonic microwaves in ferroelectric ceramics. IEEE Trans. Ultrason. Ferroelectr. Freq. Control 1998, 45, 4–10. [CrossRef] [PubMed]

Materials 2017, 10, 47

23 of 23

135. Arlt, G. Ferroelastic domain walls as powerful shear wave emitters at microwaves. Ferroelectrics 1995, 172, 95–104. [CrossRef] 136. Rojac, T.; Ursic, H.; Bencan, A.; Malic, B.; Damjanovic, D. Mobile domain walls as a bridge between nanoscale conductivity and macroscopic electromechanical response. Adv. Funct. Mater. 2015, 25, 2099–2108. [CrossRef] 137. Peng, B.; Yue, Z.; Li, L.; Peng, B.; Yue, Z.; Li, L. Evaluation of domain wall motion during polymorphic phase transition in (K, Na)NbO3 —Based piezoelectric ceramics by nonlinear response measurements. J. Appl. Phys. 2016, 109, 54107. [CrossRef] 138. Morozov, M.I.; Hoffmann, M.J.; Benkert, K.; Schuh, C. Influence of the A/B nonstoichiometry, composition modifiers, and preparation methods on properties of Li- and Ta-modified (K, Na)NbO3 ceramics. J. Appl. Phys. 2012, 112, 114107. [CrossRef] 139. Ge, J.; Qi, Y.; Wan, Y.; Tang, Y.; Wang, F.; Pan, Y.; Jiang, B.; Sun, D. Relationship between dielectric property and shrinkage in ferroelectric KNN-based ceramics. Ferroelectrics 2014, 458, 37–41. [CrossRef] 140. Huan, Y.; Wang, X.; Koruza, J.; Wang, K.; Webber, K.G.; Hao, Y.; Li, L. Inverted electro-mechanical behaviour induced by the irreversible domain configuration transformation in (K, Na)NbO3 —Based ceramics. Sci. Rep. 2016, 6, 22053. [CrossRef] [PubMed] 141. Huan, Y.; Wang, X.; Shen, Z.; Kim, J.; Zhou, H.; Li, L. Nanodomains in KNN-based lead-free piezoelectric ceramics: Origin of strong piezoelectric properties. J. Am. Ceram. Soc. 2014, 97, 701–703. [CrossRef] 142. Hayati, R.; Barzegar, A. Microstructure and electrical properties of lead free potassium sodium niobate piezoceramics with nano ZnO additive. Mater. Sci. Eng. B 2010, 172, 121–126. [CrossRef] © 2017 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC-BY) license (http://creativecommons.org/licenses/by/4.0/).