Blog

Broad-high operating temperature range and enhanced energy storage performances in lead-free ferroelectrics | Nature Communications

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Nature Communications volume  14, Article number: 5725 (2023 ) Cite this article Marine Light

Broad-high operating temperature range and enhanced energy storage performances in lead-free ferroelectrics | Nature Communications

The immense potential of lead-free dielectric capacitors in advanced electronic components and cutting-edge pulsed power systems has driven enormous investigations and evolutions heretofore. One of the significant challenges in lead-free dielectric ceramics for energy-storage applications is to optimize their comprehensive characteristics synergistically. Herein, guided by phase-field simulations along with rational composition-structure design, we conceive and fabricate lead-free Bi0.5Na0.5TiO3-Bi0.5K0.5TiO3-Sr(Sc0.5Nb0.5)O3 ternary solid-solution ceramics to establish an equitable system considering energy-storage performance, working temperature performance, and structural evolution. A giant Wrec of 9.22 J cm−3 and an ultra-high ƞ ~ 96.3% are realized in the BNKT-20SSN ceramic by the adopted repeated rolling processing method. The state-of-the-art temperature (Wrec ≈ 8.46 ± 0.35 J cm−3, ƞ ≈ 96.4 ± 1.4%, 25–160 °C) and frequency stability performances at 500 kV cm−1 are simultaneously achieved. This work demonstrates remarkable advances in the overall energy storage performance of lead-free bulk ceramics and inspires further attempts to achieve high-temperature energy storage properties.

With the continuous growth of the world population and the development of the global economy and society, worldwide energy demand keeps increasing at an alarming rate. Due to the desperate issues, it is vital to exploit a variety of clean and sustainable energy sources1,2,3,4. Compared with various current energy storage and conversion devices (e.g., lithium-ion batteries, supercapacitors, solid oxide fuel cells), electrostatic capacitors made of dielectric materials have attracted ever-increasing attention up till now owing to their benefits in terms of swift charging-discharging rates, ultrahigh power density, excellent thermal stability, and prolonged storage lifespan5,6,7,8. Nonetheless, the comparatively low recoverable energy storage density (Wrec) of current dielectric ceramic capacitors had significantly hindered their practical utilizations in sophisticated electronic components and forefront pulsed power systems. Accordingly, substantial efforts have been undertaken to synergistically boost the comprehensive energy storage characteristics of dielectric materials, especially lead-free materials, to fulfill the pressing demands of electronic devices for integration, miniaturization, and environmental friendliness9,10,11,12,13.

Currently, common-utilized dielectric capacitors developed for energy storage include thin films, polymer-based thick films, and ceramic materials1,10,13,14,15,16,17,18,19. Among the candidate dielectric materials, bulk ceramics usually have low dielectric losses, high-temperature stability, and excellent fatigue resistance, enabling them to be more suitable for applications in various operational situations, such as aerospace, hybrid electrical vehicles, and electromagnetic pulse systems. Lead-free ceramics with relaxation properties for energy storage applications, for instance, BaTiO3 (BT)20,21,22, K0.5Na0.5NbO3 (KNN)5,23,24, NaNbO3 (NN)11,25,26, Bi0.5K0.5TiO3 (BKT)8,27, and Bi0.5Na0.5TiO3 (BNT)17,28,29,30-based ceramics, have been extensively investigated in past few years. BNT-based materials exhibit intrinsic large saturation polarization (Pmax) and are attributed to hybridization between Bi 6 s and O 2p orbitals31, which surpasses commonly used lead-free relaxor ferroelectric (RFE) counterparts. In light of this, BNT-based systems have received substantial attention in the field of energy storage and have been recognized as one of the most prospective eco-friendly materials for advanced pulsed power applications4,28,32. Typical BNT-based binary or ternary solid solutions have been widely studied recently, including BNT-BT, BNT-NN, BNT-SrTiO3 (BNT-ST), and BNT-BT-NN, etc33,34,35,36. In addition, Che et al. constructed a BNT-Ag(Nb0.5Ta0.5)O3 (ANT) ceramic combined with defect engineering and realized a high Wrec of 6.6 J cm−3 at 510 kV cm−1, albeit less success in energy efficiency (ƞ less than 75%)37. Thus, there is still a manifest challenge in obtaining ultrahigh energy storage density while maintaining high efficiency over a broad operating temperature in BNT-based ceramics.

Herein, we rationally design an effective strategy to maintain high energy storage performance upon a wide working temperature range guided by the phase-field method. Specifically, the manipulations of polymorphic polar nanoregions (PNRs) by constructing morphotropic phase boundary (MPB) in Bi0.5Na0.5TiO3-Bi0.5K0.5TiO3 (BNKT) binary system and incorporated Sr(Sc0.5Nb0.5)O3 (SSN) allow us to establish an equitable system considering energy storage performance, working temperature performance, and structural evolution. A gigantic Wrec of 9.22 J cm−3, and significantly enhanced energy efficiency ƞ of 96.3% at an external electric field of 535 kV cm−1 are realized in the BNKT-20SSN ceramic by the adopted repeated rolling processing (RRP) method. Encouragingly, apart from the ultrahigh Wrec and exceptional energy efficiency ƞ mentioned above, remarkable temperature-insensitive performance (Wrec ≈ 8.46 ± 0.35 J cm−3, Δƞ/ƞ ≤ 2%, 25–160 °C) and a slight fluctuation in frequency stability (Wrec ≈ 8.63 ± 0.18 J cm−3, Δƞ/ƞ < 3%, 1–100 Hz) at 500 kV cm−1 are accomplished in the BNKT-20SSN ceramic (RRP). Simultaneously, further relationships among structural construction, energy storage performance, and working temperature performance are studied in-depth to elucidate the structure-performance connection upon temperature influence.

The phase-field method is a powerful computational method to manifest the spatiotemporal evolution of microstructures and the related physical properties. The method is extensively and popularly utilized in dielectric materials, especially energy storage-related dielectric materials, to simulate the processes in terms of grain growth, solidification, thin-film deposition, crack propagation, and breakdown5,38,39. Here, we take Bi0.5Na0.5TiO3-Bi0.5K0.5TiO3 (BNT-BKT) binary system at the morphotropic phase boundary (MPB) as the initial phase and Sr(Sc0.5Nb0.5)O3 (SSN) as the guest phase to perform a two-dimensional phase-field simulation. The physical parameters and calculation parameters were set based on the 0.92BNT-0.18BKT (BNKT) system where rich polymorphic polar nanoregions (PNRs) could be observed over the coexistences of the rhombohedral (R3c) phase and the tetragonal (P4bm) phase.

The phase-field simulation results of the phase distribution as a function of the SSN phase fraction (see Fig. 1a) show that the R phase dominates in the BNKT end, albeit considerable T phase and local polar inhomogeneity exist. We employed a color scheme using RGB colors alongside vectors to depict different polarization orientations in the five left columns (vector contours). The colors visually distinguish ferroelectric domains characterized by distinct polar directions. However, distinguishing between different R values proved challenging due to the presence of eight <111> directions. To address this issue, in the rightmost column (symmetry contours), we incorporated a color bar to depict the phase contour. It is essential to note that the organization of the color correspondence between these two types of plots is unrelated and follows different schemes. With the increase in the fraction of the guest SSN phase, the fraction of the T phase dramatically increases and the phase fraction reverses at around x = 0.3. Subsequently, the polarization distribution of the external electric field as a function of the SSN fraction was calculated, as also shown in Fig. 1a. Obviously, the long-range ferroelectric (FE) order of the host phase could be disrupted, and more abundant nanodomains are produced when x increases, permitting the material to have the relaxor behavior40. The polarization directions are more disordered with less guest phase in the initial state, but it is harder to polarize fully. When the external electric field is applied and removed, the polarization exhibits a similar arranging behavior and restores its original state at the end, proving the process is reversible. We finally selected the BNKT-20SSN (x = 0.2) phase as the target composition, notwithstanding BNKT-30SSN (x = 0.3) presents a larger polarization response because the BNKT-20SSN phase displays a relatively reasonable structural construction and a considerable polar order–disorder arrangement as a relaxor ferroelectric. The initial state polarization distribution of the BNKT-20SSN phase to temperature, see Fig. 1b, strongly indicates that temperature is irrelevant to the constructed structure and could less infect the local energy barrier. The calculated polarization response to an external electric field (along the [100] axis) of the BNKT-20SSN (x = 0.2) phase is displayed in Fig. 1c. The calculated P–E loop illustrates the polar structure evolution (O1–P1–N–P2–O2) responding quickly to an external electric field; the polarization distributions are close between O1 and O2 strongly prove the PNRs are reversible, resulting in a low remanent polarization (Pr).

a Calculated microstructural evolutions under given fields with vector contours and phase distributions with symmetry contours. b Magnification into the microstructure evolution of the BNKT-20SSN ceramic at various temperatures. c Calculated P–E hysteresis loop and microstructure evolution of the BNKT-20SSN ceramic at 25 °C. Arrows capture the polarization magnitude and orientation.

Based on the phase-field simulation results above, we selected BNKT-20SSN as the target material for further study. Calcined BNKT-20SSN powders were used as the initial materials for preparing ceramic tapes by repeated rolling processing (RRP), which was employed to generate a dense structure with minimal porosity during sintering, enhancing the breakdown strength (Eb). As expected, the average grain size of the BNKT-20SSN ceramic (RRP) with compact microstructure is substantially finer than that of the sample prepared by conventional cold isostatic pressing (Supplementary Fig. 1a, b). Correspondingly, as presented in Supplementary Fig. 2, the theoretical Eb of the BNKT-20SSN ceramic (RRP) evaluated by the Weibull distribution experiments is much higher than that of the cold isostatic pressing sample, and the Weibull modulus β value of 19 (>10) was obtained, indicating the reliability of the breakdown strength34. The magnitude of Eb is an essential factor affecting the energy storage density of dielectric materials. Figure 2a shows the room-temperature P–E hysteresis loops of the BNKT-20SSN ceramic (RRP) measured from 100 kV cm−1 to the critical electric field (535 kV cm−1) at 10 Hz. Of particular noteworthy is that the P–E loops show a near-zero Pr and negligible hysteresis during the charging/discharging process, even at a high electric field of 535 kV cm−1. In addition, the I–E curve reveals that the initial broad peak (P1) corresponds to the transformation of PNRs into a highly polarized state, while the subsequent peak (P2) signifies the relaxation of the highly polarized state, see Supplementary Fig. 3. The current value remains unchanged from its original value upon unloading the positive field, indicating the absence of an irreversible component in the field-induced phase transitions of BNKT-20SSN ceramic41,42. As can be observed from Supplementary Fig. 4, both the Pmax and ΔP increase rapidly as the electric field rises. However, the Pr value remains almost constant around zero even at the critical electric field, facilitating energy storage efficiency. As a result, a giant Wrec (9.22 J cm−3) and a prominent high-level of η (96.3%) >95% are simultaneously achieved in the BNKT-20SSN ceramic (RRP) under the critical electric field of 535 kV cm−1, which outperform the majority of recently reported lead-free bulk ceramics with admirable energy storage performance, as shown in Supplementary Fig. 5. Moreover, it was commendable that the BNKT-20SSN ceramic (RRP) demonstrates an ultrahigh energy storage performance at relatively high temperatures (~150 °C), surpassing the majority of lead-free bulk ceramics and even certain MLCCs (Multi-Layer Ceramic Capacitors), see Fig. 2b5,6,8,9,17,20,23,28,29,40,43,44,45,46,47,48,49,50,51,52,53 and Supplementary Table 1. This achievement signifies the substantial potential of BNKT-20SSN ceramic (RRP) as a promising candidate for advanced high-temperature energy storage applications.

a Room-temperature P–E loops measured till the critical electric field of the BNKT-20SSN ceramic (RRP). b Comparisons of Wrec versus η (~150 °C) between our work with some recently reported lead-free bulk ceramics and certain MLCCs (Multi-Layer Ceramic Capacitors). c Time dependence of discharge energy density under different electric fields (R = 202 Ω). d Frequency-dependent, and e temperature-dependent P–E loops at 500 kV cm−1. f Calculated P–E loops and microstructure evolution as a function of E at different temperatures. g Wrec and η as a function of temperature under 500 kV cm−1. h A comparison of energy storage performances across a wide operating temperature range between our study and other reported bulk ceramics. i Comparisons of comprehensive properties (Wtot, Wrec, η, Wd, Wrec ~150 °C and η ~150 °C) between our study and other representative ceramics with excellent energy storage comprehensive performance.

Apart from the ultrahigh Wrec and η, another crucial criterion to determine the availability of high-power pulsed electronic components is the charge/discharge performance. The overdamped discharge property for the BNKT-20SSN ceramic (RRP) under various applied electric fields was measured at room temperature using a purpose-constructed resistance–inductance–capacitance (RLC) load circuit, and the regular overdamped oscillation waveforms reflect steady discharge behavior, as illustrated in Supplementary Fig. 6. Meanwhile, it can be observed that the discharge current rapidly approaches its peak (Imax = 8.2 A, at 500 kV cm−1) and then tends to diminish swiftly. Moreover, the high discharge energy density (Wd) ~5.2 J cm−3 can be liberated in a short period of time (t0.9, 90% of Wd is released) ~244 ns at 500 kV cm−1 (Fig. 2c). According to the foregoing data, the BNKT-20SSN ceramic (RRP) exhibits excellent charge/discharge characteristics, making it a promising candidate for pulsed power applications.

To ensure steady functioning in complex settings, ideal frequency reliability, and temperature stability must be guaranteed. The frequency-dependent energy storage property of BNKT-20SSN ceramic (RRP) was investigated at 500 kV cm−1, as shown in Fig. 2d, all P–E hysteresis loops are slender with practically unchanged Pmax values at various frequencies. The corresponding results are recorded in Supplementary Fig. 7, it can be seen that the BNKT-20SSN ceramic (RRP) exhibits excellent frequency reliability (Wrec ≈ 8.63 ± 0.18 J cm−3, ƞ ≈ 94.5% ± 2.4%) across the entire frequency range (1–100 Hz). Likewise, Fig. 2e gives the temperature stability of BNKT-20SSN ceramic (RRP) at 500 kV cm−1, with the attainable Pmax value marginally decreasing as the enhancement of the random electric field raises the reaction rate of the PNRs during heating. As expected, this result fits well with the P–E curves derived from phase-field simulations, as presented in Fig. 2f. Simultaneously, the corresponding microstructural evolution in the existence of the external electric field at various temperatures demonstrates that the domain size and polarization strength decrease as temperature increases (Fig. 2f and Supplementary Fig. 8). Of particular significance is that the BNKT-20SSN ceramic (RRP) features not only a wide operating temperature range (25–160 °C) but also an unprecedently high Wrec (≈ 8.46 ± 0.35 J cm−3) and ƞ (≈96.4 ± 1.4%) in contrast to certain other cutting-edge lead-free ceramics with excellent temperature stability (Fig. 2g, h)6,8,9,18,34,40,46,47. Furthermore, the comparisons of comprehensive performances (Wtot, Wrec, η, Wd, Wrec ~ 150 °C and η ~ 150 °C) between our study and recently reported lead-free systems are summarized in Fig. 2i5,6,8,40,45,47,54. It is clear that the BNKT-20SSN ceramic (RRP) covers a vast region of the radar diagram, suggesting that overall considerable improvements of energy storage performance, particularly Wrec/η at high temperatures (around 150 °C), and excellent overdamped discharge properties can be achieved simultaneously in studied samples. Within the measured cycling numbers of 1–104 at 300 kV cm−1, the values of Pmax and Pr maintain stability as well, see Supplementary Fig. 9. In light of this, the BNKT-20SSN ceramic (RRP) outperforms other state-of-the-art lead-free ceramics in terms of overall electrical characteristics, implying a significant potential for practical application in high-performance pulsed power capacitors.

Phase-field method calculation permits an efficient strategy for the subsequent sample preparation and optimization process. As indicated by the calculations, rich polymorphic polar nanoregions (PNRs) could be observed in the BNKT-20SSN system with the coexistence of the rhombohedral (R3c) phase and the tetragonal (P4bm) phase. High-resolution neutron powder diffraction (NPD) was then conducted to verify the two-phase coexistence. The Rietveld refinement over the NPD pattern gives a successful result, suggested by χ2 of 2.87 and Rwp of 5.12%, and the fitted results of NPD data are presented in Fig. 3a and Supplementary Table 255. Direct R phase and T phase fractions by the refinement are determined to be 10.1(5) wt% and 89.9(5) wt%, respectively, over the whole pattern. The atom fractions were set as constants in the refinement due to the value shifts dramatically in each refinement cycle, indicating that the cations in BNKT-20SSN form an ideal solid solution where Sr2+ and (Sc0.5Nb0.5)4+ ions may totally diffuse into the matrix of the BNKT. Besides, no additional peaks could confirm the existence of a superlattice structure. Unlike the powder X-ray diffraction (PXRD) pattern, the NPD pattern is more sensitive to cations in the structure. The M-points (two odd numbers and one even number Miller indices), R-points (three odd numbers Miller indices, especially 111), and X-points (one odd number and two even numbers Miller indices) indicate the existence of in-plane tilting of the T phase in the structure with the Glazer notion of a0a0c+. Though all three indicator points could be observed in the NPD pattern, no peak splitting exists in either point. The weak intensities could thus come from the in-plane tilting by the T phase, the local disorder-induced out-of-plane tilting by the cubic (C) phase, and/or cation-ordering over the A/B sites. However, the X-points are often the weakest intensities among the three indicator points, and the observed X-point shows a considerable intensity than the other two, particularly the M-points, suggesting the A-site disordering by the combination of Bi3+, Na+, K+, and Sr2+.

a Constant wavelength neutron powder diffraction refinement result at room temperature. b Atomic-resolution HAADF-STEM polarization vector image along [001] direction. c–e Magnification of the marked areas in (b). f Polarization magnitude mapping, and g polarization angle mapping.

Atomic-resolution spherical aberration transmission electron microscopy was used in the scanning transmission electron microscopy (STEM) mode along with a high-angle annular dark-field (HAADF) detector in order to analyze the local structure of multiphase-nanoregion coexistence. Precise atom arrangement of the 100/010 plane was captured and fitted by the 2D Gaussian function, see Fig. 3b. The arrangement presents a direct process of the formation of polarization, the coexistence of multiphase PNRs, and nanoregion distribution. The 110-plane atom arrangement could hardly distinguish the R and T phases; the 111-plane atom arrangement shows only the overlap of A and B sites; the higher Miller indices planes projection images present insufficient resolution hindered by the Moiré fringe formation. With the processing by the 2D Gaussian function, the polarization vector could be straightforwardly described by a vector from the B-site cations center to the A-site cations corner, represented by the yellow arrows in Fig. 3b. The R and T phases could be directly identified by the long magnitude arrows, whereas the C phase shows near-zero polarization. The bond lengths calculated from the B-site cations and the A-site cations along vertical and horizontal directions are noted as c and a, and the c/a ratio was then used to distinguish the R phase and the T phase, see Supplementary Fig. 10. Apparently, the captured image presents multiphase PNRs with R and T phases coexisting in the C matrix. The size of the same polarization direction and magnitude PNRs is ~2 nm. The coexistent multiphase-nanoregion destroys the long-range ferroelectric order and produces abundant small-size PNRs, consistent with those predicted by the phase-field calculation. Moreover, the randomly distributed nanoregions of different polarization directions and magnitudes could reduce polarization anisotropy. The enhanced polarization magnitude implies that the internal field is strong. The facts stated above could all lead to a more accessible and quicker polarization response to the external electric field, higher energy efficiency, and better energy storage performance. To better understand the polarizations in the sample, color-filled 2D patterns of summarized polarization directions and magnitudes are plotted in Fig. 3f, g. The patterns are highly corresponding to the phase-field calculation, signifying the reliability and the practical value of the phase-field prediction to design the dielectric materials for energy storage.

Pure BNT and (1-x)BNKT-xSSN ceramics were prepared to understand composition-driven and temperature-driven energy storage properties further. The powder X-ray diffraction (PXRD) patterns at room temperature (Supplementary Fig. 11a) show that all samples exhibit a typical BNT-based perovskite structure without impurities. The Rietveld refinement over PXRD patterns of each component was then processed by GSAS II software to determine the phases, and the results are displayed in Fig. 4a–c and Supplementary Fig. 1255. It is widely accepted that the BNT solid solution displays a rhombohedral (R3c) symmetry at room temperature31, consistent with the refinement result. The ratio of the R phase and the T phase gradually decreases with the increase of x, and the pure T phase was determined in the 0.82BNT-0.18BKT at 25SSN ceramic, identical to those from previous studies56,57,58. Temperature-dependent PXRD patterns of the BNKT-20SSN in the temperature range of 25–200 °C (Supplementary Fig. 13) indicate no peak splitting or phase change occurs, suggesting good temperature stability. Furthermore, all elements exhibit homogeneous distribution characteristics with no segregation in BNKT-20SSN ceramic, see Supplementary Fig. 14.

The Rietveld refinement of PXRD patterns of a BNT, b BNKT, and c BNKT-20SSN. Out-of-plane PFM phase images along with amplitude after polarization with different voltages and relaxation durations, d, g BNT, e, h BNKT, and f, i BNKT-20SSN ceramics. j Temperature and frequency-dependent dielectric spectra. k Unipolar P–E hysteresis loops at 100 kV cm−1. l Calculated P–E loops at 25 °C. m Weibull distribution of the breakdown strength on samples.

The genesis of optimum energy storage properties is intimately connected to the dynamic reaction of domain structures to external electric fields18,28. Here, the dynamic domain response to an applied voltage and the relaxation behaviors of the BNT, BNKT, and BNKT-20SSN ceramics were characterized by employing piezo-response force microscopy (PFM), as displayed in Fig. 4d–i, the schematics of the applied voltage and region are similar to our previous work (Supplementary Fig. 15, 16)40. For the BNT ceramics (Fig. 4d, g), significant phase deviations and amplitude signals were observed even at a lower external voltage (10 V), and the overwhelming majority of domains did not flip after a 20 min relaxation time, demonstrating a robust FE characteristic with substantial Pmax and Pr. A similar occurrence was also seen in the BNKT ceramic (Fig. 4e, h), albeit the feedback signal (after 20 min) of the BNKT ceramic is slightly weaker than that of the BNT at an applied voltage of 10 V. On the contrary, the phase difference and amplitude signals of the BNKT-20SSN ceramic (Fig. 4f, i) are weak in the entire electrical scan zone and become more potent at higher external voltages. After a 20-min relaxation duration, the signals have significantly decreased. The facts mentioned above could all lead to a negligible Pr and hysteresis loss in the BNKT-20SSN ceramic, thus increasing the energy storage efficiency.

In order to further investigate the mechanism of relaxation polarization behavior over a wide range of temperatures, the temperature (from room temperature to 400 °C) and frequency (from 1 kHz to 1 MHz) dependent dielectric permittivity of pure BNT and (1-x)BNKT-xSSN ceramics were performed and depicted in Fig. 4j and Supplementary Fig. 17. The pure BNT ceramic possesses two typical temperature-dependent dielectric anomalous peaks, labeled as “Ts” (depolarization temperature shoulder) and “Tm” (temperature of the maximum dielectric constant) at low and high temperatures, respectively59,60,61,62. At the ambient temperature (Tamb), as displayed in Supplementary Fig. 17a, the pure BNT ceramic exhibits a non-ergodic relaxor state and behaves as a normal ferroelectric with a R3c symmetry space group31,35,63. Afterward, with the addition of BKT (Supplementary Fig. 17b), the BNKT ceramic with R and T phase coexistence at MPB displays a much higher maximum permittivity value (~ 4600 at 1 MHz) than that of the pure BNT ceramic (~2500 at 1 MHz). The Tm of the BNKT ceramic shifts to a lower temperature, and a relatively broad permittivity peak can be detected in the MPB composition64. Subsequently, with increasing the content of SSN, the maximum value of the dielectric permittivity decreases step by step, see Fig. 4j. The two dielectric anomaly humps squint towards lower temperatures, and a gradually broaden temperature plateau of the dielectric constant can be observed. The analysis results of ln(1/ε´−1/εm´) versus ln(T – Tm) for BNT and (1-x)BNKT-xSSN ceramics are presented in Supplementary Fig. 18. It is observed that the calculated values of γ for (1-x)BNKT-xSSN ceramics lie within the range of 1.86−1.99, indicating the pronounced relaxor behavior. Noteworthy, the dielectric relaxor or frequency dispersion behavior of the BNKT-20SSN ceramic could be noticed around Tamb, which allows for the formation of relaxor ferroelectric (a significant polar order–disorder arrangement) at room temperature and is beneficial to upgrading the energy storage performances32.

The bipolar P–E hysteresis loops of all components mentioned above were performed under 100 kV cm−1 at room temperature, and the results are shown in Fig. 4k. The pure BNT ceramic exhibits a canonical FE behavior with large Pmax, Pr, considerable hysteresis, and a fully saturated P–E loop. A lower coercive electric field (Ec) and a slightly improved Pmax were obtained in the BNKT (MPB) ceramic, indicating a softened FE activity. The hysteresis losses decrease dramatically, accompanied by slimmer P–E loops seen in the SSN-introduced compositions, and are consistent with those calculated from the phase-field simulation method (Fig. 4l). To better understand composition-driven breakdown strength, the values of theoretical Eb of all ceramic samples were evaluated by the Weibull distribution experiments, see Fig. 4m. With increasing the SSN concentration, the theoretical Eb values first increased and then decreased, reaching a maximum at the BNKT-20SSN ceramics. To investigate the energy storage performance of the ceramics, the P–E hysteresis loops were performed on a ferroelectric analyzer. As depicted in Supplementary Fig. 19, the room-temperature P–E loops of the BNKT-20SSN ceramic were measured from 100 kV cm−1 to the critical electric field (385 kV cm−1) at 10 Hz. As a consequence, a great Wrec of 5.23 J cm−3 along with a high energy efficiency of 90.2% are simultaneously achieved in the BNKT-20SSN ceramic, demonstrating a substantial promotional impact of combinatorial-optimization strategy guided by phase-field simulations on energy storage performance.

The selection of promising dielectric materials to combine is a foundational but inevitable procedure in material science. The presented phase-field method assistant strategy could be used to investigate dielectric material with high energy storage performance upon a wide working temperature range5,40. By learning the structural evolution, polarization distribution, and characteristics of PNRs in combinations of BNT-BKT-SSN at morphotropic phase boundary that affords possibly engineerable ceramic materials from the phase-field calculation simulation strategy, we could initiatively decide and select the optimal composition for further study to achieve the remarkable performances.

To state, various combinations in perovskite structure have been considered to acquire substantial polar nanoregions in the structure to enhance energy storage performances. The initial BNT ceramic shows a considerable ferroelectric performance, but the Pr is too large for energy storage applications31,32,37; following the BNKT binary system at the morphotropic phase boundary is an ideal fundament for massive afterward studies56,61. We intentionally choose ionic radius and valence distinguishable cation pairs as guests for targeting a solid solution and thus enable the discovery of diverse ferroelectric/energy storage performances. The phase-field method supports an efficient strategy of rational design of a targeted material with relaxor behavior by implying reduced PNRs, thermal stability, reversible response to the external electric field, high Pmax, and low Pr, all of which could effectively optimize the ferroelectric/energy storage performances. The designed BNKT-20SSN ceramic (RRP) shows extraordinary performances of giant Wrec and ultrahigh η, unprecedented temperature, and frequency stabilities. Multi-methods analysis describes the occurrence of the performances. The small, uniform, dense grain structures by the RRP process method guarantee the breakdown strength and subsequent applications. The coexisting R/T/C multiphase nanoregions confirmed by atomic-scale TEM provide small-scaled PNRs with different polarization distributions in both magnitude and orientation, suggesting a relaxor-behavior FE with large Pmax of energy storage density and small Pr of energy loss and allowing the polarization of nanoregions flipping along the same directions with lower energy and higher dynamics. The PNRs of the BNKT ceramic were confirmed with a size of ~200 nm, and the guest-invited BNKT-20SSN reduced the size to ~2 nm65. Notably, typical relaxor REs present only minor differences in polarization magnitude but only in orientation, which differs from what we reported here1,66. This phenomenon may be because of the considerable amount of SSN concentration, x = 0.20, and the fact that the A site is composed of 0.42Bi-0.34Na–0.08K–0.16 Sr, introducing considerably strong local heterogeneity and destabilizing the polarization ordering. The system retains the original PNRs of BNKT, but the PNRs are reduced and dissociated, forming a local order–disorder coexistence system. The metastable order–disorder system in BNKT-20SSN presents a lower Pmax than BNKT, but the abundant PNRs promote the performance of energy storage. Moreover, BNKT-20SSN ceramic (RRP) exhibits unprecedented temperature and frequency stabilities where the Wrec and η typically decrease when temperature increases, indicated by the decrease of Pmax and the occurrence of the non-negligible Pr. This case may come from a series of benefits. First, the average structure is stable upon temperature change to 200 °C. Second, the heterogeneity in the metastable order–disorder system may act as a pinning effect and hinder further atom thermal vibration-induced local structure change67,68. Finally, the randomly distributed PNRs also show some order–disorder, inviting a suitable local energy barrier to the external electric field and frequencies, indicated by the flat temperature-dependent dielectric response.

In summary, a rational phase-field calculation guided dielectric material designing method is presented to achieve small-scale and stable PNRs with overall well-behaved energy storage performance theoretically and experimentally. The so-calcined BNKT-20SSN guided by calculations with improved RRP techniques shows a Wrec of 9.22 J cm−3 and an ultrahigh ƞ of ~96.3% at large external electric field 535 kV cm−1, accompanied with excellent temperature-insensitive (Wrec ≈ 8.46 ± 0.35 J cm−3, ƞ ≈ 96.4 ± 1.4%, 25–160 °C) and frequency stability (Wrec ≈ 8.63 ± 0.18 J cm−3, ƞ ≈ 94.5% ± 2.4%, 1–100 Hz) at 500 kV cm−1. This research provides a paradigm for the synergistic development of lead-free dielectric materials with enhanced comprehensive energy storage capacity over a broad operating temperature range to fulfill the pressing demands of modern energy storage components.

The lead-free ceramics with the composition of (1-x)Bi0.5(Na0.82K0.18)0.5TiO3-xSr(Sc0.5Nb0.5)O3 (abbreviated as BNKT-xSSN, where x = 0, 0.05, 0.10, 0.15, 0.20, and 0.25) were fabricated via a conventional solid-state reaction method. The pre-dried oxides and carbonate powders of Bi2O3 (≥99.9%), Na2CO3 (≥99.8%), K2CO3 (≥99.5%), TiO2 (≥99.9%), SrCO3 (≥99.5%), Sc2O3 (≥99.9%) and Nb2O5 (≥99.99%) (Aladdin Chemical Reagent Co., Ltd., CN) were used as the starting materials. The raw powders were weighed according to the stoichiometric ratio, and 2 wt.% excess of Bi2O3, Na2CO3, and K2CO3 were added to compensate for the weight loss at high temperatures. The mixtures were planetarily ball-milled for 24 h at 300 rpm with ethanol in nylon jars using Y-stabilized zirconia balls as milling media. After drying, the mixtures were calcined at 800–850 °C for 4 h in an alumina crucible and ball-milled again similarly. After drying and sieving, the powders were uniaxially pressed into pellets at 10 MPa in a 12 mm diameter stainless steel cylindrical die. Then, the preformed green ceramic sheets were pressed at 300 MPa for 3 min through cold isostatic pressing. Finally, the samples were sintered at 1140–1200 °C for 4 h with a heating rate of 4 °C min−1 in closed alumina crucibles. To minimize the evaporation of the volatile elements Bi, Na, and K, the samples were embedded in powders of the same composition. For the electrical measurements, the silver paste was coated on both surfaces of the sintered samples and fired at 650 °C for 20 min to form electrodes.

The green ceramic tapes with BNKT-20SSN composition were prepared by a repeated rolling processing (RRP) method, in which the pre-sintered powders were used as raw materials. After sieving, the calcined powders of BNKT-20SSN were thoroughly mixed with polyvinyl alcohol (PVA). The extrusion force generated by the two tightly contiguous rollers could greatly improve the density of the green ceramic tape and also help to form a dense structure during sintering with minimal porosity, resulting in a significant increase in Eb. And then, the ceramic tapes prepared by RRP were punched into disks with a diameter of 12 mm and then calcined at 650 °C for 2 h to burn out the PVA. Subsequently, the green disks were embedded in the pre-sintered powders of the same composition and sintered at 1070–1140 °C for 2 h using a double crucible method, so that the volatilization of Bi, Na, and K elements during high-temperature sintering could be suppressed.

The phase purity and crystalline structure were characterized at an X-ray diffractometer (SmartLab-3 kW, Rigaku, Tokyo, Japan) with Cu Kα radiation. Neutron diffraction patterns were collected on the high-resolution powder diffractometer ECHIDNA at ANSTO over the angular range 8 ≤ 2θ/° ≤ 160, using a step size Δ2θ = 0.05° and a wavelength of 1.622 Å at room temperature. The Rietveld refinements on PXRD and NPD were analyzed using the program GSAS II55. The microstructure of the ceramics was observed by field-emission scanning electron microscopy (FE-SEM, FEI Quanta 250 FEG, Hillsboro, Oregon, USA). The Nano Measurer software was used to calculate the average grain size of the samples based on the SEM images. The HAADF atomic-scale images were acquired using an atomic-resolution STEM (aberration-corrected Titan Themis G2 microscope) and processed by 2D Gaussian fitting in MATLAB scripts to evaluate the polarization vector, magnitude, and angle maps.

For the dielectric measurements, the silver paste was coated on both surfaces of the sintered samples and fired at 650 °C for 20 min to form electrodes. Temperature and frequency-dependent dielectric permittivity and tanδ were tested by an accurate inductor–capacitance–resistance (LCR) meter (E4980A, Agilent, USA) with a 3 °C min−1 heating rate.

The room-temperature P–E loops with a triangle signal of 10 Hz and temperature and frequency-dependent P–E loops were measured on a ferroelectric analyzer (aix ACCT, TF Analyzer 2000, Aachen, Germany). For the P–E loops measurements, the sintered disks were polished to ≈35 μm and then sputtered with gold electrodes on both surfaces.

The change–discharge properties of ceramics with a thickness of ≈50 μm were investigated via a specially designed resistance, inductance, and capacitance (RLC) load circuit.

The breakdown performances were measured via a ferroelectric analyzer (aix ACCT, TF Analyzer 2000, Aachen, Germany). The statistical Eb properties for perovskite structure ceramics are generally estimated by Weibull distribution equations: Xi = ln(Ei), Yi = ln(ln(1/(1-Pi))), Pi = i/(n + 1), where Ei is the experimental breakdown strength of each sample in ascending order, Pi is the cumulative probability of dielectric breakdown, n represents the total number of samples for ceramics. The Weibull distribution modulus, which can estimate the value of Eb, was represented by the slope of the fitting line.

The surface morphology, dynamic response, and relaxation behavior of the domains under external voltages were investigated by piezoelectric force microscopy (PFM, MFP-3D, Asylum Research, USA). The surfaces of each ceramic specimen used for PFM measurement were finely polished with polycrystalline diamond polishing paste.

Phase-field simulations have been employed to investigate a single crystal undergoing the Cubic (C) to Tetragonal (T) to Rhombohedral (R) ferroelectric transition. This study encompasses varying defect doping concentrations in the range of x = 0 to 0.30, denoted as BNKT-xSSN. The comprehensive evaluation of the ferroelectric system’s total free energy is expressed as follows:

where fbulk represents the bulk free energy density,

where αij represents the coefficient, the value of which is contingent upon both concentration c and the temperature T.

fgrad represents the gradient energy density.

where G11 signifies the coefficient for gradient energy, while fcouple pertains to the dipole effect induced by doping.

\({f}_{{{{{{\rm{couple}}}}}}}=-\int {d}^{3}x{\sum}_{i={{{{\mathrm{1,2,3}}}}}}{P}_{i}(x)\cdot {\varphi }_{{loc}}(x)\) , where \({\varphi }_{{loc}}(x)\) is the dipolar field generated through doping, expected to exhibit a random distribution and remain unchanged during the cooling process. felas corresponds to the energy density associated with long-range elastic interactions, and felec pertains to the energy density attributed to electrostatic interactions. \({f}_{{elas}}=\frac{1}{2}{c}_{{ijkl}}{e}_{{ij}}{e}_{{kl}}=\frac{1}{2}{c}_{{ijkl}}({\varepsilon }_{{ij}}-{\varepsilon }_{{ij}}^{0})({\varepsilon }_{{kl}}-{\varepsilon }_{{kl}}^{0})\) , where cijkl denotes the elastic constant tensor, εij signifies the total strain, and ε0kl represents the electrostrictive stress-free strain, i.e., ε0kl = QijklPkPl. felec = fdipole + fdepola + fappl, where fdipole is the dipole-dipole interactions arising from polarization, fdepola, represents the depolarization energy density, and fappl, denoting the energy density resulting from an applied electric field. The temporal evolution of the spontaneous polarization field (P) can is determined by solving the time-dependent Ginzburg-Landau (TDGL) equation: \(\frac{{{dP}}_{i}(x,{t})}{{dt}}=-\!M\frac{\delta F}{\delta {P}_{i}(X,t)}\) , i = 1, 2, 3, where M represents the kinetic coefficient, F signifies the total free energy, and t denotes time5,69.

All data supporting this study and its findings are available within the article and its Supplementary Information. The data that support the findings of this study are available on request from the corresponding authors.

Pan, H. et al. Ultrahigh energy storage in superparaelectric relaxor ferroelectrics. Science 374, 100–104 (2021).

ADS  CAS  PubMed  Google Scholar 

Liu, C. et al. Advanced materials for energy storage. Adv. Mater. 22, E28–E62 (2010).

Yang, L. et al. Perovskite lead-free dielectrics for energy storage applications. Prog. Mater. Sci. 102, 72–108 (2019).

Yan, F. et al. Boosting energy storage performance of lead-free ceramics via layered structure optimization strategy. Small 18, 2202575 (2022).

Li, D. et al. Improved energy storage properties achieved in (K, Na)NbO3‑based relaxor ferroelectric ceramics via a combinatorial optimization strategy. Adv. Funct. Mater. 32, 2111776 (2021).

Luo, N. et al. Constructing phase boundary in AgNbO3 antiferroelectrics: pathway simultaneously achieving high energy density and efficiency. Nat. Commun. 11, 4824 (2020).

ADS  CAS  PubMed  PubMed Central  Google Scholar 

Xie, A. et al. Supercritical relaxor nanograined ferroelectrics for ultrahigh-energy-storage capacitors. Adv. Mater. 34, 2204356 (2022).

Chen, L. et al. Outstanding energy storage performance in high-hardness (Bi0.5K0.5)TiO3-based lead-free relaxors via multi-scale synergistic design. Adv. Funct. Mater. 32, 2110478 (2021).

Chen, L. et al. Giant energy-storage density with ultrahigh efficiency in lead-free relaxors via high-entropy design. Nat. Commun. 13, 3089 (2022).

ADS  CAS  PubMed  PubMed Central  Google Scholar 

Li, Q. et al. Flexible high-temperature dielectric materials from polymer nanocomposites. Nature 523, 576–579 (2015).

ADS  CAS  PubMed  Google Scholar 

Qi, H. et al. Ultrahigh energy-storage density in NaNbO3-based lead-free relaxor antiferroelectric ceramics with nanoscale domains. Adv. Funct. Mater. 29, 1903877 (2019).

Pan, Z. et al. Fatigue-free Aurivillius phase ferroelectric thin films with ultrahigh energy storage performance. Adv. Energy Mater. 10, 2001536 (2020).

Pan, H. et al. Ultrahigh-energy density lead-free dielectric films via polymorphic nanodomain design. Science 365, 578–582 (2019).

ADS  CAS  PubMed  Google Scholar 

Pan, H. et al. Giant energy density and high efficiency achieved in bismuth ferrite-based film capacitors via domain engineering. Nat. Commun. 9, 1813 (2018).

ADS  PubMed  PubMed Central  Google Scholar 

Wang, Y. et al. Gradient-layered polymer nanocomposites with significantly improved insulation performance for dielectric energy storage. Energy Storage Mater. 24, 626–634 (2020).

Guo, Y. et al. Enhancement in the energy storage performance of P(VDF-HFP)-based composites by adding PLZST inorganic nanoparticles. J. Mater. Chem. A 10, 22058–22066 (2022).

Li, J. et al. Grain-orientation-engineered multilayer ceramic capacitors for energy storage applications. Nat. Mater. 19, 999–1005 (2020).

ADS  CAS  PubMed  Google Scholar 

Yan, F. et al. Composition and structure optimized BiFeO3-SrTiO3 lead-Free ceramics with ultrahigh energy storage performance. Small 18, 2106515 (2022).

Lu, Z. et al. Superior energy density through tailored dopant strategies in multilayer ceramic capacitors. Energy Environ. Sci. 13, 2938–2948 (2020).

Hu, Q. et al. Achieve ultrahigh energy storage performance in BaTiO3-Bi(Mg1/2Ti1/2)O3 relaxor ferroelectric ceramics via nano-scale polarization mismatch and reconstruction. Nano Energy 67, 104264 (2020).

Li, W. et al. High-temperature BaTiO3-based ternary dielectric multilayers for energy storage applications with high efficiency. Chem. Eng. J. 414, 128760 (2021).

Zhao, P. et al. High‐performance relaxor ferroelectric materials for energy storage applications. Adv. Energy Mater. 9, 1803048 (2019).

Zhang, M. et al. Significant increase in comprehensive energy storage performance of potassium sodium niobate-based ceramics via synergistic optimization strategy. Energy Storage Mater. 45, 861–868 (2022).

Yang, Z. et al. Grain size engineered lead-free ceramics with both large energy storage density and ultrahigh mechanical properties. Nano Energy 58, 768–777 (2019).

Xie, A. Achieving stable relaxor antiferroelectric P phase in NaNbO3-based lead-free ceramics for energy-storage applications. J. Materiomics 8, 618–626 (2022).

Kang, R. et al. Extraordinary energy storage performance and thermal stability in sodium niobate-based ceramics modified by the ion disorder and stabilized antiferroelectric orthorhombic R phase. J. Mater. Chem. A 9, 24387–24396 (2021).

Yang, Q. et al. Excellent energy storage performance of K0.5Bi0.5TiO3-based ferroelectric ceramics under low electric field. Chem. Eng. J. 414, 128769 (2021).

Zhu, X. et al. Ultrahigh energy storage density in (Bi0.5Na0.5)0.65Sr0.35TiO3-based lead-free relaxor ceramics with excellent temperature stability. Nano Energy 98, 107276 (2022).

Li, J. et al. Multilayer lead-free ceramic capacitors with ultrahigh energy density and efficiency. Adv. Mater. 30, 1802155 (2018).

Li, C. et al. Superior energy storage performance in (Bi0.5Na0.5)TiO3-based lead-free relaxor ferroelectrics for dielectric capacitor application via multiscale optimization design. J. Mater. Chem. A 10, 9535–9546 (2022).

Zhou, X. et al. Phase structure and properties of sodium bismuth titanate lead-free piezoelectric ceramics. Prog. Mater. Sci. 122, 100836 (2021).

Yan, F. et al. Superior energy storage properties and excellent stability achieved in environment-friendly ferroelectrics via composition design strategy. Nano Energy 75, 105012 (2020).

Chu, B. et al. High-energy storage properties over a broad temperature range in La-modified BNT-based lead-free ceramics. ACS Appl. Mater. Interfaces 14, 19683–19696 (2022).

Chen, L. et al. Local diverse polarization optimized comprehensive energy-storage performance in lead-free superparaelectrics. Adv. Mater. 34, 2205787 (2022).

Qi, H. et al. Emerging antiferroelectric phases with fascinating dielectric, polarization and strain response in NaNbO3-(Bi0.5Na0.5)TiO3 lead-free binary system. Acta Mater. 208, 116710 (2021).

Liu, G. et al. Phase evolution in (1-x)(Na0.5Bi0.5)TiO3-xSrTiO3 solid solutions: a study focusing on dielectric and ferroelectric characteristics. J. Materiomics 6, 677–691 (2020).

Che, Z. et al. Phase structure and defect engineering in (Bi0.5Na0.5)TiO3-based relaxor antiferroelectrics toward excellent energy storage performance. Nano Energy 100, 107484 (2022).

Shen, Z. et al. High-throughput phase-field design of high-energy-density polymer nanocomposites. Adv. Mater. 30, 1704380 (2018).

Gao, J. et al. Dewetting of Ni silicide thin film on Si substrate: In-situ experimental study and phase-field modeling. Acta Mater. 223, 117491 (2022).

Li, D. et al. Lead-free relaxor ferroelectric ceramics with ultrahigh energy storage densities via polymorphic polar nanoregions design. Small 19, 2206958 (2022).

Wu, J. et al. Ultrahigh field-induced strain in lead-free ceramics. Nano Energy 76, 105037 (2020).

Zhang, M. et al. Phase transformations in an Aurivillius layer structured ferroelectric designed using the high entropy concept. Acta Mater. 229, 117815 (2022).

Yan, F. et al. Gradient-structured ceramics with high energy storage performance and excellent stability. Small 19, 2206125 (2023).

Yang, W. et al. Superior energy storage properties in NaNbO3-based ceramics via synergistically optimizing domain and band structures. J. Mater. Chem. A 10, 11613–11624 (2022).

Qi, H. et al. Superior energy-storage capacitors with simultaneously giant energy density and efficiency using nanodomain engineered BiFeO3-BaTiO3-NaNbO3 lead-free bulk ferroelectrics. Adv. Energy Mater. 10, 1903338 (2019).

Yuan, Q. et al. Bioinspired hierarchically structured all-inorganic nanocomposites with significantly improved capacitive performance. Adv. Funct. Mater. 30, 2000191 (2020).

Xie, A. et al. NaNbO3-(Bi0.5Li0.5)TiO3 lead-free relaxor ferroelectric capacitors with superior energy-storage performances via multiple synergistic design. Adv. Energy Mater. 11, 2101378 (2021).

Cai, Z. et al. High-temperature lead-free multilayer ceramic capacitors with ultrahigh energy density and efficiency fabricated via two-step sintering. J. Mater. Chem. A 7, 14575–14582 (2019).

Wang, D. et al. Bismuth ferrite-based lead-free ceramics and multilayers with high recoverable energy density. J. Mater. Chem. A 6, 4133–4144 (2018).

Wang, G. et al. Ultrahigh energy storage density lead-free multilayers by controlled electrical homogeneity. Energy Environ. Sci. 12, 582–588 (2019).

Zhu, L. et al. Heterovalent-doping-enabled atom-displacement fluctuation leads to ultrahigh energy-storage density in AgNbO3-based multilayer capacitors. Nat. Commun. 14, 1166 (2023).

ADS  CAS  PubMed  PubMed Central  Google Scholar 

Zhao, P. et al. Ultrahigh energy density with excellent thermal stability in lead-free multilayer ceramic capacitors via composite strategy design. J. Mater. Chem. A 9, 25914–25921 (2021).

Wang, G. et al. Fatigue resistant lead-free multilayer ceramic capacitors with ultrahigh energy density. J. Mater. Chem. A 8, 11414–11423 (2020).

Hu, D. et al. Optimization the energy density and efficiency of BaTiO3-based ceramics for capacitor applications. Chem. Eng. J. 409, 127375 (2021).

Toby, B. et al. GSAS-II: the genesis of a modern open-source all purpose crystallography software package. J. Appl. Cryst. 46, 544–549 (2013).

Jing, R. et al. Phase evolution and relaxor to ferroelectric phase transition boosting ultrahigh electrostrains in (1-x)(Bi1/2Na1/2)TiO3-x(Bi1/2K1/2)TiO3 solid solutions. J. Materiomics 8, 335–346 (2022).

Li, T. et al. Giant strain with low hysteresis in A-site-deficient (Bi0.5Na0.5)TiO3-based lead-free piezoceramics. Acta Mater. 128, 337–344 (2017).

Yin, J. et al. Ultrahigh strain in site engineering-independent Bi0.5Na0.5TiO3-based relaxor-ferroelectrics. Acta Mater. 147, 70–77 (2018).

Ma, C. et al. Creation and destruction of morphotropic phase boundaries through electrical poling: a case study of lead-free (Bi1/2Na1/2)TiO3-BaTiO3 piezoelectrics. Phys. Rev. Lett. 109, 107602 (2012).

Matsuura, M. et al. Damped soft phonons and diffuse scattering in (Bi1/2Na1/2)TiO3. Phys. Rev. B 87, 064109 (2013).

Hu, D. et al. Greatly enhanced discharge energy density and efficiency of novel relaxation ferroelectric BNT-BKT-based ceramics. J. Mater. Chem. C. 8, 591–601 (2020).

Wang, M. et al. Ultrahigh energy sorage density and efficiency in Bi0.5Na0.5TiO3-based ceramics via the domain and bandgap engineering. ACS Appl. Mater. Interfaces 13, 51218–51229 (2021).

Schütz, D. et al. Lone-pair-induced covalency as the cause of temperature- and field-induced instabilities in bismuth sodium titanate. Adv. Funct. Mater. 22, 2285–2294 (2012).

Hiruma, Y. et al. Phase transition temperature and electrical properties of (Bi1/2Na1/2)TiO3-(Bi1/2A1/2)TiO3 (A = Li and K) lead-free ferroelectric ceramics. J. Appl. Phys. 103, 084121 (2008).

Otonicar, M. et al. External-field-induced crystal structure and domain texture in (1-x)Na0.5Bi0.5TiO3-xK0.5Bi0.5TiO3 piezoceramics. Acta Mater. 127, 319–331 (2017).

Li, F. et al. Ultrahigh piezoelectricity in ferroelectric ceramics by design. Nat. Mater. 17, 349–354 (2018).

ADS  CAS  PubMed  Google Scholar 

Tuckmantel, P. et al. Local probe comparison of ferroelectric switching event statistics in the creep and depinning regimes in Pb(Zr0.2Ti0.8)O3 thin films. Phys. Rev. Lett. 126, 117601 (2021).

ADS  CAS  PubMed  Google Scholar 

Huang, Q. et al. Manipulating ferroelectric behaviors via electron-beam induced crystalline defects. Nanoscale 13, 14330–14336 (2021).

Wang, D. et al. Phase diagram of polar states in doped ferroelectric systems. Phys. Rev. B 86, 054120 (2012).

This work was supported by the National Key R&D Program of China (2021YFB3800602), the International Cooperation Project of Shaanxi Province (2021KWZ-10), the Fundamental Research Funds for the Central University, the 111 Project of China (B14040), Zhejiang Provincial Science and Technology Program under Grant LGG20F010007. The SEM work was done at the International Center for Dielectric Research (ICDR), Xi’an Jiaotong University, Xi’an, China. The authors thank Dr. Yan-Zhu Dai for her help with SEM.

Electronic Materials Research Laboratory & Multifunctional Materials and Structures, Key Laboratory of the Ministry of Education & International Center for Dielectric Research, School of Electronic Science and Engineering, Xi’an Jiaotong University, 710049, Xi’an, Shaanxi, China

Weichen Zhao, Diming Xu, Da Li, Yan Guo & Di Zhou

Australian Nuclear Science and Technology Organization, Lucas Heights, 2234, NSW, Australia

School of Physics and Information Technology, Shaanxi Normal University, 710062, Xi’an, Shaanxi, China

State Key Laboratory for Strength and Vibration of Mechanical Structures, School of Aerospace, Xi’an Jiaotong University, 710049, Xi’an, Shaanxi, China

Department of Chemistry, Zhejiang University, 310027, Hangzhou, Zhejiang, China

School of Electronic and Information Engineering, Hangzhou Dianzi University, 310018, Hangzhou, Zhejiang, China

State Key Laboratory of Electrical Insulation and Power Equipment, Xi’an Jiaotong University, 710049, Xi’an, Shaanxi, China

Frontier Institute of Science and Technology and State Key Laboratory for Mechanical Behavior of Materials, Xi’an Jiaotong University, 710049, Xi’an, Shaanxi, China

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

The work was conceived and designed by W.C.Z., D.Z., and W.C.Z. fabricated the samples, tested the energy storage, dielectric, structure, stability, and other properties, and processed related data, assisted by D.Z., D.X., and D.L. The SEM images were filmed and processed by G.Y. and T. Z. The HAADF-STEM images were filmed and processed by H.M.J. and D.X. The PFM images were filmed and processed by M.K.X. and W.F.L. The neutron diffraction data were processed and analyzed by D.X. and M.A. Temperature-dependent XRD data were processed and analyzed by D.E.S. The manuscript was drafted by W.C.Z. and revised by D.Z. and D.X. The phase-field simulations were performed by D.W. and discussed with D.Z. and W.C.Z. All authors participated in the data analysis and discussions.

Correspondence to Diming Xu, Dong Wang or Di Zhou.

The authors declare no competing interests.

Nature Communications thanks Xin Chen, Haixue Yan and the other anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Zhao, W., Xu, D., Li, D. et al. Broad-high operating temperature range and enhanced energy storage performances in lead-free ferroelectrics. Nat Commun 14, 5725 (2023). https://doi.org/10.1038/s41467-023-41494-1

DOI: https://doi.org/10.1038/s41467-023-41494-1

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

Journal of Materials Science: Materials in Electronics (2024)

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Nature Communications (Nat Commun) ISSN 2041-1723 (online)

Broad-high operating temperature range and enhanced energy storage performances in lead-free ferroelectrics | Nature Communications

High rate battery Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.