Enhanced Piezoelectricity and Thermal Stability of Electrostrain Performance in BiFeO3-Based Lead-Free Ceramics

BiFeO3–based ceramics possess an advantage over large spontaneous polarization and high Curie temperature, and are thus widely explored in the field of high–temperature lead–free piezoelectrics and actuators. However, poor piezoelectricity/resistivity and thermal stability of electrostrain make them less competitive. To address this problem, (1 − x) (0.65BiFeO3–0.35BaTiO3)–xLa0.5Na0.5TiO3 (BF–BT–xLNT) systems are designed in this work. It is found that piezoelectricity is significantly improved with LNT addition, which is contributed by the phase boundary effect of rhombohedral and pseudocubic phase coexistence. The small–signal and large–signal piezoelectric coefficient (d33 and d33*) peaks at x = 0.02 with 97 pC/N and 303 pm/V, respectively. The relaxor property and resistivity are enhanced as well. This is verified by Rietveld refinement, dielectric/impedance spectroscopy and piezoelectric force microscopy (PFM) technique. Interestingly, a good thermal stability of electrostrain is obtained at x = 0.04 composition with fluctuation η = 31% (Smax'−SRTSRT×100%), in a wide temperature range of 25–180 °C, which is considered as a compromise of negative temperature dependent electrostrain for relaxors and the positive one for ferroelectric matrix. This work provides an implication for designing high–temperature piezoelectrics and stable electrostrain materials.


Introduction
As a mechanical-electricity conversion functional material, piezoelectric ceramics have wide applications in the defense, industrial, and medical fields, etc. [1][2][3]. Due to extraordinary piezoelectric properties and electromechanical coupling effects near morphotropic phase boundary (MPB), Pb(Zr, Ti)O 3 -based (PZT) ceramics have been a hot research topic [4,5]. However, with an increasingly serious environmental concern, lead-based materials are gradually replaced by lead-free materials. A series of lead-free materials such as (Bi 0.5 Na 0.5 )TiO 3 (BNT-) [6][7][8], BiFeO 3 (BF-) [9][10][11], (K 0.5 Na 0.5 )NbO 3 (KNN-) based [12,13] ceramics have been widely developed. For example, Liu et al. reported a non-textured BNT-based ceramic with a large electrostrain value of~0.7% at room temperature (RT), but it exhibited relatively large strain hysteresis and a decrease in electrostrain as temperature increased [14]; Due to polymorphic phase transition effect in KNN-based ceramics, the macroscopic performance usually exhibited a temperature sensitivity, although compositionally graded multilayer composite and layered distribution of dopants strategy largely

Materials and Methods
We obtained (1 − x) (0.65BiFeO 3 -0.35BaTiO 3 )-xLa 0.5 Na 0.5 TiO 3 (BF-BT-xLNT, x = 0~0.06) ceramics by conventional solid-state synthesis using Bi 2 O 3 , Fe 2 O 3 , TiO 2 , BaCO 3 , La 2 O 3, and Na 2 CO 3 . Excessive 2%mol Bi 2 O 3 was added for compensation. All raw materials were mixed with alcohol and ball-milled for 15 h and then calcined at 750 • C for 6 h twice. Before the second ball-milling step, 1% mol MnO 2 was added to reduce the leakage current. Then, powders were mixed with 8 wt.% PVA binder and pressed into the discs with a diameter of 10 mm. After removing the binder, the discs were sintered at 1040 • C for 3 h. For electrical measurements, a silver paste was coated on both sides of the polished samples and fired at 560 • C for 10 min.
The crystal structure was measured by X-ray diffractometer (XRD, Rigaku Smart-lab). A scanning electron microscope (SEM, Regulus 8230; Hitachi Co., Tokyo, Japan) was used to observe the sample microstructures. Before observing the microstructure, all sample surfaces are polished smooth and hot corroded at 950 • C for 30 min. A Raman microscope was used to acquire Raman spectra (Horiba Jobin-Yvon HR800, France). Dielectric/impedance properties were acquired by Wayne Kerr 6500B impedance analyzer (Wayne Kerr Electronic Instrument Co., Shenzhen, China). Polarization hysteresis (P-E) loops and electrostrain (S-E) curves were measured using a ferroelectric measuring system (Precision LC, Radiant Technologies, Inc. Albuquerque, NM, USA) at a frequency of 1 Hz. The sample was poled for 15 min under an electric field (E = 60 kV/cm) at 120 • C, and then piezoelectric coefficients d 33 were measured with a quasi-static d 33 m (YE2730A, China). The domain structure was characterized by the PFM technique (Asylum Research).

Results and Discussion
Room temperature (RT) XRD patterns of BF-BT-xLNT ceramics (x = 0~0.06) are shown in Figure 1a-g, respectively. All samples have a pure perovskite structure without a second phase, indicating that LNT is completely dissolved into BF-BT matrix and forms a solid solution. To get an in-depth understanding of phase structure and content evolution with changing LNT, Rietveld refinement is performed and R3cH (R phase) and Pm3m (Pc phase) space group models are exploited [30,31]. Low fitted values of R wp , R p, and χ 2 indicate that fitting results are reliable. Also, locally magnified (111) diffraction peaks are displayed in their insets. Obviously, a wide (111) peak for the x = 0 sample indicates the existence of the R phase. The (111) peak gradually becomes narrow and sharp with an increase in LNT content, indicating that the R phase is gradually substituted by a Pc phase, and finally evolves into a single Pc phase at x = 0.05, as displayed in Figure 1h. The refined lattice parameter is plotted in Figure 1i. Table 1 also shows the lattice parameters and R-factors obtained by Rietveld refinement for better understanding. Obviously, with LNT addition, the lattice parameter generally exhibits a downward shift, which is mainly attributed to smaller ionic radii of Na + and La 3+

Results and Discussion
Room temperature (RT) XRD patterns of BF-BT-xLNT ceramics (x = 0~0.06) are shown in Figure 1a-g, respectively. All samples have a pure perovskite structure without a second phase, indicating that LNT is completely dissolved into BF-BT matrix and forms a solid solution. To get an in-depth understanding of phase structure and content evolution with changing LNT, Rietveld refinement is performed and R3cH (R phase) and Pm3 ̅ m (Pc phase) space group models are exploited [30,31]. Low fitted values of Rwp, Rp, and χ 2 indicate that fitting results are reliable. Also, locally magnified (111) diffraction peaks are displayed in their insets. Obviously, a wide (111) peak for the x = 0 sample indicates the existence of the R phase. The (111) peak gradually becomes narrow and sharp with an increase in LNT content, indicating that the R phase is gradually substituted by a Pc phase, and finally evolves into a single Pc phase at x = 0.05, as displayed in Figure 1h. The refined lattice parameter is plotted in Figure 1i. Table 1 also shows the lattice parameters and R-factors obtained by Rietveld refinement for better understanding. Obviously, with LNT addition, the lattice parameter generally exhibits a downward shift, which is mainly attributed to smaller ionic radii of Na + and La 3+ (CN = 12, RNa + = 1.39 Å , RLa 3+ = 1.36 Å ) than that of Ba 2+ and Bi 3+ (CN = 12, RBa 2+ = 1.61 Å , RBi 3+ = 1.45 Å ) [32].       Raman spectra technique is a powerful tool to detect phase transition at short-range scales, and Raman spectra for BF-BT-xLNT ceramics are performed, as shown in Figure  3a. Generally, the wavenumber range in 50~1000 cm −1 is divided into three vibrational modes, namely A-site (50~200 cm −1 ), B-O bond (200~400 cm −1 ), and BO6 octahedral vibration (400~1000 cm −1 ) [33,34]. To clearly demonstrate vibrational mode changes, raw Raman spectra are fitted by the Lorentzian function, and a series of deconvoluted Raman peaks are delineated, as shown in Figure 3b. Here, two representative peaks (G and Raman spectra technique is a powerful tool to detect phase transition at short-range scales, and Raman spectra for BF-BT-xLNT ceramics are performed, as shown in Figure 3a. Generally, the wavenumber range in 50~1000 cm −1 is divided into three vibrational modes, namely A-site (50~200 cm −1 ), B-O bond (200~400 cm −1 ), and BO 6 octahedral vibration (400~1000 cm −1 ) [33,34]. To clearly demonstrate vibrational mode changes, raw Raman spectra are fitted by the Lorentzian function, and a series of deconvoluted Raman peaks are delineated, as shown in Figure 3b. Here, two representative peaks (G and H bands) are Nanomaterials 2023, 13, 942 5 of 14 selected to analyze the wavenumber and FWHM (full width at half maximum) evolution to detect phase transition, as shown in Figure 3c,d. Notably, two discontinuous changes in wavenumber and FWHM are observed (as highlighted by shadow), which strongly suggests that phase transition occurs [35]. Since the Raman shift is interrelated to crystal stress and polarization, the first abrupt change of G and H bands probably corresponds to distorted local stress and the polarization field [36]. With an increase of heterovalent ionic proportion (Na + and La 3+ ) and corresponding local random field, the second vibration change may be related to the ferroelectric-to-relaxor (FR) phase transition, which will be discussed infra in detail. H bands) are selected to analyze the wavenumber and FWHM (full width at half maximum) evolution to detect phase transition, as shown in Figure 3c,d. Notably, two discontinuous changes in wavenumber and FWHM are observed (as highlighted by shadow), which strongly suggests that phase transition occurs [35]. Since the Raman shift is interrelated to crystal stress and polarization, the first abrupt change of G and H bands probably corresponds to distorted local stress and the polarization field [36]. With an increase of heterovalent ionic proportion (Na + and La 3+ ) and corresponding local random field, the second vibration change may be related to the ferroelectric-to-relaxor (FR) phase transition, which will be discussed infra in detail. Temperature dependence of dielectric constant (εr) of BF-BT-xLNT ceramics (x = 0~0.06) at 1 kHz~500 kHz are shown in Figure 4a-g, respectively. For x = 0 compositions, it shows a relaxor-like behavior near 300 °C, which is related to a dipolar relaxation caused by •• hopping [37]. The dielectric drift at high temperatures is caused by a large conductivity due to an increased •• motion, as highlighted in Figure 4a [38]. As x increases to 0.01 and 0.02, the dielectric peak becomes sharp, and εr at Tm (the temperature for dielectric maxima) also increases ( Figure 4h). This is also observed from an obvious hump at the imaginary part of the dielectric constant (ε′′) curve, as shown in the insets of Figure 4b,c. Notably, the high-temperature dielectric drift is significantly suppressed, indicating the resistivity is markedly improved with LNT addition, which is also observed from the dielectric loss (tanδ) in Figure 4i. As to x = 0.03~0.06, diffusive phase transition and frequency dispersion are clearly observed, exhibiting a relaxation property [39,40]. The dielectric curves for x = 0~0.06 samples at 100 kHz are collected in Figure  4h. The peak εr value peaks at x = 0.01 and 0.02 and then sharply decreases for further increasing LNT content, accompanied by a wide and diffuse dielectric shape. This is more clearly indicated by εr/εm versus T/Tm curves in Figure 4j. Therefore, an addition of an LNT component first improves the dielectric performance with an appropriate proportion of R and Pc phase. With an excessive LNT content, εr is significantly suppressed with a dominant relaxor Pc phase. In order to quantitatively describe relaxor properties Temperature dependence of dielectric constant (ε r ) of BF-BT-xLNT ceramics (x = 0~0.06) at 1 kHz~500 kHz are shown in Figure 4a-g, respectively. For x = 0 compositions, it shows a relaxor-like behavior near 300 • C, which is related to a dipolar relaxation caused by V .. o hopping [37]. The dielectric drift at high temperatures is caused by a large conductivity due to an increased V .. o motion, as highlighted in Figure 4a [38]. As x increases to 0.01 and 0.02, the dielectric peak becomes sharp, and ε r at T m (the temperature for dielectric maxima) also increases (Figure 4h). This is also observed from an obvious hump at the imaginary part of the dielectric constant (ε ) curve, as shown in the insets of Figure 4b,c. Notably, the high-temperature dielectric drift is significantly suppressed, indicating the resistivity is markedly improved with LNT addition, which is also observed from the dielectric loss (tanδ) in Figure 4i. As to x = 0.03~0.06, diffusive phase transition and frequency dispersion are clearly observed, exhibiting a relaxation property [39,40]. The dielectric curves for x = 0~0.06 samples at 100 kHz are collected in Figure 4h. The peak ε r value peaks at x = 0.01 and 0.02 and then sharply decreases for further increasing LNT content, accompanied by a wide and diffuse dielectric shape. This is more clearly indicated by ε r /ε m versus T/T m curves in Figure 4j. Therefore, an addition of an LNT component first improves the dielectric performance with an appropriate proportion of R and Pc phase. With an excessive Nanomaterials 2023, 13, 942 6 of 14 LNT content, ε r is significantly suppressed with a dominant relaxor Pc phase. In order to quantitatively describe relaxor properties for this system, two parameters of ∆T relax [T m (100 kHz) − T m (1 kHz)] and ∆T span (temperature span corresponds to ε r /ε m = 0.8 in Figure 4j) are adopted [41]. Obviously, ∆T relax and ∆T span exhibit an increasing trend as LNT content increases, as plotted in Figure 4k. Also, a modified Curie-Weiss law is used to denote the relaxor degree γ: , where C is Curie-Weiss constant, ε m is dielectric maxima and the parameter γ is adopted to reveal the relaxor degree [42,43]. As shown in Figure 4l, γ value steadily increases from 1.73 for x = 0 to 1.91 for x = 0.06 composition. Therefore, an enhanced relaxor properties are obtained with LNT addition.  Figure 4j) are adopted [41]. Obviously, ΔTrelax and ΔTspan exhibit an increasing trend as LNT content increases, as plotted in Figure 4k. Also, a modified Curie-Weiss law is used to denote the relaxor degree γ: where C is Curie-Weiss constant, εm is dielectric maxima and the parameter γ is adopted to reveal the relaxor degree [42,43]. As shown in Figure 4l, γ value steadily increases from 1.73 for x = 0 to 1.91 for x = 0.06 composition. Therefore, an enhanced relaxor properties are obtained with LNT addition. To determine the resistivity and conduction mechanism for BF-BT-xLNT ceramics, complex impedance spectra (CIS) are measured over 280-380 °C with an interval of 20 °C, as shown in Figure 5a-g. All ceramics exhibit a single semicircle complex impedance shape, which is related to a single relaxation mechanism of bulk response [44]. The data are well fitted by two parallel R-CPE equivalent circuits, as shown in the inset of Figure  5e, which represent the grain boundary and bulk (grain) contributions, respectively. As To determine the resistivity and conduction mechanism for BF-BT-xLNT ceramics, complex impedance spectra (CIS) are measured over 280-380 • C with an interval of 20 • C, as shown in Figure 5a-g. All ceramics exhibit a single semicircle complex impedance shape, which is related to a single relaxation mechanism of bulk response [44]. The data are well fitted by two parallel R-CPE equivalent circuits, as shown in the inset of Figure 5e, which represent the grain boundary and bulk (grain) contributions, respectively. As known, bulk resistance (R b ) often dominates ferroelectric ceramics, and an extrapolated intercept on Z axis corresponds to R b in CIS. Obviously, R b gradually decreases with increasing temperature, indicating a negative temperature coefficient behavior. The R b value is obtained by fitting the experimental CIS, as shown in Figure 5a-g. For a better understanding, the R b values are also listed in Table 2. Interestingly, the LNT addition greatly increases the R b of the BF-BT matrix (Figure 5h,i), which is consistent with a depressed tanδ. In addition, bulk conductivity (σ b ) and activation energy (E a ) are calculated by the following formulas [45]: where h is sample thickness, S is electrode area, k B is Boltzmann constant, and T is the measured temperature. Figure 5i depicts the σ b as a function of inverse of temperature, and the fitted value of E a is obtained between 0.93 and 1.24 eV, which is close to the E a of V .. o (~1 eV) [46]. Therefore, V .. o dominates the high-temperature conductivity and leads to a large leakage current, which is mainly due to the volatilization of bismuth and the reduction of Fe 3+ [47].
Polarization hysteresis (P-E) and corresponding current density (J-E) loops of BF-BT-xLNT ceramics at RT are measured at 1 Hz, as shown in Figure 6a,b. Low saturated/remanent polarization (P m /P r ) with an incomplete domain switching is observed for x = 0 composition. Interestingly, for x = 0.01 and 0.02 samples, saturated P-E loops and a sharp J-E peak are observed, featuring a normal ferroelectric. A sharp J-E peak usually indicates strong ferroelectricity with a fast domain switching behavior under E [26]. As to x = 0.03 and 0.04, splitting J-E peaks are observed, indicating an emergence of an intermediate state with nonergodic/ergodic and ferroelectric phase coexistence. Meanwhile, P m , P r, and coercive field (E c ) decrease strikingly, as plotted in Figure 6c. Slant P-E loops and J-E platform with low P r and E c value indicates a pure relaxor phase for x = 0.05 and 0.06 samples, which also accords with pure Pc phase for both compositions. In addition, bipolar and unipolar electrostrain (S-E) curves of BF-BT-xLNT ceramics are shown in Figure 6d,e, respectively. For x ≤ 0.03 samples, typical butterfly-shape bipolar S-E curves are present, and they gradually evolve into a sprout-shape [decrease in negative strain (S neg )] with increasing LNT content. In general, larger S neg indicates more non-180 • domain switching under E, and this predicts improvement of piezoelectric performance. For x = 0.04~0.06 samples, an almost zero-S neg indicates a gradual FR phase transition. Clearly, positive strain (S pos ) and S neg peaks x = 0.02 sample, and the normalized d * 33 (S max /E max ) is calculated as well, as plotted in Figure 6f. The d 33 and d * 33 value simultaneously achieve optimal value (97 pC/N and 303 pm/V) at x = 0.02, which is contributed by phase boundary effect with proper proportion of R and Pc phase content. Excessive addition of LNT degrades piezoelectricity and ferroelectricity and thus d * 33 decreases. To measure the temperature stability of piezoelectric performance for these compositions, d 33 is collected at elevated annealing temperature, as shown in Figure 7a, which is considered a crucial benchmark for evaluating practical high-temperature performances. Apparently, x = 0.02 sample not only possesses a peak d 33 value but also maintains a wide temperature span (≤240 • C). Also, frequency-dependent ε r for poled x = 0.02 sample is also measured after different annealing temperatures, as shown in Figure 7b. Typical resonance and anti-resonance peaks are observed for sufficiently poled samples. As the temperature increases, the peak gradually degrades, indicating depolarization steadily occurs [48]. These peaks dampen gradually and finally vanish as the temperature increased to~240 • C. This is also indicated by the contour plot, as indicated in Figure 7c. To investigate the depolarization mechanism on the local structure of representative x = 0.02 and 0.03 ceramics, in situ Raman spectra is performed, as shown in Figure 7d-g. All Raman spectra are , where I m (ω) is Raman intensity and n(ω, T) = 1/ exp ω kT − 1 is Bose-Einstein temperature factor [49]. The evolution of wavenumber and FWHM of selected D and G modes are plotted in Figure 7e,f,h,i. The corresponding discontinuities are highlighted by the shadow. The discontinuity occurs around T d~2 40 • C and 140 • C for x = 0.02 and 0.03 samples, respectively, which indicates an FR phase transition and also accords with the T d [50].
pressed tanδ. In addition, bulk conductivity (σb) and activation energy (Ea) are calculated by the following formulas [45]: where h is sample thickness, S is electrode area, kB is Boltzmann constant, and T is the measured temperature. Figure 5i depicts the σb as a function of inverse of temperature, and the fitted value of Ea is obtained between 0.93 and 1.24 eV, which is close to the Ea of •• (~1 eV) [46]. Therefore, •• dominates the high-temperature conductivity and leads to a large leakage current, which is mainly due to the volatilization of bismuth and the reduction of Fe 3+ [47].   formance. For x = 0.04~0.06 samples, an almost zero-Sneg indicates a gradual FR phase transition. Clearly, positive strain (Spos) and Sneg peaks x = 0.02 sample, and the normalized * ( / ) is calculated as well, as plotted in Figure 6f. The d33 and * value simultaneously achieve optimal value (97 pC/N and 303 pm/V) at x = 0.02, which is contributed by phase boundary effect with proper proportion of R and Pc phase content. Excessive addition of LNT degrades piezoelectricity and ferroelectricity and thus * decreases. To measure the temperature stability of piezoelectric performance for these compositions, d33 is collected at elevated annealing temperature, as shown in Figure 7a, which is considered a crucial benchmark for evaluating practical high-temperature performances. Apparently, x = 0.02 sample not only possesses a peak d33 value but also maintains a wide temperature span (≤240 °C). Also, frequency-dependent εr for poled x = 0.02 sample is also measured after different annealing temperatures, as shown in Figure 7b. Typical resonance and anti-resonance peaks are observed for sufficiently poled samples. As the temperature increases, the peak gradually degrades, indicating depolarization steadily occurs [48]. These peaks dampen gradually and finally vanish as the temperature increased to~240 °C. This is also indicated by the contour plot, as indicated in Figure  7c. To investigate the depolarization mechanism on the local structure of representative x = 0.02 and 0.03 ceramics, in situ Raman spectra is performed, as shown in Figure 7d  Temperature dependent unipolar S-E curves of BF-BT-xLNT ceramics are shown in Figure 8a-g. For x ≤ 0.02 compositions with normal ferroelectric character, they exhibit a positive temperature dependence of electrostrain, which mainly originates from thermally activated domain switching [30,51]. For x = 0.03 composition, the electrostrain shows a positive temperature dependence below T d~1 40 • C, and finally the electrostrain performance remains stable. For x = 0.04 sample with less concentration of ferroelectric domains (R/Pc = 0.25/0.75), the electrostrain shows better temperature stability, which is probably due to a compromise of thermally activated domain switching and agitation. As LNT content further increases, the temperature stability of electrostrain is maintained but strain value decreases. Here, the change rate of electrostrain with temperature is defined as η = S max −S RT S RT × 100%, where S max is unipolar electrostrain peak value in the test temperature range of 25-180 • C, S RT is the electrostrain at RT. Evolution of unipolar S-E with temperature for x = 0~0.06 compositions are summarized in Figure 8h. Obviously, the x = 0.04 component exhibits a large electrostrain and strong temperature stability (E = 60 kV/cm). Within 25-180 • C, x = 0.04 sample presents a smallest η value of 31%, indicating an excellent temperature stability of electrostrain performance. A comparison of electrostrain performance with some representative lead-free piezoelectrics are present in Figure 8i. The electrostrain value fluctuates more than 100% for some BF-based ceramics, showing poor temperature stability [22,52,53]. Also, the electrostrain change exhibits a parabolic-like or monotonous decrease for BNT-and KNN-based ceramics, depending on phase transition temperature [54,55]. The enhanced temperature stability of electrostrain performance for x = 0.04 composition is due to the synergistic effect of negative temperature dependent electrostrain for relaxors with Pc phase and the positive one for ferroelectric matrix with R phase, which strongly suggest the special proportion of relaxor and ferroelectric phase can produce thermally stable electrostrain performance and thus is highly desirable for high-temperature of actuator applications. stability between this work with some representative BF- [22,52,53], BNT- [54], and KNN- [55] based ceramics.
It is well known that evolution in macroscopic performance is accompanied by changes in domain structure, and PFM images can convey domain structure information at a local scale. Figure 9a-c shows the PFM amplitude images for x = 0, 0.03, and 0.06 Figure 8. (a-g) Temperature dependent unipolar S-E curves at E = 60 kV/cm for BF-BT-xLNT ceramics; (h) evolution of unipolar S-E with temperature; (i) comparison of electrostrain thermal stability between this work with some representative BF- [22,52,53], BNT- [54], and KNN- [55] based ceramics.
It is well known that evolution in macroscopic performance is accompanied by changes in domain structure, and PFM images can convey domain structure information at a local scale. Figure 9a-c shows the PFM amplitude images for x = 0, 0.03, and 0.06 samples, respectively. The highlighted area in the amplitude image represents the piezoelectric response strength [56,57]. Consistent with the piezoelectric performance, x = 0 and 0.03 samples exhibit large bright areas, whereas x = 0.06 has only few bright areas, as indicated by the blue dashed box. Also, as shown in phase images of Figure 9d-f, broad strip-like long-range ordered macrodomains are clearly observed in the x = 0 and 0.03, and thin stripshaped short-range ordered nanodomains are observed in the dim area of x = 0.03 and 0.06. This indicates that both large-sized macrodomains and small-sized nanodomains coexist in the ferroelectric and nonergodic phase composition (x = 0.03 as an example). For relaxor phase in x = 0.06, small-sized nanodomains occupy almost the entire area. Therefore, largesized macrodomains with nanodomains blending have a positive effect on the piezoelectric effect, whereas the composition with increased concentration of nanodomains are beneficial for temperature-stable electrostrain performance.

Conclusions
BF-BT-xLNT (x = 0~0.06) ceramics are obtained via conventional solid-state synthesis. Rietveld refinement, Raman spectroscopy, and dielectric analysis show the system undergoes from R-to-Pc phase transition that is driven by LNT addition. Dielectric and impedance spectra show that the addition of LNT reduces the •• concentration and greatly increases the resistivity of BF-BT ceramics. At the same time, the piezoelectric performance is optimized at x = 0.02 composition (d33 = 97 pC/N), which is contributed by phase boundary effect and enhanced resistivity for efficiently poling. Notably, temperature stability of electrostrain is obtained for x = 0.04 composition, which is due to the synergistic effect of negative temperature dependent electrostrain for relaxor nanodomains with Pc phase and the positive one for ferroelectric bulk domains with R phase. This is certified by PFM images that nanodomains emerges from ferroelectric matrix with increasing doping content. It can be seen that the addition of LNT not only improves the resistivity and piezoelectric properties of BF-BT ceramics, but also enhances the temperature stability of electrostrain. Therefore, this work has implications for the design and application of high-performance piezoelectrics and temperature-stable electrostrain materials.

Conclusions
BF-BT-xLNT (x = 0~0.06) ceramics are obtained via conventional solid-state synthesis. Rietveld refinement, Raman spectroscopy, and dielectric analysis show the system undergoes from R-to-Pc phase transition that is driven by LNT addition. Dielectric and impedance spectra show that the addition of LNT reduces the V .. o concentration and greatly increases the resistivity of BF-BT ceramics. At the same time, the piezoelectric performance is optimized at x = 0.02 composition (d 33 = 97 pC/N), which is contributed by phase boundary effect and enhanced resistivity for efficiently poling. Notably, temperature stability of electrostrain is obtained for x = 0.04 composition, which is due to the synergistic effect of negative temperature dependent electrostrain for relaxor nanodomains with Pc phase and the positive one for ferroelectric bulk domains with R phase. This is certified by PFM images that nanodomains emerges from ferroelectric matrix with increasing doping content. It can be seen that the addition of LNT not only improves the resistivity and piezoelectric