Characterization of Pb-free KNbO 3-and ( Na , Bi ) TiO 3-based piezoelectric single-crystals using X-ray and neutron diffraction

In view of serious environmental concerns with traditional Pb-based piezoelectrics, the search for new Pb-free alternatives has intensified recently. A thorough investigation of structure-property relationships in Pb-free piezoelectrics is desired in order to design new material compositions with high electromechanical properties that can be operated over a broader range of conditions. Recent availability of high-quality single crystals has not only opened the possibility for achieving multifold enhancements in the electromechanical properties in this new class of materials, but has also provided opportunities to undertake fundamental studies on their structure-property relationships. In the following pages, I review some of the recent X-ray and neutron diffraction studies of Pb-free piezoelectric single crystals, which have provided novel insights into their multiscale stimuli-induced structural mechanisms, including phase transitions, nanoscale structural ordering, lattice instability, and domain re-orientation. Opportunities and challenges for future progress in this area of study are discussed.


Introduction
Piezoelectricity is the coupling between dielectric and elastic displacements in a material that are induced under the effects of an applied mechanical stress or electric field [1]. The piezoelectric effect is necessarily observed in materials with non-centrosymmetric crystal structures. Large piezoelectric effects are observed in ferroelectric materials, for which each crystallographic unit cell also exhibits a spontaneous and reversible polarization or dipole moment under practical conditions [2]. Owing to the large electro-mechanical coupling phenomena in ferroelectric materials, they are of interest for a broad array of applications including precision actuators, impact and load sensors, telecommunication, and biomedical imaging [3][4][5]. With growing innovation in sectors such as internet of things (IoT), mechatronics, and decentralized energy harvesting, a large jump in the usage of ferroelectric materials as smart electromechanical transducers can be foreseen in the future.
Traditionally, the industry for ferroelectrics in electromechanical applications has been dominated by materials with high Pb-content [6]. For example, the most popular industrial piezoelectrics, such as Pb(Zr x Ti 1−x )O 3 (PZT) and PbMg 1/3 Nb 2/3 O 3 -PbTiO 3 (PMN-PT), contain up to 60% Pb by weight. Pb is known to be toxic to the functioning of biological life-forms. Although Pb is usually stable in oxide form, it can eventually leak into the environment upon disposal over the long-term, and therefore poses danger to the ecosystem. Consequently, a larger scale adoption of Pb-based ferroelectrics is foreseen to pose serious environmental danger. With such foresight, more stringent regulations on the use of Pb in electronic materials are being imposed by several agencies worldwide [7][8][9][10]. In view of this, the development of Pb-free ferroelectrics with electromechanical properties at par with or exceeding those of Pb-based ferroelectrics is of utmost importance [11]. However, a limiting factor for broader investigation of Pb-free ferroelectric single crystals could be difficulties in synthesizing crystals of adequate dimensions and quality. In this respect, recent advancements on synthesis of good quality Pb-free piezoelectric single crystals of complex compositions with high functional properties is an exciting development. Notable mentions in this category include crystals of Li-doped potassium sodium niobate using the Bridgman method [24], (1−x)BaTiO3−xCaTiO3 using floating zone technique [25], (Ba(1−x),Cax)(Ti(1−y)Zry)O3 by Chochralski method [26], (Ba(1−x),Cax)(Ti(1−y)Zry)O3 by the flux-method [27], (Na0.5Bi0.5)TiO3-BaTiO3 and (K,Na)(Nb,Ta)O3 using the top-seeded solution growth (TSSG) method [23,[28][29][30], and (Na0.5Bi0.5)TiO3-BaTiO3-and (K,Na)(Nb)O3-based crystals using the seed-free solid state crystal growth technique [31][32][33]. The developments with regard to solid state crystal growth are particularly noteworthy, since they provide a cheaper alternative to growing crystals from melt. An exhaustive review on the art of converting polycrystalline ceramics to single-crystalline form for Pb-free ferroelectrics is given by Kang et al. [34].

Characterization of Multiscale Structural Mechanisms Using Diffraction
The above mentioned recent developments towards synthesis of Pb-free ferroelectric single crystals are not only exciting from the perspective of achieving higher properties, but they can also facilitate better identification of the fundamental physical mechanisms responsible for electromechanical responses. The macroscopic electromechanical responses in a ferroelectric material, including both stress-induced polarization and electric-field-induced strains, are the result of structural changes that are operative at multiple length scales-these include atomic displacements and distortion of crystallographic unit cell, short-range ordering within nanoscale polar clusters, and motion of ferroelectric/ferroelastic domain boundaries ( Figure 2) [5,35,36]. A clear identification of the prevalence of each of these structural mechanisms in materials of different compositions and their respective roles towards enhancement of macroscopic properties can lead to faster development of new Pb-free ferroelectrics. High-energy X-ray and neutron diffraction (or scattering in general) experiments of materials under applied external stimuli has been a preferred method of investigation for identification of microscopic structural mechanisms in a wide variety of materials. Reasons for such a preference include easy sample preparation, the option of gaining structural information from bulk material volume, ease of designing sample environments, and relatively faster data collection. In situ diffraction experiments have been broadly applied in the past However, a limiting factor for broader investigation of Pb-free ferroelectric single crystals could be difficulties in synthesizing crystals of adequate dimensions and quality. In this respect, recent advancements on synthesis of good quality Pb-free piezoelectric single crystals of complex compositions with high functional properties is an exciting development. Notable mentions in this category include crystals of Li-doped potassium sodium niobate using the Bridgman method [24], (1−x)BaTiO 3−x CaTiO 3 using floating zone technique [25], (Ba (1−x) ,Ca x )(Ti (1−y) Zr y )O 3 by Chochralski method [26], (Ba (1−x) ,Ca x )(Ti (1−y) Zr y )O 3 by the flux-method [27], (Na 0.5 Bi 0.5 )TiO 3 -BaTiO 3 and (K,Na)(Nb,Ta)O 3 using the top-seeded solution growth (TSSG) method [23,[28][29][30], and (Na 0.5 Bi 0.5 )TiO 3 -BaTiO 3and (K,Na)(Nb)O 3 -based crystals using the seed-free solid state crystal growth technique [31][32][33]. The developments with regard to solid state crystal growth are particularly noteworthy, since they provide a cheaper alternative to growing crystals from melt. An exhaustive review on the art of converting polycrystalline ceramics to single-crystalline form for Pb-free ferroelectrics is given by Kang et al. [34].

Characterization of Multiscale Structural Mechanisms Using Diffraction
The above mentioned recent developments towards synthesis of Pb-free ferroelectric single crystals are not only exciting from the perspective of achieving higher properties, but they can also facilitate better identification of the fundamental physical mechanisms responsible for electromechanical responses. The macroscopic electromechanical responses in a ferroelectric material, including both stress-induced polarization and electric-field-induced strains, are the result of structural changes that are operative at multiple length scales-these include atomic displacements and distortion of crystallographic unit cell, short-range ordering within nanoscale polar clusters, and motion of ferroelectric/ferroelastic domain boundaries ( Figure 2) [5,35,36]. A clear identification of the prevalence of each of these structural mechanisms in materials of different compositions and their respective roles towards enhancement of macroscopic properties can lead to faster development of new Pb-free ferroelectrics. High-energy X-ray and neutron diffraction (or scattering in general) experiments of materials under applied external stimuli has been a preferred method of investigation for identification of microscopic structural mechanisms in a wide variety of materials. Reasons for such a preference include easy sample preparation, the option of gaining structural information from bulk material volume, ease of designing sample environments, and relatively faster data collection. In situ diffraction experiments have been broadly applied in the past to identify the microscopic structural mechanisms operative under external electric fields or stress in both single-crystal and polycrystalline forms of Pb-based ferroelectrics [37]. In particular, diffraction experiments in single-crystalline materials were instrumental in identifying subtler changes in materials structure, which are not readily observable from experiments with polycrystalline ceramics, such as small distortions from ground-state crystal symmetry, changes in nanoscale ordering, and domain reorientation pathways [38][39][40]. to identify the microscopic structural mechanisms operative under external electric fields or stress in both single-crystal and polycrystalline forms of Pb-based ferroelectrics [37]. In particular, diffraction experiments in single-crystalline materials were instrumental in identifying subtler changes in materials structure, which are not readily observable from experiments with polycrystalline ceramics, such as small distortions from ground-state crystal symmetry, changes in nanoscale ordering, and domain reorientation pathways [38][39][40]. Recent advancements with regard to available flux and resolution for diffraction experiments has broadened the scope and precision of structural insights that can be gained from in situ experiments. Indeed, advanced diffraction experiments on polycrystalline ceramics have elucidated several electric field and/or stress-induced structural mechanisms that are unique to Pb-free piezoelectrics. A review of these developments can be found in references [41,42]. Nevertheless, similar reviews on single-crystal diffraction experiments of Pb-free piezoelectrics are absent. The author is of the opinion that such a review on the study of structure-property correlations of Pb-free piezoelectrics using single-crystal diffraction is warranted in view of recent developments and potential opportunities.

Crystallographic Phase Transition
External stimuli induced changes in the original crystal structure of a ferroelectric crystal constitutes the intrinsic component of its electromechanical response. Without consideration to local structural disorders, this involves changes in the relative atomic positions within the unit cell that defines the average long-range crystallographic structure. It is debatable whether the alteration in atomic arrangements constitutes a structural distortion or true phase transition, which are different for specific cases [20]. Nevertheless, it is important to identify such structural mechanisms to understand material behavior from a fundamental perspective and to explore possible avenues for enhancement in electromechanical properties. Diffraction is widely used as an experimental tool for characterization of stimuli-induced crystallographic phase transitions (or distortions) in ferroelectrics. The crystallographic structural changes can be characterized from data collected in a diffraction experiment by monitoring changes in some essential features, such as hkl diffraction peak position to determine changes in lattice planes spacing and splitting of hkl reflections to identify changes in crystallographic symmetry. In addition, change in intensities of hkl reflections can reflect various crystallographic and microstructural changes such as relative positions of different ions, Recent advancements with regard to available flux and resolution for diffraction experiments has broadened the scope and precision of structural insights that can be gained from in situ experiments. Indeed, advanced diffraction experiments on polycrystalline ceramics have elucidated several electric field and/or stress-induced structural mechanisms that are unique to Pb-free piezoelectrics. A review of these developments can be found in references [41,42]. Nevertheless, similar reviews on single-crystal diffraction experiments of Pb-free piezoelectrics are absent. The author is of the opinion that such a review on the study of structure-property correlations of Pb-free piezoelectrics using single-crystal diffraction is warranted in view of recent developments and potential opportunities.

Crystallographic Phase Transition
External stimuli induced changes in the original crystal structure of a ferroelectric crystal constitutes the intrinsic component of its electromechanical response. Without consideration to local structural disorders, this involves changes in the relative atomic positions within the unit cell that defines the average long-range crystallographic structure. It is debatable whether the alteration in atomic arrangements constitutes a structural distortion or true phase transition, which are different for specific cases [20]. Nevertheless, it is important to identify such structural mechanisms to understand material behavior from a fundamental perspective and to explore possible avenues for enhancement in electromechanical properties. Diffraction is widely used as an experimental tool for characterization of stimuli-induced crystallographic phase transitions (or distortions) in ferroelectrics. The crystallographic structural changes can be characterized from data collected in a diffraction experiment by monitoring changes in some essential features, such as hkl diffraction peak position to determine changes in lattice planes spacing and splitting of hkl reflections to identify changes in crystallographic symmetry. In addition, change in intensities of hkl reflections can reflect various crystallographic and microstructural changes such as relative positions of different ions, state of microstructural defects, as well as internal strains, which can be distinguished based on quantitative modeling [36,37,[41][42][43][44].
The scattered intensity is essentially a Fourier transformation of the microscopic arrangement of the atoms inside a material. This can be represented as, where N is the total number of atoms and f m is the scattering factor of an atom m that is located at the lattice site at the location → R m . The scattering factor is different for different elements, depending on the interaction of the X-rays with the electron clouds of the atoms or the interaction of the neutrons with the atomic nuclei. It can be readily seen from Equation (1) that the intensities of the diffracted beams are very sensitive to the arrangement of atoms in a crystal. The diffraction pattern from polycrystalline ceramics constitutes a summation of scattering intensities from different oriented crystallites within the sample, which is presented either as a 1-dimensional line profile with multiple hkl diffraction peaks or concentric circles of hkl reflections, depending on specific experimental geometry. Diffraction from polycrystalline ceramic samples are often adequate to identify and characterize many of changes as listed above. However, in other cases, single-crystal diffraction experiments are necessary for characterization of subtler changes in crystal symmetry, which can otherwise remain elusive. For example, high-resolution single-crystal diffraction has been particularly helpful in clearly identifying monoclinic structural distortions and polarization reorientation mechanism in Pb-based relaxor ferroelectric materials [38]. Although such experiments are not yet commonplace for Pb-free ferroelectric materials, examples from some recent studies are reviewed below.
Kreisel et al. used high-energy X-ray diffraction of Na 0.5 Bi 0.5 TiO 3 single crystals to show a series of phase transitions induced under the application of mechanical pressure which were identified from splitting of diffraction peaks, such as shown in Figure 3 [45]. A rhombohedral-to-monoclinic transition was observed at lower pressures, while a further transition to tetragonal phase could be observed under a pressure of 11 GPa. Although a similar phase transition was also shown from diffraction measurements of polycrystalline ceramic samples, single-crystal diffraction provided a more definitive description of structural mechanisms involving correlated symmetry changes and octahedral tilting.
In yet another interesting experiment, Ogino et al. determined the polarization rotation pathway for the average monoclinic crystal structure in a Na 0.5 Bi 0.5 TiO 3 -BaTiO 3 single crystal sample [30]. Based on detailed analysis of diffraction spot positions, and further assisted by DFT calculations, the authors showed that the polarization vector rotates by ∼2 • in the monoclinic a-c plane when an electric field E of 70 kV/cm is applied along <001> pseudo-cubic direction. state of microstructural defects, as well as internal strains, which can be distinguished based on quantitative modeling [36,37,[41][42][43][44]. The scattered intensity is essentially a Fourier transformation of the microscopic arrangement of the atoms inside a material. This can be represented as, where N is the total number of atoms and fm is the scattering factor of an atom m that is located at the lattice site at the location m R  . The scattering factor is different for different elements, depending on the interaction of the X-rays with the electron clouds of the atoms or the interaction of the neutrons with the atomic nuclei. It can be readily seen from Equation (1) that the intensities of the diffracted beams are very sensitive to the arrangement of atoms in a crystal. The diffraction pattern from polycrystalline ceramics constitutes a summation of scattering intensities from different oriented crystallites within the sample, which is presented either as a 1-dimensional line profile with multiple hkl diffraction peaks or concentric circles of hkl reflections, depending on specific experimental geometry. Diffraction from polycrystalline ceramic samples are often adequate to identify and characterize many of changes as listed above. However, in other cases, single-crystal diffraction experiments are necessary for characterization of subtler changes in crystal symmetry, which can otherwise remain elusive. For example, high-resolution single-crystal diffraction has been particularly helpful in clearly identifying monoclinic structural distortions and polarization reorientation mechanism in Pb-based relaxor ferroelectric materials [38]. Although such experiments are not yet commonplace for Pb-free ferroelectric materials, examples from some recent studies are reviewed below.
Kreisel et al. used high-energy X-ray diffraction of Na0.5Bi0.5TiO3 single crystals to show a series of phase transitions induced under the application of mechanical pressure which were identified from splitting of diffraction peaks, such as shown in Figure 3 [45]. A rhombohedral-to-monoclinic transition was observed at lower pressures, while a further transition to tetragonal phase could be observed under a pressure of 11 GPa. Although a similar phase transition was also shown from diffraction measurements of polycrystalline ceramic samples, single-crystal diffraction provided a more definitive description of structural mechanisms involving correlated symmetry changes and octahedral tilting.
In yet another interesting experiment, Ogino et al. determined the polarization rotation pathway for the average monoclinic crystal structure in a Na0.5Bi0.5TiO3-BaTiO3 single crystal sample [30]. Based on detailed analysis of diffraction spot positions, and further assisted by DFT calculations, the authors showed that the polarization vector rotates by ∼2° in the monoclinic a-c plane when an electric field E of 70 kV/cm is applied along <001> pseudo-cubic direction.

Nanoscale Structural Disorders
In many ferroelectrics, the optimum electromechanical properties often correlate with the presence of nanoscale structural disorders. This means that while translational periodicity of the crystal lattice is maintained on average, there are local deviations from the average structure because of energetically degenerate quasi-random structural features such as atomic displacements, octahedral tilts, and/or elastic distortions. Instead of a perfect long-range order, these structural features are correlated only over short length scales within nanoscale domains. The presence of such nanoscale domains and their relation to large functional properties has been reasonably explained for Pb-based ferroelectrics, although there still exist significant knowledge gaps [5,46]. Recent investigations are also bringing forth the potential central role that such nanoscale domains can play towards enhancing the electromechanical properties of new Pb-free ferroelectrics. One promising strategy could be the engineering of local structural disorder in such a way so as to induce large structural fluctuations under critical operating conditions. However, in order to successfully execute such an approach, it is a pre-requisite that the types of structural disorder in different Pb-free ferroelectrics, as well as how they respond to external stimuli such as electric fields or mechanical stress, are clearly understood. It needs to be mentioned here that due to differences in chemical bonding environments, the types of local structural disorders in the new Pb-free ferroelectrics are expected to be substantially different than those studied earlier for Pb-based ferroelectrics. For example, A-site displacements due to a more covalent Pb-O (A-O) bond drive ferroelectric distortion in Pb-based ferroelectrics. However, Ba-O and K-O bond has a more ionic character and therefore, <111> type Ti and Nb displacements are the dominant factors for ferroelectricity in BaTiO 3 and KNbO 3 [47,48]. It is thus warranted that external-stimuli-induced changes in nanoscale atomic ordering are investigated anew in Pb-free ferroelectrics.
The measurement and analyses of experimental signatures from nanoscale correlated disorders in materials are very different from those used in classical crystallography to determine long-range ordered structures [49]. Using tools of classical crystallography, the long-range atomic structure can be determined accurately by noting the position and the area of well-marked intense peaks in the diffraction pattern known as Bragg peaks. Instead, the correlated atomic disorder within nanoscale domains needs to be determined from the diffuse scattered intensities in between the Bragg peaks [49][50][51]. Diffraction experiments of high-quality single crystals enables measurement of such diffuse scattering signatures in reciprocal space in much more details than what can be achieved in polycrystalline powder samples. In this regard, a brief introduction to the origin of structured diffuse scattering for single-crystal diffraction is provided below.
In the presence of local disorders, the above Equation (1) can be modified as where → u m and → u n refer to small displacements from the average positions of → R m and → R n [52,53]. Essentially, Equation (2) takes into account both chemical disorder on different atomic sites ( f m and f n ), as well as disorder in atomic displacements. If such disorder occurs in a truly random fashion, then it results in incoherent scattering, which contributes to the background. However, if there is short-range order present, either with regard to atomic site occupancies or atomic displacements, then well-defined diffuse scattering patterns are observed in addition to the Bragg peaks. Here, it can be seen that the diffuse scattering is also a function of → Q or the diffraction wave vector. Therefore, in order to better model the nanoscale ordered structures that give rise to patterned diffuse scattering within the different hkl reciprocal lattice planes, measurement of the same should be taken with clear distinction of → Q direction. This is only achievable with diffraction from single-crystal samples.
Below, I have reviewed some recent single-crystal diffraction studies in Pb-free ferroelectrics that have shed light on the kind of structural disorders present in these materials. In addition, examples of in situ experiments are provided that are enabling real-time investigation of how local structures evolve under external stimuli.

Atomic Displacement Disorder
Although not readily appreciated, atomic displacement disorder is broadly prevalent in solid solutions of Pb-free ferroelectrics such as BaTiO 3 and KNbO 3 . These two prototypical compounds had been traditionally described as belonging to the displacive type ferroelectrics, in which the equilibrium positions of the B-site cation shifts off-center along a particular direction which matches with the average polarization for the respective phase, such as [001] in the tetragonal phase [2]. However, subsequent investigations overwhelmingly demonstrated that these exhibit characteristics partial to both displacive and order-disorder type ferroelectrics.
According to the order-disorder model, originally proposed by Comes et al., the B-site cation in BaTiO 3 and KNbO 3 is displaced in all phases along the <111> body diagonal directions, as illustrated in Figure 4 [54,55]. In the high-temperature cubic phase, displacements of B-site cation along all the <111> directions are equally probable, leading to no net polarization. With decreasing temperature, the degeneracy of the B-site cation displacements along the different body-diagonals is broken. In the ferroelectric tetragonal phase, B-site cation displacements are favored for only 4 of the 8 <111> directions, thereby leading to a net electrical polarization along [001] direction. At still lower temperatures in the ferroelectric orthorhombic phase, only 2 of the 8 <111> directions are favored, leading to a net polarization along [110]. In the ferroelectric rhombohedral ground-state phase at even lower temperatures, the structure is fully ordered with all B-site cation displacements along single [111] direction. The different correlations among the <111> B-atom displacements can be determined from the diffuse scattering patterns. 1-dimensional correlation of the <111> B-atom displacements along <100> crystallographic directions gives rise to sheets of diffuse scattering intensity in the {100} reciprocal lattice planes, such as shown in Figure 4c. Based on the types of atomic displacements in the different average crystallographic phase, one can therefore expect to have three sheets of diffuse scattering in the cubic phase, two such sheets in the tetragonal phase, and only one such sheet in the orthorhombic phase; since the atomic displacements are correlated over long-distances in the rhombohedral phase, no such diffuse scattering sheets are visible.

Atomic Displacement Disorder
Although not readily appreciated, atomic displacement disorder is broadly prevalent in solid solutions of Pb-free ferroelectrics such as BaTiO3 and KNbO3. These two prototypical compounds had been traditionally described as belonging to the displacive type ferroelectrics, in which the equilibrium positions of the B-site cation shifts off-center along a particular direction which matches with the average polarization for the respective phase, such as [001] in the tetragonal phase [2]. However, subsequent investigations overwhelmingly demonstrated that these exhibit characteristics partial to both displacive and order-disorder type ferroelectrics.
According to the order-disorder model, originally proposed by Comes et al., the B-site cation in BaTiO3 and KNbO3 is displaced in all phases along the <111> body diagonal directions, as illustrated in Figure 4 [54,55]. In the high-temperature cubic phase, displacements of B-site cation along all the <111> directions are equally probable, leading to no net polarization. With decreasing temperature, the degeneracy of the B-site cation displacements along the different body-diagonals is broken. In the ferroelectric tetragonal phase, B-site cation displacements are favored for only 4 of the 8 <111> directions, thereby leading to a net electrical polarization along [001] direction. At still lower temperatures in the ferroelectric orthorhombic phase, only 2 of the 8 <111> directions are favored, leading to a net polarization along [110]. In the ferroelectric rhombohedral ground-state phase at even lower temperatures, the structure is fully ordered with all B-site cation displacements along single [111] direction. The different correlations among the <111> B-atom displacements can be determined from the diffuse scattering patterns. 1-dimensional correlation of the <111> B-atom displacements along <100> crystallographic directions gives rise to sheets of diffuse scattering intensity in the {100} reciprocal lattice planes, such as shown in Figure 4c. Based on the types of atomic displacements in the different average crystallographic phase, one can therefore expect to have three sheets of diffuse scattering in the cubic phase, two such sheets in the tetragonal phase, and only one such sheet in the orthorhombic phase; since the atomic displacements are correlated over long-distances in the rhombohedral phase, no such diffuse scattering sheets are visible.  Based on the Scherrer equation, the thickness of the diffuse scattering sheets is inversely proportional to the correlation lengths of the atomic displacements ordering. The dimensions of the diffuse scattering sheets can be determined from 2-dimensional scattering patterns using an area detector, such as shown in Figure 5, thereby providing a semi-quantitative means to determine the length scales of atomic correlations in a material. Using single-crystal diffuse scattering measurements with laboratory X-rays, the 1-dimensional correlations among B-atom displacements in BaTiO 3 and KNbO 3 were earlier estimated to be between 50 Å and 100 Å [54,55]. More recent measurements with high-energy X-rays provide estimates of the order of 2 nm [56,57], which is closer to that predicted from theoretical calculations [58]. In addition to 1-dimensional correlations, 2-dimensional correlations in the form of nanoscale polar clusters among B-atom off-center displacements are also possible [59,60].
It is further proposed that the addition of different atomic substitutions can modify short-range correlations within these polar clusters with correlated B-site atomic displacements [61][62][63]. While there have been some microscopy studies on the modifying effects of solid solution additives on short-range correlations in BaTiO 3 and KNbO 3 -based ferroelectric compositions, single-crystal diffuse scattering studies of these materials are however yet to be undertaken. There is therefore a lot of room to gain new knowledge about local atomic displacements ordering in these materials. Based on the Scherrer equation, the thickness of the diffuse scattering sheets is inversely proportional to the correlation lengths of the atomic displacements ordering. The dimensions of the diffuse scattering sheets can be determined from 2-dimensional scattering patterns using an area detector, such as shown in Figure 5, thereby providing a semi-quantitative means to determine the length scales of atomic correlations in a material. Using single-crystal diffuse scattering measurements with laboratory X-rays, the 1-dimensional correlations among B-atom displacements in BaTiO3 and KNbO3 were earlier estimated to be between 50 Å and 100 Å [54,55]. More recent measurements with high-energy X-rays provide estimates of the order of 2 nm [56,57], which is closer to that predicted from theoretical calculations [58]. In addition to 1-dimensional correlations, 2-dimensional correlations in the form of nanoscale polar clusters among B-atom off-center displacements are also possible [59,60].
It is further proposed that the addition of different atomic substitutions can modify short-range correlations within these polar clusters with correlated B-site atomic displacements [61][62][63]. While there have been some microscopy studies on the modifying effects of solid solution additives on short-range correlations in BaTiO3 and KNbO3-based ferroelectric compositions, single-crystal diffuse scattering studies of these materials are however yet to be undertaken. There is therefore a lot of room to gain new knowledge about local atomic displacements ordering in these materials. Most remarkably, in situ diffraction with synchrotron X-rays have shown that the nanoscale correlations for atomic displacements in both BaTiO3 and KNbO3 are highly susceptible to applied pressure and electric fields. The susceptibility of B-site atomic displacements on applied pressure was studied by Ravy et al. from in situ high-energy X-ray diffraction experiments [56]. They showed that the intensity of the diffuse scattering lines, as described above, are strong functions of the magnitude of applied pressures. Such behavior can be described as a result of changes in both magnitude and correlations among B-site atomic displacements. Similarly, the effect of applied electric fields on B-site atomic displacement correlations in single-crystal KNbO3 were shown by Pramanick et al. It was shown that application of a modest electric field of 300 V/mm in the room-temperature orthorhombic phase of KNbO3 along [001] can drastically change the length scale of Nb displacements from ∼2 nm to ∼7 nm [57]. This was determined from changes in the thickness of diffuse scattering lines as a function of applied electric fields, such as Figure 6. Most remarkably, in situ diffraction with synchrotron X-rays have shown that the nanoscale correlations for atomic displacements in both BaTiO 3 and KNbO 3 are highly susceptible to applied pressure and electric fields. The susceptibility of B-site atomic displacements on applied pressure was studied by Ravy et al. from in situ high-energy X-ray diffraction experiments [56]. They showed that the intensity of the diffuse scattering lines, as described above, are strong functions of the magnitude of applied pressures. Such behavior can be described as a result of changes in both magnitude and correlations among B-site atomic displacements. Similarly, the effect of applied electric fields on B-site atomic displacement correlations in single-crystal KNbO 3 were shown by Pramanick et al. It was shown that application of a modest electric field of 300 V/mm in the room-temperature orthorhombic phase of KNbO 3 along [001] can drastically change the length scale of Nb displacements from ∼2 nm to ∼7 nm [57]. This was determined from changes in the thickness of diffuse scattering lines as a function of applied electric fields, such as Figure 6.

Octahedral Tilt Disorder
In addition to off-center atomic displacements, tilting of oxygen octahedra constitutes another defining feature for crystal structure of perovskite ferroelectrics. The various relations between octahedral tilting arrangements and average crystal symmetry were originally explained by Glazer [64]. The perovskite crystal lattice can be viewed as a network of corner linked BO3 octahedra with interstitial A ions, such as illustrated in Figure 7a. The octahedra can be titled (or rotated) around the crystallographic a, b, or c axes, which are identified as placeholders in the Glazer notation. If the unit cell dimensions along two axes are the same, then the same letter can be used for two space holders. The rotation of the adjacent octahedra along a tilting axis can be in phase (marked by superscript + ) or antiphase (marked by superscript − ). For no tilting, a superscript of 0 is used. For example, the octahedral tilting patterns for a 0 a 0 c + and a − b − c − are shown in Figure 7b and c, respectively.

Octahedral Tilt Disorder
In addition to off-center atomic displacements, tilting of oxygen octahedra constitutes another defining feature for crystal structure of perovskite ferroelectrics. The various relations between octahedral tilting arrangements and average crystal symmetry were originally explained by Glazer [64]. The perovskite crystal lattice can be viewed as a network of corner linked BO 3 octahedra with interstitial A ions, such as illustrated in Figure 7a. The octahedra can be titled (or rotated) around the crystallographic a, b, or c axes, which are identified as placeholders in the Glazer notation. If the unit cell dimensions along two axes are the same, then the same letter can be used for two space holders. The rotation of the adjacent octahedra along a tilting axis can be in phase (marked by superscript + ) or antiphase (marked by superscript − ). For no tilting, a superscript of 0 is used. For example, the octahedral tilting patterns for a 0 a 0 c + and a − b − c − are shown in Figure 7b and c, respectively. The effect of oxygen octahedral tilting becomes particularly important for perovskite ferroelectrics with strong hybridization between s and p orbitals of the neighboring oxygen and A-site ions, such as in PZT. For Pb-free ferroelectrics, this is exemplified by the system Na0.5Bi0.5TiO3 (NBT) and its solid solutions. The changes in the average crystal structure in this system comprise both re-orientation of the polarization direction, as well as the octahedral tilt arrangements [65,66]. Below the TC ∼840 K, Na0.5Bi0.5TiO3 assumes a ferroelectric/ferrielectric phase with average polarization direction parallel to [001] and a 0 a 0 c + octahedral tilt arrangement. Upon cooling below T∼570 K, the system transitions to a ferroelectric phase with average rhombohedral symmetry, which is defined by polarization parallel to <111> and a − a − a − octahedral tilt arrangement. The lower temperature phase is also identified in some publications as having monoclinic symmetry, instead of rhombohedral, although that viewpoint is debated in other places [66][67][68][69][70][71]. In any case, it is now understood that Na0.5Bi0.5TiO3 and its solid solutions exhibit considerable nanoscale disorder and therefore assertions of any long-range crystal symmetry for this system is going to be open to interpretations. Single-crystal X-ray and neutron diffraction of Na0.5Bi0.5TiO3-based ferroelectrics has helped to clarify many of the subtleties in this system with regard to structural disorders.
The diffuse scattering pattern for the Na0.5Bi0.5TiO3 system shows several interesting features such as described below [45,[72][73][74][75][76][77]. Single-crystal diffraction revealed presence of diffuse scattering streaks, which are asymmetrically extended from the Bragg peaks along the lower-angular <001> directions. The peculiar extension of the streaks, as well as their intensities, are hkl-dependent. Based on detailed analyses of both X-ray and neutron diffuse scattering patterns from NBT single crystals, a general model could be proposed whereby the diffuse streaks are caused by the inclusion of nanoscale tetragonal plate-like structures within a rhombohedral matrix. Since the tetragonal and rhombohedral phases have different octahedral tilting arrangements, such an inclusion can be viewed as stacking faults separating adjacent rhombohedral regions. The reason for such tetragonal plate-like nanostructures is not entirely clear, although models based on atomic segregation similar to Guinier-Preston zones (GPZ) and memory effect leading to persistence of local regions of higher temperature phase have been proposed [45,76]. The asymmetric nature of the streaks has been tentatively explained as due to the slight expansion of the inclusions with respect to the matrix phase; nevertheless, the issue needs to be investigated further [76]. In addition to the asymmetric <001> diffuse scattering streaks, ridges of diffuse scattering were also observed with orientations along <110>, which appear to be correlated with the presence of random local fields similar to some Pb-based ferroelectrics [73]. Additionally, Huang scattering due to presence of long-range strain fields, supposedly caused by the presence of nanoscale inclusions, has also been noted [76]. All the The effect of oxygen octahedral tilting becomes particularly important for perovskite ferroelectrics with strong hybridization between s and p orbitals of the neighboring oxygen and A-site ions, such as in PZT. For Pb-free ferroelectrics, this is exemplified by the system Na 0.5 Bi 0.5 TiO 3 (NBT) and its solid solutions. The changes in the average crystal structure in this system comprise both re-orientation of the polarization direction, as well as the octahedral tilt arrangements [65,66]. Below the T C ∼840 K, Na 0.5 Bi 0.5 TiO 3 assumes a ferroelectric/ferrielectric phase with average polarization direction parallel to [001] and a 0 a 0 c + octahedral tilt arrangement. Upon cooling below T ∼570 K, the system transitions to a ferroelectric phase with average rhombohedral symmetry, which is defined by polarization parallel to <111> and a − a − a − octahedral tilt arrangement. The lower temperature phase is also identified in some publications as having monoclinic symmetry, instead of rhombohedral, although that viewpoint is debated in other places [66][67][68][69][70][71]. In any case, it is now understood that Na 0.5 Bi 0.5 TiO 3 and its solid solutions exhibit considerable nanoscale disorder and therefore assertions of any long-range crystal symmetry for this system is going to be open to interpretations. Single-crystal X-ray and neutron diffraction of Na 0.5 Bi 0.5 TiO 3 -based ferroelectrics has helped to clarify many of the subtleties in this system with regard to structural disorders.
The diffuse scattering pattern for the Na 0.5 Bi 0.5 TiO 3 system shows several interesting features such as described below [45,[72][73][74][75][76][77]. Single-crystal diffraction revealed presence of diffuse scattering streaks, which are asymmetrically extended from the Bragg peaks along the lower-angular <001> directions. The peculiar extension of the streaks, as well as their intensities, are hkl-dependent. Based on detailed analyses of both X-ray and neutron diffuse scattering patterns from NBT single crystals, a general model could be proposed whereby the diffuse streaks are caused by the inclusion of nanoscale tetragonal plate-like structures within a rhombohedral matrix. Since the tetragonal and rhombohedral phases have different octahedral tilting arrangements, such an inclusion can be viewed as stacking faults separating adjacent rhombohedral regions. The reason for such tetragonal plate-like nanostructures is not entirely clear, although models based on atomic segregation similar to Guinier-Preston zones (GPZ) and memory effect leading to persistence of local regions of higher temperature phase have been proposed [45,76]. The asymmetric nature of the streaks has been tentatively explained as due to the slight expansion of the inclusions with respect to the matrix phase; nevertheless, the issue needs to be investigated further [76]. In addition to the asymmetric <001> diffuse scattering streaks, ridges of diffuse scattering were also observed with orientations along <110>, which appear to be correlated with the presence of random local fields similar to some Pb-based ferroelectrics [73]. Additionally, Huang scattering due to presence of long-range strain fields, supposedly caused by the presence of nanoscale inclusions, has also been noted [76]. All the diffuse scattering features were observed to be dependent on composition for Na 0.5 Bi 0.5 TiO 3 -BaTiO 3 (NBT-BT) solid solutions, which indicated a gradual transition from a shorter-to-longer range ordering of polar atomic displacements and octahedral tilting [77]. These observations indicate that ferroelectrics based on Na 0.5 Bi 0.5 TiO 3 are complex nanostructured oxides with hierarchical structures, which require extensive studies for their comprehensive understanding.
Most interestingly, similar to what is described in Section 2.2.1. A, the nanoscale structures showed strong susceptibility for changes under the application of external stimuli. Kreisel et al. showed from in situ high-energy X-ray diffraction experiments that application of hydrostatic pressure up to 2.8 GPa first leads to decrease in the intensity of the asymmetric diffuse scattering peaks in Na 0.5 Bi 0.5 TiO 3 single crystals, and finally to their eventual disappearance (Figure 8a) under 3.7 GPa [45]. Similar decrease in intensities of L-shaped asymmetric diffuse scattering peaks could also be observed in Na 0.5 Bi 0.5 TiO 3 single crystals under applied electric fields (Figure 8b) [78,79]. Such observations indicated that the structural differences between the matrix and the nanoscales inclusions could be irreversibly overcome by application of mechanical pressure or electric fields, ultimately leading to a phase transition in the short-range order associated with nanoscale inclusions. diffuse scattering features were observed to be dependent on composition for Na0.5Bi0.5TiO3-BaTiO3 (NBT-BT) solid solutions, which indicated a gradual transition from a shorter-to-longer range ordering of polar atomic displacements and octahedral tilting [77]. These observations indicate that ferroelectrics based on Na0.5Bi0.5TiO3 are complex nanostructured oxides with hierarchical structures, which require extensive studies for their comprehensive understanding. Most interestingly, similar to what is described in Section II. A, the nanoscale structures showed strong susceptibility for changes under the application of external stimuli. Kreisel et al. showed from in situ high-energy X-ray diffraction experiments that application of hydrostatic pressure up to 2.8 GPa first leads to decrease in the intensity of the asymmetric diffuse scattering peaks in Na0.5Bi0.5TiO3 single crystals, and finally to their eventual disappearance (Figure 8a) under 3.7 GPa [45]. Similar decrease in intensities of L-shaped asymmetric diffuse scattering peaks could also be observed in Na0.5Bi0.5TiO3 single crystals under applied electric fields (Figure 8b) [78,79]. Such observations indicated that the structural differences between the matrix and the nanoscales inclusions could be irreversibly overcome by application of mechanical pressure or electric fields, ultimately leading to a phase transition in the short-range order associated with nanoscale inclusions. The above-mentioned changes in nanoscale structural correlations are significant since they have a direct effect on electromechanical susceptibility, as discussed in the following section. These initial studies on prototypical compounds have demonstrated that external field-or stress-induced changes in atomic displacements correlations should be studied in more detail in the solid solutions of these compounds. The availability of high quality single crystals will certainly be a facilitating factor for undertaking such studies.

Lattice Instability
The influence of nanoscale structural correlations towards macroscopic electromechanical properties of ferroelectrics can be better appreciated if fundamental relationships between the two are firmly established in addition to the already apparent empirical relationships. However, establishing clear physical relationships between characteristics of nanoscale structural domains and electromechanical susceptibility is not straightforward and only a few examples of investigative studies exists in this area, such as reference [46]. The critical linkage between short-range atomic The above-mentioned changes in nanoscale structural correlations are significant since they have a direct effect on electromechanical susceptibility, as discussed in the following section. These initial studies on prototypical compounds have demonstrated that external field-or stress-induced changes in atomic displacements correlations should be studied in more detail in the solid solutions of these compounds. The availability of high quality single crystals will certainly be a facilitating factor for undertaking such studies.

Lattice Instability
The influence of nanoscale structural correlations towards macroscopic electromechanical properties of ferroelectrics can be better appreciated if fundamental relationships between the two are firmly established in addition to the already apparent empirical relationships. However, establishing clear physical relationships between characteristics of nanoscale structural domains and electromechanical susceptibility is not straightforward and only a few examples of investigative studies exists in this area, such as reference [46]. The critical linkage between short-range atomic ordering within domains and macroscopic electromechanical susceptibility can be found through investigation of lattice vibration modes or phonons, which informs about the various inherent or evolving lattice instabilities. Inelastic neutron (and more recently, X-ray) scattering provides an excellent means for characterization of phonons in piezoelectric crystals. While strictly speaking this falls outside the scope of their current review, which is focused on diffraction or elastic scattering, a brief mention is made here of the characterization of phonons in Pb-free piezoelectric crystals in recognition of their perceived importance.
The central role of lattice dynamics towards the electromechanical properties of crystalline materials was originally recognized by Cochran [80]. From a basic viewpoint, it can be understood that the intrinsic susceptibility for mechanical deformation of a ferroelectric crystal is related to the elastic instability of the lattice [20], which is furthermore related to the propagation characteristics of long-wavelength acoustic phonon modes [81]. In gauging the effects of nanoscale ordered domains on electromechanical properties, one should therefore look at the correlations between structures and the propagation of acoustic phonon modes. Such investigations could indicate new pathways for designing new compositions of Pb-free ferroelectrics, as indicated from the following examples.
In prototypical KNbO 3 , it was demonstrated much earlier that the short-range correlations among the off-centered Nb displacements have a strong effect on the propagation of acoustic phonons; however, such effects could only be observed for phonons with shorter wavelengths or longer wavevectors (q > 0.2) [82]. In an interesting demonstration, using high-energy X-ray single-crystal diffraction and inelastic neutron scattering measurement, it was shown that when the correlation length for Nb-displacements is increased from ∼2 nm to ∼7 nm, the propagation of even longer wavelength (or shorter wavevector) acoustic phonons are damped in a transverse direction (see Figure 9) [57]. This indicated how controlling the dimensions of the nanoscale domains that are related to correlated atomic displacements could enhance the susceptibility for electromechanical shear deformation in this ferroelectric material. A similar effect of electric-field induced enhanced lattice instability was also demonstrated for BaTiO 3 single crystal [83]. ordering within domains and macroscopic electromechanical susceptibility can be found through investigation of lattice vibration modes or phonons, which informs about the various inherent or evolving lattice instabilities. Inelastic neutron (and more recently, X-ray) scattering provides an excellent means for characterization of phonons in piezoelectric crystals. While strictly speaking this falls outside the scope of their current review, which is focused on diffraction or elastic scattering, a brief mention is made here of the characterization of phonons in Pb-free piezoelectric crystals in recognition of their perceived importance. The central role of lattice dynamics towards the electromechanical properties of crystalline materials was originally recognized by Cochran [80]. From a basic viewpoint, it can be understood that the intrinsic susceptibility for mechanical deformation of a ferroelectric crystal is related to the elastic instability of the lattice [20], which is furthermore related to the propagation characteristics of long-wavelength acoustic phonon modes [81]. In gauging the effects of nanoscale ordered domains on electromechanical properties, one should therefore look at the correlations between structures and the propagation of acoustic phonon modes. Such investigations could indicate new pathways for designing new compositions of Pb-free ferroelectrics, as indicated from the following examples.
In prototypical KNbO3, it was demonstrated much earlier that the short-range correlations among the off-centered Nb displacements have a strong effect on the propagation of acoustic phonons; however, such effects could only be observed for phonons with shorter wavelengths or longer wavevectors ( [82]. In an interesting demonstration, using high-energy X-ray single-crystal diffraction and inelastic neutron scattering measurement, it was shown that when the correlation length for Nb-displacements is increased from ∼2 nm to ∼7 nm, the propagation of even longer wavelength (or shorter wavevector) acoustic phonons are damped in a transverse direction (see Figure 9) [57]. This indicated how controlling the dimensions of the nanoscale domains that are related to correlated atomic displacements could enhance the susceptibility for electromechanical shear deformation in this ferroelectric material. A similar effect of electric-field induced enhanced lattice instability was also demonstrated for BaTiO3 single crystal [83]. The dampening of the phonons, observed from a broadening of the inelastic neutron scattering peak, indicate increased lattice instability and hence greater susceptibility to mechanical deformation under field. After Pramanick et al. [57] Reprinted from reference [57], with permission from John Wiley and Sons.
In the case of Na0.5Bi0.5TiO3, a recent inelastic neutron scattering study by Luo et al. indicated dynamic tilting modes of the oxygen octahedra, which are possibly coupled to the dimensions of the nanoscale domain structures [77]. Detailed in situ studies might show how short-range ordering of octahedral tilts and their dynamic vibrational modes evolve under applied electric fields or  In the case of Na 0.5 Bi 0.5 TiO 3 , a recent inelastic neutron scattering study by Luo et al. indicated dynamic tilting modes of the oxygen octahedra, which are possibly coupled to the dimensions of the nanoscale domain structures [77]. Detailed in situ studies might show how short-range ordering of octahedral tilts and their dynamic vibrational modes evolve under applied electric fields or mechanical pressure, and therefore better elucidate the role of nanoscale domains on macroscopic electromechanical properties.

Domain Reorientation
In addition to the structural changes described above, re-orientation of microscopic domains can also significantly contribute to the macroscopic electromechanical response of a ferroelectric crystal [35,43]. It is therefore important to be able to characterize the nature and extent of domain reorientation in ferroelectric crystals under the application of an external stimulus. Single-crystal diffraction can provide clear and detailed insights into domain reorientation pathways, which are not readily obtainable from experiments on polycrystalline ceramics.
The basis for characterization of domain reorientation from single-crystal diffraction is explained below. In order to maintain continuity across a domain boundary, the two adjacent ferroelectric domains in a ferroelectric crystal are misoriented by a certain angular range, which is defined by crystallographic considerations. For example, in tetragonal crystals, the misorientation angle is θ ∼ 2 tan −1 (c/a) between the ferroelectric domains separated by non-180 • domain walls [84], where c and a are the crystallographic unit cell parameters, as schematically illustrated in Figure 10. This effect causes splitting of the diffraction peaks coming from the different domain variants in the crystal. Based on the relative intensities of the symmetry-related diffraction spots, the volume fraction for the different domain variants can be estimated. This is illustrated below for the case of domain switching phenomena in Na 0.5 Bi 0.5 TiO 3 single crystals. mechanical pressure, and therefore better elucidate the role of nanoscale domains on macroscopic electromechanical properties.

Domain Reorientation
In addition to the structural changes described above, re-orientation of microscopic domains can also significantly contribute to the macroscopic electromechanical response of a ferroelectric crystal [35,43]. It is therefore important to be able to characterize the nature and extent of domain reorientation in ferroelectric crystals under the application of an external stimulus. Single-crystal diffraction can provide clear and detailed insights into domain reorientation pathways, which are not readily obtainable from experiments on polycrystalline ceramics.
The basis for characterization of domain reorientation from single-crystal diffraction is explained below. In order to maintain continuity across a domain boundary, the two adjacent ferroelectric domains in a ferroelectric crystal are misoriented by a certain angular range, which is defined by crystallographic considerations. For example, in tetragonal crystals, the misorientation angle is ( ) a c / tan 2 1 − θ between the ferroelectric domains separated by non-180° domain walls [84], where c and a are the crystallographic unit cell parameters, as schematically illustrated in Figure  10. This effect causes splitting of the diffraction peaks coming from the different domain variants in the crystal. Based on the relative intensities of the symmetry-related diffraction spots, the volume fraction for the different domain variants can be estimated. This is illustrated below for the case of domain switching phenomena in Na0.5Bi0.5TiO3 single crystals. Using high-energy synchrotron X-rays, Kitanaka et al. studied the microscopic mechanism for polarization switching in Na0.5Bi0.5TiO3 single crystals under electric fields applied parallel to <100> and <111> directions [85]. Based on the differences in intensities of the symmetry-related diffraction spots, such as shown in Figure 11, it was demonstrated that polarization switching in both cases was dominated by step-wise non-180° (71° and/or 109°) domain switching processes. A preference for 71° reorientation of the polarization direction was furthermore asserted based on theoretical calculations using DFT.
(a) Using high-energy synchrotron X-rays, Kitanaka et al. studied the microscopic mechanism for polarization switching in Na 0.5 Bi 0.5 TiO 3 single crystals under electric fields applied parallel to <100> and <111> directions [85]. Based on the differences in intensities of the symmetry-related diffraction spots, such as shown in Figure 11, it was demonstrated that polarization switching in both cases was dominated by step-wise non-180 • (71 • and/or 109 • ) domain switching processes. A preference for 71 • reorientation of the polarization direction was furthermore asserted based on theoretical calculations using DFT. polarization switching in Na0.5Bi0.5TiO3 single crystals under electric fields applied parallel to <100> and <111> directions [85]. Based on the differences in intensities of the symmetry-related diffraction spots, such as shown in Figure 11, it was demonstrated that polarization switching in both cases was dominated by step-wise non-180° (71° and/or 109°) domain switching processes. A preference for 71° reorientation of the polarization direction was furthermore asserted based on theoretical calculations using DFT. A similar investigation of polarization switching behavior in K0.25Na0.75NbO3 single-crystals was undertaken by Deng et al. [86]. The macroscopic symmetry for this crystal was determined to be of monoclinic C type, with the spontaneous polarization direction close to <054>. Based on both single-crystal X-ray diffraction and polarized light microscopy, it was determined that polarization switching is again dominated by non-180° domain reorientations processes. This was evident from change in relative intensities of the symmetry-related diffraction intensities of (101)pc reflections from a (101) oriented crystal for different amplitudes of applied electric fields.

Opportunities and Challenges
Above, I have reviewed some recent developments with regard to the characterization of Pb-free piezoelectric materials using single-crystal X-ray and neutron diffraction. These studies have broadened our understanding of structure-property relationships in this new class of materials and further opened up new avenues of investigation. In fact, they have shown that the science of correlating complex materials chemistry with local bonding environments, multiscale structural disorders, and ultimately, electromechanical properties under applied stimuli for the new Pb-free piezoelectric materials is still in its infancy. In this regard, there exists huge opportunities to gain new knowledge in this area. As reviewed above, single-crystal diffraction can be a very effective characterization technique for this purpose, especially under in situ experimental conditions. It needs to be mentioned here that care must be taken while analyzing single-crystal diffraction data for characterization of multi-scale structural changes. If measurements are specifically focused near single ZC, it might lead to incorrect conclusions. For example, additional spots can appear near a ZC due to either change in symmetry due to a crystallographic phase transition (or distortion) or non-180° domain reorientation. A comprehensive survey of a broad region of reciprocal space is therefore necessary. This is facilitated by two significant developments in terms of instrumental capabilities. With the advent of 3rd generation synchrotron and high-energy X-rays, large areas of reciprocal space can be surveyed in minutes. This is due to the fact that for small X-ray wavelengths (∼0.1 Å), the Ewald sphere is essentially flat, and therefore reflections within a single hkl plane can be obtained in a single shot [52]. Secondly, neutron Laue diffraction is gaining in popularity with the advent of spallation sources and time-of-flight Laue diffractometers. Modern instruments, such as TOPAZ at the Spallation Neutron Source (SNS), have extensive area-detector coverage that are A similar investigation of polarization switching behavior in K 0.25 Na 0.75 NbO 3 single-crystals was undertaken by Deng et al. [86]. The macroscopic symmetry for this crystal was determined to be of monoclinic C type, with the spontaneous polarization direction close to <054>. Based on both single-crystal X-ray diffraction and polarized light microscopy, it was determined that polarization switching is again dominated by non-180 • domain reorientations processes. This was evident from change in relative intensities of the symmetry-related diffraction intensities of (101) pc reflections from a (101) oriented crystal for different amplitudes of applied electric fields.

Opportunities and Challenges
Above, I have reviewed some recent developments with regard to the characterization of Pb-free piezoelectric materials using single-crystal X-ray and neutron diffraction. These studies have broadened our understanding of structure-property relationships in this new class of materials and further opened up new avenues of investigation. In fact, they have shown that the science of correlating complex materials chemistry with local bonding environments, multiscale structural disorders, and ultimately, electromechanical properties under applied stimuli for the new Pb-free piezoelectric materials is still in its infancy. In this regard, there exists huge opportunities to gain new knowledge in this area. As reviewed above, single-crystal diffraction can be a very effective characterization technique for this purpose, especially under in situ experimental conditions. It needs to be mentioned here that care must be taken while analyzing single-crystal diffraction data for characterization of multi-scale structural changes. If measurements are specifically focused near single ZC, it might lead to incorrect conclusions. For example, additional spots can appear near a ZC due to either change in symmetry due to a crystallographic phase transition (or distortion) or non-180 • domain reorientation. A comprehensive survey of a broad region of reciprocal space is therefore necessary. This is facilitated by two significant developments in terms of instrumental capabilities. With the advent of 3rd generation synchrotron and high-energy X-rays, large areas of reciprocal space can be surveyed in minutes. This is due to the fact that for small X-ray wavelengths (∼0.1 Å), the Ewald sphere is essentially flat, and therefore reflections within a single hkl plane can be obtained in a single shot [52]. Secondly, neutron Laue diffraction is gaining in popularity with the advent of spallation sources and time-of-flight Laue diffractometers. Modern instruments, such as TOPAZ at the Spallation Neutron Source (SNS), have extensive area-detector coverage that are capable of three-dimensional → Q-space mapping from a stationary single crystal [87]. In a time-of-flight neutron diffractometer, the scattering intensities from the crystal appear at the detector at different times of flight, which correspond to different wavelengths, and therefore to different wave vectors → Q. Consequently, in addition to resolving the scattered intensities in two dimensions with the use of area detectors, scattering intensities are also recorded by the arrival times of the neutrons at the detector. In this manner, the three-dimensional scattering patterns can be obtained very effectively without the need for any rotation of the crystal during data collection. These advancements in instrumental capabilities also open the possibility to conduct time-resolved experiments, which will enable characterization of transient states during application of external stimuli. In addition, care must be taken to distinguish between elastic diffuse scattering, which informs about the static local disorder in a material and the inelastic scattering which informs about the dynamic structural states. This is often not possible for X-ray instruments, where the scattering intensities are integrated over all possible energy-transfers. Modern neutron scattering instruments at spallation sources, such as CORELLI at SNS [88], affords the possibility for such a distinction between elastic and inelastic diffuse scattering. Finally, one must also be careful to account for intensity changes due to changes in internal strains and extinction effects, while analyzing diffraction data from single-crystal samples.
It goes without saying that the availability of good-quality single-crystals is a must for undertaking comprehensive structural measurements. This means that the sample should be of adequate dimension for a specific experimental purpose. Measurement of inelastic signals and diffuse scattering often needs larger samples as compared to measurement of Bragg diffraction peaks. Additionally, the samples should be properly annealed to reduce internal stress related effects. At the very least, the sample preparation procedures must be clearly explained. Additionally, if the experimental goal is to obtain diffraction data for precise structural refinement, then it is better that the single-crystal sample contains a single domain or twin variant. However, this requirement could be relaxed by comprehensive data collection over large area of the reciprocal space and properly indexing the peaks from different domain variants using advanced software [87].
With continuing improvements in the areas of single-crystal growth techniques and capabilities of scattering instruments, it should be possible to significantly advance our understanding of the fundamental physico-chemical principles that drive large enhancements in Pb-free piezoelectrics, and therefore, to accelerate their discovery process.