Enhanced Energy Storage Performance in Na0.5Bi0.5TiO3-Based Relaxor Ferroelectric Ceramics via Compositional Tailoring

Owing to the high power density, excellent operational stability and fast charge/discharge rate, and environmental friendliness, the lead-free Na0.5Bi0.5TiO3 (NBT)-based relaxor ferroelectrics exhibit great potential in pulsed power capacitors. Herein, novel lead-free (1−x)(0.7Na0.5Bi0.5TiO3-0.3Sr0.7Bi0.2TiO3)-xBi(Mg0.5Zr0.5)O3 (NBT-SBT-xBMZ) relaxor ferroelectric ceramics were successfully fabricated using a solid-state reaction method and designed via compositional tailoring. The microstructure, dielectric properties, ferroelectric properties, and energy storage performance were investigated. The results indicate that appropriate Bi(Mg0.5Zr0.5)O3 content can effectively enhance the relaxor ferroelectric characteristics and improve the dielectric breakdown strength by forming fine grain sizes and diminishing oxygen vacancy concentrations. Therefore, the optimal Wrec of 6.75 J/cm3 and a η of 79.44% were simultaneously obtained in NBT-SBT-0.15BMZ at 20 °C and 385 kV/cm. Meanwhile, thermal stability (20–180 °C) and frequency stability (1–200 Hz) associated with the ultrafast discharge time of ~49.1 ns were also procured in the same composition, providing a promising material system for applications in power pulse devices.


Introduction
With the great need for green energy, electrical energy storage has become an important technology to store the energies generated from the wind, solar energy, tide, and other sources [1,2]. Nowadays, numerous energy storage devices, such as batteries, electrochemical bilayer capacitors, and dielectric capacitors, have been extensively investigated in recent decades, with dielectric capacitors being regarded as a more promising category for high power industries owing to their ultrahigh power density, excellent operational stabilities, superior charge-discharge performance (less than 100 ns), fatigue resistance, and impressive mechanical properties [3].
Theoretically, to evaluate the energy storage performance (ESP) of a dielectric material, the energy storage density (W) and energy storage efficiency (η) are two major parameters, which can be deduced based on the polarization-electric field (P-E) hysteresis loops, and expressed as follows [4].
Bi(M 1 M 2 )O 3 , NN-Bi(M 1 M 2 )O 3, etc., due to its significant impact on the improvement of DBS and lowering of the P r [16]. By doping Bi(Mg 1/3 Ta 2/3 )O 3 , the average DBS in (1−x) (0.6Na 0.5 Bi 0.5 TiO 3 -0.4Sr 0.7 Bi 0.2 TiO 3 )-xBa(Mg 1/3 Ta 2/3 )O 3 ceramics increases from 230 kV/cm to 570 kV/cm [17], indicating that Bi(M 1 M 2 )O 3 dopants can significantly increase the ESP. In this work, a compositional tailoring strategy was employed to attain superior ESP. Doping with SBT and Bi(Mg 0.5 Zr 0.5 )O 3 (BMZ) can not only enhance the relaxation degree by breaking the long-range ferroelectric order and forming the PNRs, but also improve the DBS. To achieve this strategical objective, the (1−x)(0.7Na 0.5 Bi 0.5 TiO 3 -0.3Sr 0.7 Bi 0.2 TiO 3 )-xBi(Mg 0.5 Zr 0.5 )O 3 (NBT-SBT-xBMZ) (x = 0, 0.09, 0.12, 0.15 and 0.18) ceramics were successfully fabricated using a solid-state reaction method. An enhanced ESPs, e.g., a high W rec of 6.75 J/cm 3 and a η of 79.44% along with an ultrafast discharge time of~49.1 ns, were obtained in the sample of x = 0.15. Furthermore, excellent stabilities of energy storage density versus temperature and frequency were procured for the ceramics with the same composition, providing a potential material candidate for practical applications in advanced power pulse devices.
The polycrystalline structure of the NBT-SBT-xBMZ was measured using X-ray diffractometry (XRD, Rigaku Ultima IV diffractometer) with the Cu Kα radiation. The microtopography of various specimens was characterized using a scanning electron microscope (SEM, Hitachi S-3400N) and the grain size distributions were statistically analyzed employing ImageJ software. Temperature-dependent dielectric properties and impedance properties were procured using a dielectric property measurement system (DMS-500, Partulab, Wuhan, China and E4980A, Agilent, Santa Clara, CA, USA) at 1, 10, 100 kHz, and 1 MHz from 30 • C to 350 • C with a heating rate of 3 • C/min. A Sawyer-Tower circuit (The Precision Premier II, Radiant Technologies, Inc., Albuquerque, NM, USA) was employed to measure the P-E hysteresis loops and current-electric field (I-E) loops. A DC high voltage with a rising rate of 0.1 kV/s was supplied to the sample using a high voltage signal source (Model 610E, Trek, Inc., Lockport, New York, NY, USA). The charge and discharge performance of the ceramics was measured using a well-designed charge/discharge circuit. The ceramics were polished down to 0.1 mm in thickness and sputtered with gold electrodes of 2 mm in diameter on each side were prepared for all electrical measurements. During the electrical measurements, the samples were put into silicone oil to avoid electric breakdown.

Phase Structure and Microstructure
The room-temperature XRD patterns for NBT-SBT-xBMZ specimens are presented in Figure 1a. One can observe that all compositions display a pure polycrystalline perovskite structure without a trace of a secondary phase, indicating that multi-cations have completely entered into the perovskite lattice sites, e.g., Mg 2+ , Zr 4+ , and Sr 2+ sites. In addition, the peak near 46 • shifts towards low diffraction angles, disclosing the lattice expansion due to the substitution of ions with a larger ionic radius (here, R(Mg 2+ ) = 0.072 nm and R(Zr 4+ ) = 0.072 nm) for a smaller one (R(Ti 4+ ) = 0.0605 nm). An asymmetric single-peak around 45 • gradually evolves to a weak split-peak, as shown in Figure 1b, implying the multi-phase coexistence, which can be further confirmed using the Rietveld refinement (see Figure 1c-g). Similar to other NBT-based relaxor ferroelectrics [16,18], rhombohedral phase (R phase, space group R3c), and tetragonal phase (T phase, space group P4bm) coexist for all specimens doped with BMZ at room temperature. As shown in Table 1, all the reliability factors of the weighted pattern (R wp ) and chi-square (χ 2 ) are less than 7.5% and 4, respectively, suggesting a good agreement between the measured results and the refined results using R phase and T phase two-phase mode. Furthermore, when the BMZ content increases, there is a rise in T phase fraction and a drop in R phase. To further explore the relationship between the BMZ content and the phase fraction evolved, Raman spectroscopy (see Figure S1) was employed. Based on the published data [19,20], the Raman spectrum of NBT-SBT-xBMZ ceramics can be divided into four regions within 100~1000 cm −1 in terms of the corresponding vibration modes. The regions ranging from 50 to 200, from 200 to 400, from 400 to 700, and >700 cm −1 are related to the vibration mode of A-site cations, B-O bond, oxygen octahedron (BO 6 ), and A1 + E bonds, respectively [20]. Compared with the sample of x = 0, other samples doped with BMZ present more diffuse Raman peaks in the range of 50~200 and 200~400 cm −1 , implying an enhanced disorder in both A and B sites. Meanwhile, the multi-ionic doping also leads to the increase of oxygen octahedral tilting, which can be verified by the Raman peaks, ranged from 400 to 700 cm −1 [21]. All results indicate that the short-range disorder has been enhanced via multi-ionic doping, resulting in the induced PNRs and further improved ESP [16].  The surface morphologies for NBT-SBT-xBMZ specimens are shown in Figure 2a-e. All specimens exhibit distinct grain boundary and dense grains without visible pores. The corresponding average grain sizes (AGS) shown in Figure 2f-j are statistically calculated using the ImageJ program. One can observe that the AGS for NBT-SBT-xBMZ specimens first decreases from 3.54 ± 1.37 µm (x = 0) to 1.85 ± 0.53 µm (x = 0.15), and then increases to 2.20 ± 1.37 µm when x further increases (x = 0.18 in this case), which indicates that the grain boundary migration can be inhibited and the AGS can be effectively decreased at appropriate BMZ content. This behavior can be attributed to the variation of lattice strain energy (∆G strain ), which can be determined by Equation (4) [22].
where N A is the Avogadro constant, M the Young's modulus, r d the doped ionic radius and r 0 the optimal radius of the lattice site. Owing to the large doped ionic radius (R(Mg 2+ ) = 0.072 nm and R(Zr 4+ ) = 0.072 nm), the change in ∆G strain resulting from the radius mismatch becomes increasingly large and then inhibits the grain boundary migration for appropriate BMZ contents, leading to a decrease in grain size. The fine grain size and compact structure might be beneficial to the ESP.

Dielectric Properties
The dielectric properties as a function of temperature and frequency for NBT-SBT-xBMZ ceramics measured in the temperature range from 30 to 340 °C at 1 kHz, 10 kHz, 100 kHz, and 1 MHz are shown in Figures 3a-d and S2. Compared with the pure NBT-SBT with a weak relaxor ferroelectric characteristic, the specimens doped with BMZ reveal stronger relaxor ferroelectric characteristics, i.e., more distinct frequency dispersion and a

Dielectric Properties
The dielectric properties as a function of temperature and frequency for NBT-SBT-xBMZ ceramics measured in the temperature range from 30 to 340 • C at 1 kHz, 10 kHz, 100 kHz, and 1 MHz are shown in Figures 3a-d and S2. Compared with the pure NBT-SBT with a weak relaxor ferroelectric characteristic, the specimens doped with BMZ reveal stronger relaxor ferroelectric characteristics, i.e., more distinct frequency dispersion and a broader phase transition temperature regime. In addition, two anomalies, including the low-temperature shoulder (T s ) and high-temperature peak (T m ), are found in all specimens. Accordingly, the low-temperature one is caused by the thermal evolution of R3c and P4bm PNRs, while the high-temperature one is related to the phase transition from R3c to P4bm PNRs and the thermal evolution of P4bm PNRs [23]. As revealed in Figure 3e, the T m s shift towards low temperatures from 329 to 266 • C with the increasing BMZ content, implying higher T phase content may be obtained at a certain temperature, which is consistent with the results of XRD refinement shown in Table 1. Especially for the specimens doped with BMZ, the T m s decrease at a sluggish rate. Meanwhile, the ε m s at T m s gradually decrease, and the permittivity curves become smooth, which is conducive to obtaining a high DBS and excellent temperature stability [24]. and P4bm PNRs, while the high-temperature one is related to the phase transition from R3c to P4bm PNRs and the thermal evolution of P4bm PNRs [23]. As revealed in Figure 3e, the Tms shift towards low temperatures from 329 to 266 °C with the increasing BMZ content, implying higher T phase content may be obtained at a certain temperature, which is consistent with the results of XRD refinement shown in Table 1. Especially for the specimens doped with BMZ, the Tms decrease at a sluggish rate. Meanwhile, the s at Tms gradually decrease, and the permittivity curves become smooth, which is conducive to obtaining a high DBS and excellent temperature stability [24].
Generally, the diffuseness for some relaxor ferroelectrics can be estimated using the modified Curie-Weiss law expressed in Equation (5) where T is temperature above Tm, the dielectric constant at T, C the Curie constant. The exponent , the value of which is between 1 (normal ferroelectrics) and 2 (relaxor ferroelectrics), is used to describe the diffuseness degree. As presented in Figure 3f, the diffuseness degree significantly increases from 1.38 (x = 0) to 2.00 (x = 0.18) by doping BMZ, resulting from the breaking of the correlation of adjacent dipoles and constructing of PNRs [25]. Note that the for the specimens doped with BMZ are 1.87, 1.90, 1.94, and 2.00, respectively, indicating that the behaviors close to relaxor ferroelectrics are expected in all samples, which plays a critical role in decreasing Pr, improving η, and enhancing the thermal stabilities, thus an enhanced ESP can be attained.   Generally, the diffuseness for some relaxor ferroelectrics can be estimated using the modified Curie-Weiss law expressed in Equation (5),

Ferroelectric and ESP
where T is temperature above T m , ε the dielectric constant at T, C the Curie constant. The exponent γ, the value of which is between 1 (normal ferroelectrics) and 2 (relaxor ferroelectrics), is used to describe the diffuseness degree. As presented in Figure 3f, the diffuseness degree significantly increases from 1.38 (x = 0) to 2.00 (x = 0.18) by doping BMZ, resulting from the breaking of the correlation of adjacent dipoles and constructing of PNRs [25]. Note that the γ for the specimens doped with BMZ are 1.87, 1.90, 1.94, and 2.00, respectively, indicating that the behaviors close to relaxor ferroelectrics are expected in all samples, which plays a critical role in decreasing P r , improving η, and enhancing the thermal stabilities, thus an enhanced ESP can be attained. Figure 4a presents the bipolar P-E hysteresis loops of NBT-SBT-xBMZ ceramics (x = 0, 0.09, 0.12, 0.15, and 0.18) at 150 kV/cm and 100 Hz. The sample of x = 0 exhibits strong ferroelectric characteristics with a large P m (41.48 µC/cm 2 ) and P r (27.65 µC/cm 2 ), which is unfavorable for ESP. Note that the P-E loops for the specimens doped with BMZ become slimmer. Specifically, the P m slightly decreases while the P r significantly decreases at the same electric field of 150 kV/cm. We propose a possible reason for the decrease in P m at the same electric field with the increase of BMZ. The rise in weakly polar T phase fraction is accompanied by a drop in a strongly polar R phase, which reduces the P m to a certain extent. This can be verified by the Rietveld refined results listed in Table 1. The current densityelectric field (I-E) loops at 150 kV/cm and 100 Hz are depicted in Figure 4b. Two prominent current peaks correspond to the switching of ferroelectric domains [26]. It should be noted that four current peaks were obtained in the samples doped with BMZ instead of two peaks, indicating a reversible field-induced phase transition, i.e., four compositions are in the ergodic phase. Furthermore, current peaks are gradually smooth with the increase of BMZ due to the closing to ideal relaxor ferroelectric characteristics (γ = 2), which is consistent with the results of Raman spectra and dielectric properties as a function of temperature and frequency. slimmer. Specifically, the Pm slightly decreases while the Pr significantly decreases at the same electric field of 150 kV/cm. We propose a possible reason for the decrease in Pm at the same electric field with the increase of BMZ. The rise in weakly polar T phase fraction is accompanied by a drop in a strongly polar R phase, which reduces the Pm to a certain extent. This can be verified by the Rietveld refined results listed in Table 1. The current density-electric field (I-E) loops at 150 kV/cm and 100 Hz are depicted in Figure 4b. Two prominent current peaks correspond to the switching of ferroelectric domains [26]. It should be noted that four current peaks were obtained in the samples doped with BMZ instead of two peaks, indicating a reversible field-induced phase transition, i.e., four compositions are in the ergodic phase. Furthermore, current peaks are gradually smooth with the increase of BMZ due to the closing to ideal relaxor ferroelectric characteristics (γ = 2), which is consistent with the results of Raman spectra and dielectric properties as a function of temperature and frequency. The unipolar P-E loops measured at 100 Hz under various electric fields are shown in Figure 5a. In accordance with the dielectric properties, all specimens show the nature of typical relaxor ferroelectrics, i.e., all P-E loops are slim with small Prs. As mentioned above, the ΔP and DBS play a dominant role in affecting the ESP. The corresponding ΔPs procured by P-E loops and statistic DBS fitted using a Weibull distribution are illustrated in Figure 5b,c. Similar to the previous report [27], the delayed polarization saturation is attained via doping BMZ, leading to a larger Wrec. Moreover, with the increasing applied electric field, the ΔPs continuously increases, resulting in the improved ESP, i.e., the Wrec gradually increases and the η maintains ~80% for all compositions (see Figures S3 and S4). In addition, as revealed in Figure 5c, the DBS first increases (from 165.92 kV/cm @ x = 0 to 393.39 kV/cm @ x = 0.15), then decreases with the further increase of x (358.05 kV/cm @ x = 0.18), of which the tendency is opposite with the grain size. Hence, the highest value of DBS in x = 0.15 might be caused by the fewer vacancies, associated with the fine grain sizes and the dense structure. Furthermore, the slopes (β) of Weibull distribution are over 15, indicating the excellent stability of sample quality and high reliability of statistic DBS. Based on Equations (1)-(3), the relevant ESPs are indirectly procured, and shown in Figure  5d. Owing to the largest ΔP (50.58 μC/cm 2 ) and the highest DBS (393.39 kV/cm) simultaneously achieved in the sample of x = 0.15, the optimal Wrec of 6.75 J/cm 3 with a η of 79.44% The unipolar P-E loops measured at 100 Hz under various electric fields are shown in Figure 5a. In accordance with the dielectric properties, all specimens show the nature of typical relaxor ferroelectrics, i.e., all P-E loops are slim with small P r s. As mentioned above, the ∆P and DBS play a dominant role in affecting the ESP. The corresponding ∆Ps procured by P-E loops and statistic DBS fitted using a Weibull distribution are illustrated in Figure 5b,c. Similar to the previous report [27], the delayed polarization saturation is attained via doping BMZ, leading to a larger W rec . Moreover, with the increasing applied electric field, the ∆Ps continuously increases, resulting in the improved ESP, i.e., the W rec gradually increases and the η maintains~80% for all compositions (see Figures S3 and S4). In addition, as revealed in Figure 5c, the DBS first increases (from 165.92 kV/cm @ x = 0 to 393.39 kV/cm @ x = 0.15), then decreases with the further increase of x (358.05 kV/cm @ x = 0.18), of which the tendency is opposite with the grain size. Hence, the highest value of DBS in x = 0.15 might be caused by the fewer vacancies, associated with the fine grain sizes and the dense structure. Furthermore, the slopes (β) of Weibull distribution are over 15, indicating the excellent stability of sample quality and high reliability of statistic DBS. Based on Equations (1)-(3), the relevant ESPs are indirectly procured, and shown in Figure 5d. Owing to the largest ∆P (50.58 µC/cm 2 ) and the highest DBS (393.39 kV/cm) simultaneously achieved in the sample of x = 0.15, the optimal W rec of 6.75 J/cm 3 with a η of 79.44% were also obtained for the same composition. Meanwhile, as revealed in Figure 5, the sample of x = 0.15 is superior to other reported lead-free ceramics [10,[28][29][30][31][32][33][34][35][36][37][38][39][40][41][42] in both W rec and E, meaning a promising material system was successfully fabricated in this work. sample of x = 0.15 is superior to other reported lead-free ceramics [10,[28][29][30][31][32][33][34][35][36][37][38][39][40][41][42] in both Wrec and E, meaning a promising material system was successfully fabricated in this work. Generally, the DBS of a material is determined by the intrinsic and extrinsic factors, in which the extrinsic factors are of dominance due to the ultrahigh intrinsic DBS (usually more than 1000 kV/cm), which outclasses the extrinsic ones. To explore the main reason for the optimal DBS occurred in the sample of x = 0.15 instead of 0.18, a complex impedance analyses were carried out. Figure 6a displays the Nyquist plots (Z vs. Z ) measured at 100~2 MHz and 500 • C for NBT-SBT-xBMZ ceramics (x = 0, 0.09, 0.12, 0.15, 0.18). It can be observed that a separate semicircle form presents in all impedance spectra, which indicates the same conductivity behaviors [43]. Specifically, for electroceramics, an equivalent circuit based on two parallel R and C elements in series contains grain resistance R1, grain capacitance C1, grain boundary resistance R2, and grain boundary capacitance C2, can be utilized to model the impedance results, as shown in the inset of Figure 6a [16]. The fitted results for the resistances and capacitances of grain and grain boundary are shown in Table 2. It is observed that the resistance of grain boundary outclasses that of the grain, indicating the dominant contribution of the grain boundary. Meanwhile, the maximal R2 occurs in the sample of x = 0.15. Obviously, the semicircles become tremendously enlarged by doping BMZ, implying a significant enhancement in resistance. Furthermore, as revealed in Figure 6a, the semicircles first become large and then small with the increase of BMZ, which is consistent with the trend of resistances and Weibull distributions as a function of x. The activation energy (E a ) for the conductivity can be calculated based on the Arrhenius relation (Equation (6)) [44]

Ferroelectric and ESP
where σ is the conductivity, σ 0 the pre-constant, k B the Boltzmann constant, T the absolute temperature. The E a of NBT-SBT-xBMZ ceramics is depicted in Figure 6b. One can see that the values of E a for all specimens are close to 1, indicating the conductivity at high temperatures is related to the oxygen vacancies [45]. Note that E a increases from 1.16 (@ x = 0) to 1.56 (@ x = 0.15) and then decreases to 1.41 (@ x = 0.18), revealing a similar tendency for oxygen vacancies with different BMZ contents. For the vacancies that exist in ceramics, firstly, both the oxygen vacancies and the bismuth vacancies exist. Secondly, the oxygen vacancy concentration is usually determined by the oxygen partial pressure, dopants, cooling rate from the sintering temperature, etc. Since these conditions are almost the same for four samples, the oxygen vacancy concentration can be regarded as almost the same with different BMZ contents. Thirdly, the incorporation of BMZ may produce Bi 3+ vacancies as a result of the evaporation of Bi 3+ ions during the sintering at high temperatures, which will compensate for the oxygen vacancies since their electrical signs are opposite. The bismuth vacancy concentration increases with the BMZ content. Fourthly, except for x = 0.18, which might contain too large bismuth vacancy concentration to be completely compensated by the oxygen vacancies, the oxygen vacancies dominate before x = 0.15, and the resistance increases with the increasing BMZ content, i.e., x value. At x = 0.15, the oxygen vacancy concentration is almost the same as that of bismuth vacancy, resulting in the largest impedance shown in Figure 6. Therefore, the increased impedance and the decreased vacancy concentrations are conducive to the large DBS of the ceramics.   Beyond that, the AGS also makes a critical impact on the DBS and is inversely proportional to the DBS according to Equation (7) [1].
where is an exponent ranging from 0.2 to 0.4. As mentioned above, the sample of x = 0.15 has the smallest AGS, demonstrating that a higher DBS might be achieved with the same composition. It can be concluded that the fine grain sizes, dense microstructure, smooth current peak, higher resistance, increased Ea, and diminished oxygen vacancy concentration collectively contribute to the improved DBS.
It is essential to investigate the stabilities of frequency and temperature for sample of x = 0.15, which can help confirm the applicable working temperature range and frequency range. The frequency-dependent and temperature-dependent unipolar P-E loops are shown in Figure 7a,c. Meanwhile, the corresponding ESP calculated based on Equations (1)-(3) are depicted in Figure 7b,d, respectively. All specimens show slim P-E loops at various frequencies and temperatures. Moreover, with the reduction of frequency, the P-E loops become slimmer, resulting in an increase in η and Wrec. As a result, only a variation of 12.6% in Wrec can be observed in the frequency range of 1-200 Hz, indicating good frequency stability. Interestingly, the P-E loops become significantly slim when the temperature exceeds ~60 °C, which is closed to the Ts for the sample of x = 0.15. It can be attributed to the ergodic relaxor ferroelectric state, in which the PNRs formed by doping BMZ can convert into long-range ferroelectric order under the application of a strong enough electric field, and convert back to the ergodic state, resulting in the maintained large Pm and dramatically decreased Pr [46]. Especially, despite a large variation in η from 75.6% (20 °C) Beyond that, the AGS also makes a critical impact on the DBS and is inversely proportional to the DBS according to Equation (7) [1].
where a is an exponent ranging from 0.2 to 0.4. As mentioned above, the sample of x = 0.15 has the smallest AGS, demonstrating that a higher DBS might be achieved with the same composition. It can be concluded that the fine grain sizes, dense microstructure, smooth current peak, higher resistance, increased E a , and diminished oxygen vacancy concentration collectively contribute to the improved DBS.
It is essential to investigate the stabilities of frequency and temperature for sample of x = 0.15, which can help confirm the applicable working temperature range and frequency range. The frequency-dependent and temperature-dependent unipolar P-E loops are shown in Figure 7a,c. Meanwhile, the corresponding ESP calculated based on Equations (1)- (3) are depicted in Figure 7b,d, respectively. All specimens show slim P-E loops at various frequencies and temperatures. Moreover, with the reduction of frequency, the P-E loops become slimmer, resulting in an increase in η and W rec . As a result, only a variation of 12.6% in W rec can be observed in the frequency range of 1-200 Hz, indicating good frequency stability. Interestingly, the P-E loops become significantly slim when the temperature exceeds~60 • C, which is closed to the T s for the sample of x = 0.15. It can be attributed to the ergodic relaxor ferroelectric state, in which the PNRs formed by doping BMZ can convert into long-range ferroelectric order under the application of a strong enough electric field, and convert back to the ergodic state, resulting in the maintained large P m and dramatically decreased P r [46]. Especially, despite a large variation in η from 75.6% (20 • C) to 90.92% (110 • C), only ±13.2% variation is detected in a temperature range from 20 to 180 • C, implying a good thermal stability. One can conclude that the sample of x = 0.15 possesses good thermal and frequency stabilities simultaneously. Beyond that, the T s can be shifted to lower temperatures via compositional tailoring, providing a method to obtain the superior ESP near room temperature.
to 90.92% (110 °C), only ±13.2% variation is detected in a temperature range from 20 to 180 °C, implying a good thermal stability. One can conclude that the sample of x = 0.15 possesses good thermal and frequency stabilities simultaneously. Beyond that, the Ts can be shifted to lower temperatures via compositional tailoring, providing a method to obtain the superior ESP near room temperature. Moreover, the charge-discharge properties are also of importance to practical applications. Theoretically, the current density (CD), power density (PD), discharge speed, and discharged energy storage density (Wd) can be procured according to Equations (8)-(10) [47]. Especially the discharge time (t0.9) is evaluated by the time when Wd is up to 90% of the maximal value under an electric field. Moreover, the charge-discharge properties are also of importance to practical applications. Theoretically, the current density (C D ), power density (P D ), discharge speed, and discharged energy storage density (W d ) can be procured according to Equations (8)-(10) [47]. Especially the discharge time (t 0.9 ) is evaluated by the time when W d is up to 90% of the maximal value under an electric field.
where I max is the maximal current peak in undamped i(t)~t curves, S the effective electrode area, V the effective volume of sample, E the applied electric field, R the load resistance (100 Ω in this work) and i(t) the discharge current in overdamped curves. Therefore, the charge-discharge properties for the sample of x = 0.15 were evaluated for its performance in a circuit. The undamped discharge current curves under various electric fields (20-180 kV/cm) are displayed in Figure 8a. The inset shows the variation of I max , suggesting that the highest I max reaches 17 A under 180 kV/cm. Meanwhile, as revealed in Figure 8b, a high C D of 541.40 A/cm 2 and P D of 48.73 MW/cm 3 are simultaneously achieved at the same electric field. Figure 8c presents the overdamped discharge waveforms under various electric fields (20-180 kV/cm).  Based on the classic RLC circuit, the discharge process can be modeled using a secondorder differential equation of current as a function of time. Basically, the solution of the equation is composed of an exponential function and a sinusoidal function. The former determines the amplitude of the discharged current, while the latter determines the periodical oscillation of the discharged current, depending on the load, i.e., the resistance. If the resistance is smaller, e.g., R << 2 √ L/C, the decay of the discharged current is not prominent, the oscillation is dominant. This process is usually called under-damped process. However, if the resistance is larger, the decay of discharged current becomes significant, and the oscillation of discharged current becomes unimportant, but it does exist to some extent. This is what we observed in Figure 8c [48].
One can see that t 0.9 slightly increases along with the increase of electric field and reaches~49.1 ns under 180 kV/cm, implying an ultrafast discharge time. In addition, a W d of 1.16 J/cm 3 was calculated based on the overdamped discharge waveform according to Equation (10). As presented in Figure 8e, the discharge time t 0.9 of~49.1 ns in our work is superior among the published materials [13,28,29,[49][50][51], meaning its great potential in advanced power pulse devices.

Conclusions
In conclusion, novel lead-free (1−x)(0.7Na 0.5 Bi 0.5 TiO 3 -0.3Sr 0.7 Bi 0.2 TiO 3 )-xBi(Mg 0.5 Zr 0.5 )O 3 (NBT-SBT-xBMZ) relaxor ceramics were successfully fabricated using a solid-state reaction method and designed via compositional tailoring. Appropriate Bi(Mg 0.5 Zr 0.5 )O 3 content can effectively enhance the relaxor ferroelectric behaviors and improve the DBS by forming fined grain sizes, increasing resistance, and diminishing oxygen vacancy concentrations. Therefore, the optimal W rec of 6.75 J/cm 3 and a η of 79.44% were simultaneously obtained in NBT-SBT-0.15BMZ at 20 • C and 385 kV/cm. Meanwhile, a good thermal stability (20-180 • C) and frequency stability (1-200 Hz) associated with the ultrafast discharge time of~49.1 ns were also obtained in the same composition, providing a promising material system for applications in power pulse devices.
Supplementary Materials: The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ma15175881/s1, Figure S1: Raman spectra for NBT-SBT-xBMZ ceramics; Figure S2: Dielectric constant and loss as a function of temperature and frequency for NBT-SBT-0BMZ ceramics in the temperature range of 30 to 340 • C; Figure S3: Unipolar P-E loops as a function of electric field at 100 Hz for (a) x = 0.09; (b) x = 0.12; (c) x = 0.15; (d) x = 0.18; Figure S4: The corresponding ESP as a function of electric field for P-E loops shown in Figure S3.

Conflicts of Interest:
The authors declare no conflict of interest.