High energy storage properties of (Nb0.5La0.5)4+ complex-ion modified (Ba0.85Ca0.15)(Zr0.10Ti0.90)O3 ceramics
Abstract
In this work, a traditional solid-state method was adopted to prepare dense Ba0.85Ca0.15[(Zr0.1Ti0.9)1−x–(Nb0.5La0.5)x]O3(BCZT–xNL) lead-free relaxation ferroelectrics with excellent energy storage performance (ESP). Remarkably, a large breakdown field strength of BDS (∼410kV/cm) resulting from decreased grain size by double ions doping at the B-site was achieved in BCZT–xNL ceramics. Especially, the BCZT–0.12NL ceramic displays a good recoverable energy density (Wrec) of 3.1J/cm3, a high efficiency of 86.7% (η) and a superior fatigue resistance, as well as a superior charge–discharge performance (CD∼1623.14A/cm2, PD∼162.31MW/cm3, WD∼1.30J/cm3, t0.9∼49.6ns) and good thermal stability. Besides, upon NL doping, the FE to RFE phase transition results in a dielectric behavior traversing the relaxation phase (γ∼1.84) accompanied by frequency dispersion, which is beneficial to improve ESP. The superior ESP in this work indicates that BCZT–0.12NL ceramics are a promising candidate used in energy storage capacitors.
1. Introduction
The electronics industry and new-generation information industry technologies such as artificial intelligence, the Internet of Things, and 5G are rapidly developing and spreading. The modernization of information science is driving the rapid development of energy storage materials such as batteries and dielectric capacitors toward higher power density and larger energy storage density.1,2 Due to its advantages of high energy storage density, large power density, extremely fast energy storage speed and miniaturization, environmental protection dielectric ceramic capacitors have attracted extensive attention in pulse power weapons, wind variable propeller system, spacecraft and other high power pulse power systems.3,4,5 Therefore, the development of dielectric ceramic capacitors with excellent energy storage density and efficiency is imminent. It is widely recognized that the evaluation of energy storage performance (ESP) can be achieved by measuring the hysteresis loops of polarization–electric field (P–E), which can be calculated as follows6,7:
In this work, a new RFE Ba0.85Ca0.15[(Zr0.1Ti0.9)(1−x)–(Nb0.5La0.5)x]O3(BCZT-xNL) was successfully prepared by the traditional solid-state preparation method. The double ion (Nb0.5La0.5)4+ replaces the (Zr0.10Ti0.90)4+ ion with oxygen and cation vacancies to compensate for the charge balance. Charge compensation makes BCZT ceramics able to withstand larger external electric fields without being broken down.19,24 As shown in Fig. 1, the introduction of large radius ions inhibits the growth of grains, decreases the size of grains and increases the number of grain boundaries, which helps to improve the BDS of BCZT and thus improve the electrostatic properties of ceramics. The experimental results show that BCZT–0.12NL ceramics have large breakdown field strength (∼410kV/cm), good energy storage density (∼3.1J/cm3) and efficiency (∼86.7%). Consequently, NL doping is a new strategy to improve the ESP of BCZT ceramics.

Fig. 1. Strategy of BCZT to obtain high ESP.
2. Experimental Procedures
Ba0.85Ca0.15[(Zr0.1Ti0.9)1−x–(Nb0.5La0.5)x]O3 (x=0.0, 0.025, 0.06, 0.09, 0.12, 0.15) ceramics were prepared using the traditional solid-state method. Reagent grade raw materials of BaCO3 (99%), CaCO3 (99.99%), ZrO2 (99.99%), TiO2 (99.99%), Nb2O5 (99.99%) and La2O3 (99.5%) with an accurate stoichiometric ratio were mixed by ball milling for 24h in anhydrous ethanol. After drying, the mixtures were pre-sintered for 2h at 1200∘C. Next, the powder added with PVA was ground for 2h and pressed into a 10mm disks at 6MPa, then the cold isostatic pressing technology was used to hold the pressure at 220MPa for 5min. Finally, the disks were sintered at 1400–1450∘C for 4h in air to obtain BCZT–xNL ceramics. The manufactured round flake ceramic samples were ground, polished, and tested after being electrodeposited. The specification of the samples for dielectric test is 8mm in diameter, 0.6mm in thickness, and the upper and lower surfaces are coated with electrodes; the specification of the samples for ferroelectric test is 8mm in diameter, 0.15mm in thickness, and the electrode area is 0.0314cm2.
The crystal structure of ceramic samples was characterized by X-ray diffractometer (MiniFlex 600, Rigaku, Japan) with Cu Kα1 radiation (λ=1.540562Å) in a 2θ range of 20–80∘. The surface micromorphology of ceramics was observed by cold field emission scanning electron microscope (SU-8020). The dielectric properties’ temperature dependence was tested by Dielectric/Impedance spectrometer (TH2827A, Tonghui, TZDM-200-300) in the frequency range from 1kHz to 200kHz. Impedance spectra were obtained by high temperature dielectric impedance measurements (DMS-2000, BALAB). The P–E loops were measured via a ferroelectric analyzer (Model11610E, Radiant, America). The charge–discharge performances of BCZT–xNL ceramics were characterized via a pulse charge–discharge instrument (CFD-001, Shanghai Tongguo Intelligent Technology Co. Ltd., China).
3. Results and Analyses
Figure 2 illustrates the XRD patterns of BCZT–xNL (x=0.0, 0.025, 0.06, 0.09, 0.12, 0.15) ceramics in a 2θ range of 20–80∘. The XRD patterns exhibit a single perovskite phase without secondary phase formation, indicating that (Nb0.5 La0.5)4+ entirely diffuses into BCZT lattice and forms a stable solid solution.25 With the increase of (Nb0.5La0.5)4+ content, the positions of (200) peaks shift slightly toward lower angles. This indicates a distortion of the lattice, which may be attributed to the larger ionic radius of (Nb0.5La0.5)4+ with 6-coordination (0.835Å), compared to (Zr0.10Ti0.90)4+ (0.692Å).26 To further investigate the crystal structure of ceramics, XRD data of BCZT–xNL (x=0, 0.12) ceramics were refined by GSAS (General Structure Analysis System) software, and the data are presented in Table 1. Pure BCZT ceramic has the coexistence of rhombohedral (R) and tetragonal phases (T) for MPB composition near room temperature.27 As shown in Fig. 2(c), pure BCZT ceramics show the structure of Rhombohedral (R3m) and Tetragonal (P4mm), with the weighing fraction accounting for 69.50% and 30.50%. With the introduction of (Nb0.5La0.5)4+, when x=0.12, the phase structure of the ceramic sample becomes a pseudo-cubic phase, as shown in Fig. 2(d). Finally, according to the finishing results, the relative densities of components x=0 and 0.12 are calculated, and the calculation results are shown in Table 1. The calculation formula is as follows :

Fig. 2. (a) XRD patterns of BCZT-xNL ceramics, (b) XRD curves near 38–46∘, (c)–(d) BCZT-xNL (x=0, 0.12) refinement results.

Fig. 3. (a)–(f) Scanning electron micrographs of BCZT-xNL ceramics after thermal corrosion, (g) Weibull distribution of BDS of BCZT-xNL ceramics, and (f) bulk density, average grain size and BDS of BCZT-xNL ceramics.
Sample | x=0 | x=0.12 | |||
---|---|---|---|---|---|
Crystal system | Atoms | Rhombohedral | Tetragonal | Atoms | Cubic |
Space group | R3m | P4mm | Pm3m | ||
Cell parameters | a=b | 4.627Å | 4.024Å | a=b=c | 4.0193Å |
c | 4.003Å | 3.997Å | |||
c∕a | 0.866 | 0.993 | c∕a | 1 | |
α=β=90∘, γ=120∘ | α=β=γ=90∘ | α=β=γ=90∘ | |||
Weigh fraction | 69.50% | 30.50% | 100% | ||
R-factors | Rp | 7.27% | Rp | 6.04% | |
Rwp | 9.72% | Rwp | 8.84% | ||
χ2 | 1.32 | χ2 | 1.07 | ||
Relative density | 93.18% | 95.52% |
As illustrated in Figs. 4(a)–4(f), the changes of dielectric constant and dielectric loss of BCZT-xNL ceramics with temperature were measured at different frequencies (1, 10, 50, 100 and 200kHz). Pure BCZT ceramics show two obvious phase transition peaks at 32∘C and 82.5∘C, namely the phase transition from rhombic to tetragonal (TR−T) and from tetragonal to cubic (TC), as shown in Fig. 4(a).27,35 The ceramic phase structure is altered by the transformation of two-phase transition peaks into one peak after doping with (Nb0.5La0.5)4+. At room temperature, when x=0.12, the paraelectric structure of the pseudo-cubic phase is observed, which is consistent with the refinement results. With the increase of doping amount, dielectric relaxation behavior with frequency dispersion effect becomes more obvious, which is due to the doping of double ions disrupting the long-range ordered structure inside the BCZT ceramics and forming short-range ordered structure polar nanoregion (PNR).36,37,38,39 In Fig. 4(a), the dielectric properties of BCZT–xNL (x=0.0, 0.025, 0.06, 0.09, 0.12, 0.15) ceramic at 1kHz and −160–200∘C were compared. With the increase of x value, the maximum permittivity dropped sharply, and the dielectric peak widened, TC moved to the low temperature region. Numerous studies have shown that dielectric with medium dielectric constant is beneficial to high energy density.40,41,42,43 In addition, dielectric loss also decreases with the decrease of doping and is stable over a wider temperature range, which facilitates the reduction of thermal losses and the possibility of thermal breakdowns.44,45 Diffusion factor γ is an important parameter describing the degree of relaxation of ferroelectric ceramics, which is generally calculated by modified Curie–Weiss law46,47 :

Fig. 4. (a)–(f) Temperature dependence of the dielectric constant and dielectric loss for BCZT-xNL ceramics with x=(a) 0, (b) 0.025, (c) 0.06, (d) 0.09, (e) 0.12, (f) 0.15.

Fig. 5. (a) Temperature dependence of the dielectric constant and dielectric loss for BCZT-xNL ceramics at 1kHz, (b) ln(1∕εr−1∕εm) versus ln(T−Tm) for BCZT-xNL (x=0, 0.025, 0.06, 0.09, 0.12) ceramics.
The microstructure of polycrystalline dielectric ceramics has a significant influence on their electrical properties.50,51 Complex impedance analysis (CIA) is a major means to study the internal electrical properties of ceramics. Figure 6 shows the Nyquist plots as well as the fitted curves for the BCZT–xNL ceramics. The semicircular arcs of all samples decrease with increasing temperature, which indicates a decrease in resistivity.9,52 The complex impedance profile of ideal polycrystalline ceramics shows three semicircular arcs53 and only the high-frequency section arcs contributed by the grains and the low-frequency section arcs contributed by the grain boundaries when the frequency range of the test is not sufficiently low, as shown in Figs. 6(a) and 6(b). In order to obtain the values of Rg and Rgb, two equivalent circuits consisting of a resistor in parallel with a capacitive phase element (CPE) in series are used to fit. Generally speaking, grain boundaries of ceramics have stronger insulation and capacitance than grains, mainly due to the existence of nonstoichiometric distributed oxygen and dangling bonds on the grain boundaries, which act as charge carrier traps and barrier layers.54,55 When the (Nb0.5La0.5)4+ is doped, the two arcs gradually merge into one and grain boundary conduction dominates while containing mixed, ionic and electronic conduction.56 BCZT–0.15NL ceramics show arcs at low frequencies, which can be caused by ionic spatial polarization and oxygen vacancies being blocked at the electrodes.57,58 The grain resistance (Rg) and grain boundary resistance (Rgb) were obtained by fitting the data for all samples within 440–500∘C using Z-View software. In all cases, the fitted data matched well with the test data. The conductivity of the BCZT–xNL ceramics is in accordance with the Arrhenius equation59 :

Fig. 6. BCZT-xNL ceramic Nyquist plot at different temperatures along with the fitting results using (a) x=0.00, (b) x=0.025, (c) x=0.06, (d) x=0.09, (e) x=0.12, (f) x=0.15.

Fig. 7. (a) Grain conductivity and (b) grain boundary conductivity as a function of 1000∕T, the variation of the fitted activation energy with x in the inset figure.
Figure 8(a) illustrates the bipolar P–E circuits of BCZT–xNL ceramics at 100Hz, 90kV/cm. It can be seen that, as expected by the design, the P–E curves become elongated with NL doping, and the BCZT–xNL ceramics transform from typical ferroelectrics to relaxation ferroelectrics, which is in agreement with the dielectric results.10,16 Figure 8(b) shows the unipolar P–E loops of BCZT–xNL ceramic under critical breakdown. It can be preliminarily seen that the P–E loop of BCZT–0.12NL ceramic is thinner and has the largest BDS. The variations of Pm, Pr and BDS with component of x in BCZT ceramics are shown in Fig. 8(c). As the amount of x increases, the Pm and Pr continue to decrease, while the BDS continues to increase. The decrease in Pr is mainly due to the construction of relaxation ferroelectrics and may also be related to the introduction of NL to produce PNRs. The increase in Pr at x=0.15 may be due to the generation of niobium vacancies as a result of over-doping. Niobium vacancies, as negatively charged defects, tend to pair with oxygen vacancies to form defective dipoles, which can be aligned with spontaneous polarization in the presence of an applied electric field.60,61 The removal of the applied electric field resulted in the formation of high residual polarization in BCZT–0.15NL ceramics. The ESP of ceramics are calculated by Eqs. (2) and (3), and the results are shown in Fig. 8(d). It can be seen that BCZT–0.12NL ceramic has the best ESP, reaching an energy storage density (Wrec) of ∼3.11J/cm3, an energy storage efficiency (η) of ∼86.7%, and a breakdown field strength of ∼410kV/cm.

Fig. 8. (a) Bipolar P–E loops of BCZT-xNL ceramics at 90kV/cm, (b) unipolar P–E loops of BCZT-xNL ceramics at breakdown field strength, (c) Pm, Pr, and BDS of BCZT-xNL ceramics, and (d) energy storage density and efficiency computed from the (b) plot.
Considering the practical applications, the frequency stability and anti-fatigue characteristics of ferroelectric properties of BCZT–0.12NL ceramic were investigated at 150kV/cm as shown in Fig. 9. In the range 5–150Hz, the Wrec fluctuation is within 11.4%, and the η fluctuation is within 10.8%. Under the 105 cycles test, the Wrec fluctuation is only 3.3%, and the η fluctuation is only 2.7%. Therefore, BCZT–0.12NL ceramics have good frequency stability and excellent anti-fatigue characteristics. It has practical applications value in pulse power systems and high energy storage capacitors.

Fig. 9. (a), (b) Frequency stability of BCZT–0.12NL ceramic ESP, (c), (d) Fatigue stability of BCZT–0.12NL ceramic ESP.
Power density and current density are important evaluation indexes of ceramic capacitors in practical applications. Figure 10 displays the discharge behavior of 100–200kV/cm of BCZT–0.12NL ceramic measured in a DC electric field. The power density (CD) and current density (PD) of ceramics can be calculated by using the I–t curve under underdamped conditions (load resistance is 0):

Fig. 10. BCZT–0.12NL ceramic (a) underdamped impulse discharge current profile, (b) corresponding CD and PD.
To better represent the storage capacity of the ceramic, the overdamped (nonzero load) discharge current waveform is shown in Figs. 11(a)–11(c). The discharge energy density (WD) of dielectric storage materials is as follows62,63,64 :

Fig. 11. (a)–(c) Overdamped pulse performance for BCZT–0.12NL ceramic, (d)–(f) variable temperature overdamped pulse performance for BCZT–0.12NL ceramic.
4. Conclusions
In this paper, BCZT–xNL ceramics were successfully synthesized by the traditional solid-state method. XRD showed a single pseudo-cubic phase perovskite structure and SEM images showed a compact uniform microstructure. The introduction of the B-site double ion (Nb0.5La0.5)4+ inhibited the grain growth. The small grain size (0.70μm) and more grain boundaries result in a large breakdown field strength (BDS∼410kV/cm). The wide dielectric peak is accompanied by frequency dispersion and strong relaxation behavior (γ∼1.84). The FE to RFE phase transition leads to low Pr, thin hysteresis loop, high energy storage efficiency (∼86.7%) and high energy storage density (3.1J/cm3). At the same time, it has good frequency stability and excellent fatigue resistance. In addition, the maximum CD (∼1623.14A/cm2), PD (∼162.31MW/cm3), WD (∼1.301J/cm3) and extremely low discharge speed t0.9 (∼49.6ns) are suitable for applications in pulse power capacitors.
Acknowledgment
This work was supported by the National Science Foundation of China (NSFC) (Grant Nos. 52272119 and 52202143). The authors would also like to thank the Natural Science Basic Research Plan in the Shaanxi Province of China (Grant No. 2022JQ338), Young Talent Fund of University Association for Science and Technology in Shaanxi, China (20230415), the Fundamental Research Funds for the Central Universities (Grant No. GK202401009), Key Research and Development Program of Shaanxi Provincial Science and Technology Department (Grant No. 2023-YBGY-162) and the Fundamental Innovation Project in School of Materials Science and Engineering (SNNU).
ORCID
Zhanhui Peng https://orcid.org/0000-0003-1000-258X
Xiaolian Chao https://orcid.org/0000-0003-2957-1192
Zupei Yang https://orcid.org/0000-0001-5096-2134