Processing math: 100%
World Scientific
Skip main navigation

Cookies Notification

We use cookies on this site to enhance your user experience. By continuing to browse the site, you consent to the use of our cookies. Learn More
×

System Upgrade on Tue, May 28th, 2024 at 2am (EDT)

Existing users will be able to log into the site and access content. However, E-commerce and registration of new users may not be available for up to 12 hours.
For online purchase, please visit us again. Contact us at customercare@wspc.com for any enquiries.

High energy storage properties of (Nb0.5La0.5)4+ complex-ion modified (Ba0.85Ca0.15)(Zr0.10Ti0.90)O3 ceramics

    https://doi.org/10.1142/S2010135X24400101Cited by:4 (Source: Crossref)
    This article is part of the issue:

    Abstract

    In this work, a traditional solid-state method was adopted to prepare dense Ba0.85Ca0.15[(Zr0.1Ti0.9)1x–(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 (CD1623.14A/cm2, PD162.31MW/cm3, WD1.30J/cm3, t0.949.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 (PE), which can be calculated as follows6,7:

    Wtot=PmPrEdP,(1)
    Wrec=PmPrEdP,(2)
    η=WrecWtot×100%,(3)
    where the maximum polarization intensity (Pm), remnant polarization (Pr) and electric field intensity (E) directly determine the value of Wrec. However, there is a contradiction between achieving large ΔP (PmPr) and large breakdown field strength (BDS).8,9,10 Therefore, it’s necessary to get a balance between ΔP and BDS to optimize the ESP of dielectric ceramics.11 Relaxor ferroelectrics (RFEs) own large Pm and small Pr accompanied with moderate BDS. It has been shown that the BDS of RFEs can be effectively increased by optimizing the doping strategy of the component design.12,13,14 Recently, researches of the ESP of BCZT-based ceramics have been widely concerned.10,15,16,17 Zhao et al. successfully constructed a “relaxation ferroelectric” by introducing (Al0.5Nb0.5)4+ at the B site of calcite structure, realizing improved Wrec of 0.70J/cm3 in BNBT–xAN ceramics and 1.41J/cm3 in BNKT–xAN ceramics.18,19 Zhou et al. introduced the (Ni13Nb23)4+ with large ionic radius into BCZT and controlled the grain size to submicron scale to increase BDS.20 Chu et al. reported that Nb2O5 as a good firing aid can promote domain wall movement, reduce porosity and increase ceramic densities.21 Kumar et al. investigated Nb2O5-doped Bi0.5(Na0.5K0.5)0.5TiO3 ceramics and reported that the ceramic with 0.4%Nb exhibits relaxation behavior with respect to the degree of diffusivity of the phase transition.22 Similarly, the strong relaxation behavior was observed in the Ba0.90Ca0.10Zr0.10Ti0.90O3 system (γ1.85).23 Further, compared to single-ion doping, double-ion substitution is more likely to achieve high Wrec and η at low to medium field strengths. The local random field induced by charge and ion radius imbalance breaks the remotely ordered ferroelectric state to produce PNRs, which makes it easier to construct relaxing ferroelectrics. Under the action of electric field, nanodomains are more likely to invert than micro-domains, which leads to a decrease in Pr and an increase in η. Therefore, the construction of RFEs by B-site double ion substitution is a feasible and effective strategy to improve the ESP of BCZT-based ceramic.

    In this work, a new RFE Ba0.85Ca0.15[(Zr0.1Ti0.9)(1x)–(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.

    Fig. 1. Strategy of BCZT to obtain high ESP.

    2. Experimental Procedures

    Ba0.85Ca0.15[(Zr0.1Ti0.9)1x–(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 1200C. 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–1450C 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 PE 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 :

    ρb=m1m1m2(ρ0ρl)+ρl,(4)
    where ρl is the density of air (0.0012g/cm3), m1 is the mass of the sample in air, m2 is the wet weight of the sample in distilled water, the ρ0 is the density of distilled water at the test temperature. The theoretical density (ρth) of ceramic samples is calculated as follows :
    ρth=ZMVNA,(5)
    where M is the molar mass of the compound; V is the cell volume of the compound; Z is the number of molecules in the unit cell; NA is Avogadro’s constant. According to the bulk density (ρb) and theoretical density (ρth) of the ceramic sample, its relative density (ρr) can be calculated :
    ρr=ρbρth×100%.(6)
    Figure 3 displays the SEM micrographs of NL co-doped BCZT ceramic samples after thermal corrosion recorded. The grain of all ceramic samples was dense and homogeneous micromorphology. Large densities are a key factor in obtaining high breakdown field strengths in ceramics, and large pores and many cracks lead to large leakage currents and thus low breakdown field strengths. Cold isostatic pressure technology can make ceramic samples in all directions by 220Mpa uniform pressure. The sintering, before the use of physical pressure to tighten the ceramic powder particles, can effectively reduce the sintering of large pores that cannot be a discharged phenomenon, thus greatly reducing the pores and cracks. Relative density can reflect the porosity and compactness inside the ceramic. The closer the relative density is to 1, the lower the porosity and the better the compactness. According to Table 1, BCZT–0.12NL has the best densification. Besides, the grain size distribution is recorded by Nano Measurer software, as shown in Figs. 3(a)–3(f). Compared with pure BCZT (22.95μm), the grain size of doped ceramic sample is significantly reduced, which is caused by the large ionic radius of the (Nb0.5La0.5)4+ B-site complex-ion.15,16,28 The low diffusivity of La3+ as a rare earth ion retards the mass transfer of particles during ceramic sintering.29,30 The reduced grain size and high density would contribute to the enhanced dielectric breakdown field strength, which is a key factor benefiting the high ESP. The results of the Weibull distribution show the change tendency of BDS, as Fig. 3(g). Weibull plots, which need eight sets of sample breakdown data are used to interpret BDS data.31,32 the average value of BDS can be expressed by the following formula19,33:
    Xi=ln(Ei),(7)
    Yi=ln(ln(1i1+n)).(8)
    where Xi and Yi correspond to the x- and y-axes of the Weibull distribution, Ei is BDS of each specific at the ith time, i is the rank of specimens, and n is the sum of specimens. The average BDS determined by the fitting lines increases with the (Nb0.5La0.5)4+ addition. This may be due to the complex-ion substitution in B-site, which leads to charge imbalance, and oxygen vacancy compensation charge balance, which leads to greater breakdown field strength.19,20,24 In addition, according to the work of Tunkasiri, the effect of grain size on BDS can be expressed by EGα, where E is the field strength, G is the grain size, and α is a constant.34 As a result, the large compactness, small grains and charge compensation enable the doped BCZT ceramics to obtain a higher breakdown electric field. The results show that BCZT–0.12NL ceramics have the best density ρb up to 5.72g/cm3, BDS up to 400kV/cm, and an average grain size of 0.7μm.

    Fig. 2.

    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.

    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.

    Table 1. Refinement results of BCZT-xNL (x=0, 0.12) ceramic XRD data.

    Samplex=0x=0.12
    Crystal systemAtomsRhombohedralTetragonalAtomsCubic
    Space groupR3mP4mmPm3m
    Cell parametersa=b4.627Å4.024Åa=b=c4.0193Å
    c4.003Å3.997Å
    ca0.8660.993ca1
    α=β=90, γ=120α=β=γ=90α=β=γ=90
    Weigh fraction69.50%30.50%100%
    R-factorsRp7.27%Rp6.04%
    Rwp9.72%Rwp8.84%
    χ21.32χ21.07
    Relative density93.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 32C and 82.5C, namely the phase transition from rhombic to tetragonal (TRT) 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–200C 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 :

    1ε1εr=(TTm)γC,(9)
    where C represents the modified Curius coefficient, γ is the diffusion factor, taking a value between 1 and 2, the closer to 2 the stronger the relaxation. Figure 5(b) shows the γ value obtained from the calculation of the permittivity at 1kHz. With the addition of the doping, the γ factor increases and remains above 1.8. It can be seen that the appropriate addition of (Nb0.5La0.5)4+ B-site complex-ion can effectively lead to the generation of PNRs, which can increase the relaxation and ESP.15,48,49

    Fig. 4.

    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.

    Fig. 5. (a) Temperature dependence of the dielectric constant and dielectric loss for BCZT-xNL ceramics at 1kHz, (b) ln(1εr1εm) versus ln(TTm) 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–500C 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 :

    σ=σ0exp(EaKT),(10)
    where σ0 is a constant, σ stands for the conductivity value, Ea is the activation energy, K is the Boltzmann’s constant and T is the Kelvin temperature. Figures 7(a) and 7(b) show the variation curves of grain conductivity and grain boundary conductivity with 1000T, respectively. It can be seen that Eg (1.99eV) and Egb (2.96eV) are maximum when x=0.12. Therefore, the doping of ions leads to grain reduction and increase in grain boundaries. Carriers need to penetrate past the grain boundaries for transport within the ceramic, which leads to an increase in activation energy, which helps increase BDS to some extent.

    Fig. 6.

    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.

    Fig. 7. (a) Grain conductivity and (b) grain boundary conductivity as a function of 1000T, the variation of the fitted activation energy with x in the inset figure.

    Figure 8(a) illustrates the bipolar PE circuits of BCZT–xNL ceramics at 100Hz, 90kV/cm. It can be seen that, as expected by the design, the PE 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 PE loops of BCZT–xNL ceramic under critical breakdown. It can be preliminarily seen that the PE 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.

    Fig. 8. (a) Bipolar PE loops of BCZT-xNL ceramics at 90kV/cm, (b) unipolar PE 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.

    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 It curve under underdamped conditions (load resistance is 0):

    CD=ImaxS,(11)
    PD=EImax2S,(12)
    where Imax is the maximum current in the It curve, S is the electrode area of the test sample, and E is the field voltage. As illustrated in Fig. 10(a), the current underdamping increases with the electric field increase. Figure 10(b) is calculated from Fig. 10(a). When the external electric field is 200kV/cm, the Imax (50.97A), CD (1623.14A/cm2) and PD (162.31MW/cm3) of BCZT–0.12NL ceramics reach the maximum value. The high power density indicates that BCZT–0.12NL ceramics have great application potential in pulsed power systems.

    Fig. 10.

    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 :

    WD=Ri(t)2dtV,(13)
    where R (300Ω) is the load resistance, i is the current, t is the discharge time and V is the single volume of the test sample. t0.9 is defined as the time required when WD reaches 90%, and is generally used to measure the speed of discharge.65,66,67 Similarly, with the E increases, WD increases from 0.42J/cm3 to 1.30J/cm3, and t0.9 eventually remains at 49.6ns. Finally, the temperature stability of the overdamped charge–discharge performance of BCZT–0.12NL ceramic (30–190C) was explored, as shown in Figs. 11(d)–11(f). The results show that the WD and t0.9 of BCZT–0.12NL ceramic overdamping fluctuate moderately in the temperature range of 30–190C with good temperature stability.

    Fig. 11.

    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 (BDS410kV/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