Open Access Article
This Open Access Article is licensed under a Creative Commons Attribution-Non Commercial 3.0 Unported Licence

Effect of Na2+ and Bi3+ substitution in lead-free BZT-Pr multifunctional ceramics for energy storage and thermal energy conversion

I. Zouari*ab, A. Dahrib, Z. Sassic, L. Seveyratc, N. Abdelmoulab and H. Khemakhemb
aLaboratory of Nanomaterials and Systems for Renewable Energies (LaNSER) LR15CRTEn05, Research and Technology Center of Energy (CRTEn), Techno-Park Borj Cedria, Bp 95, Hammam-Lif, Tunis, 2050, Tunisia. E-mail: ilhem.zouari@crten.mesrs.tn
bLaboratory of Multifunctional Materials and Applications (LaMMA), LR16ES18, Faculty of Sciences of Sfax, University of Sfax, BP 1171, 3000 Sfax, Tunisia
cUniv Lyon, INSA-LYON, LGEF, EA682, Villeurbanne, F-69621, France

Received 10th January 2025 , Accepted 23rd March 2025

First published on 28th March 2025


Abstract

For this study, Ba1−yxPr2y/3y/3(Na0.5Bi0.5)xTi0.9Zr0.1O3 (y = 0.5%; x = 0.05, 0.1 and 0.15) (abbreviated as BZT-Pr–xNBT) ceramics were prepared by a solid state reaction method. The samples were structurally characterized by X-ray diffraction (XRD) and the software Fullprof, which showed the atomic positions, bond distances, and bond angles of these materials in the tetragonal structure. The crystallite size, determined by both the Williamson–Hall and Debye–Scherrer methods, increases with increasing Na0.5Bi0.5TiO3 (abbreviated as NBT) content in Ba1−yPr2y/3y/3Ti0.9Zr0.1O3 (y = 0.5%) (abbreviated as BZT-Pr) ceramic. The ultraviolet-visible (UV-Vis) diffuse reflectance spectrum indicated that the BZT-Pr–xNBT ceramics exhibited direct band gap values in the range of 2.90–2.78 eV. Moreover, this study investigates the ferroelectric, energy storage, and pyroelectric energy harvesting properties of BZT-Pr–xNBT ceramics with varying NBT contents (x = 0.05, 0.10 and 0.15). The results reveal that NBT substitution significantly influences the ferroelectric characteristics and enhances recoverable energy storage density (Wrec) and energy storage efficiency (η). Temperature-dependent scaling relations for the coercive field (EC) and remnant polarization (Pr) were derived for all samples. For BZT-Pr–0.15NBT ceramic, the scaling relation for EC and Pr (ECT−0.089, PrT−0.233) highlights improved thermal stability in this composition. Pyroelectric energy harvesting under the Olsen cycle identified BZT-Pr–0.05NBT as the most efficient composition. These findings emphasize the crucial impact of NBT substitution in optimizing the ferroelectric and pyroelectric properties of BZT-Pr–xNBT ceramics.


1. Introduction

Ferroelectric materials have gained significant technological and commercial importance over the past few decades, attracting growing interest due to their remarkable properties.1–6 These include hysteresis, which enables their use in non-volatile memory;7,8 high dielectric constants,9,10 making them ideal for capacitors; and excellent high-temperature capacitance stability, crucial for automotive electronic systems such as multilayer ceramic capacitors.11,12 Additionally, their narrow energy band gaps make them suitable for photocatalytic and optoelectronic applications,13,14 while their high pyroelectric coefficients are valuable for infrared detection.15,16 Furthermore, their exceptional energy conversion capabilities play a critical role in powering micro-electro-mechanical systems, demonstrating their versatility in a wide range of applications.17–20 In all these fields and applications, lead-free materials have garnered increased attention in recent bibliographic studies due to significant environmental and health concerns.21–25 Specifically, in the field of thermal energy conversion (pyroelectric energy harvesting) and energy storage, BCT–BZT ceramics,26,27 Sr–BCTZ,28–30 Sn–BSTZ ceramics,31 and NBT–6BT–Sr compositions32 are promising candidates.

In our previous works, we have already studied some physical properties of the Ba1−yxPr2y/3y/3(Na0.5Bi0.5)xTi0.9Zr0.1O3 (y = 0.5%; x = 0.05, 0.1 and 0.15) (abbreviated as BZT-Pr–xNBT) system.33 We have determined the lattice parameters of the compositions in their tetragonal structure, and we have observed an improved dielectric temperature stability with broadened curves due to a core–shell grain structure in these ceramics. Specifically, the dielectric study revealed that the BZT-Pr–0.15NBT composition meets the EIA X7R specification. In addition, a vibrational and photoluminescence study was conducted on these compositions. The analysis of hysteresis cycles, which demonstrate the ferroelectric properties of these materials, was limited to determining the remanent polarization and coercive field values for each composition.

In the present study, we will make a complementary study of the physical properties of this BZT-Pr–xNBT (x = 0.05, 0.1 and 0.15) system. Indeed, using X-ray diffraction (XRD) data, the atomic positions, bond distances, bond angles, and crystallite size have been determined. In addition, using UV-Vis analysis spectra, the optical band gaps have been calculated. Moreover, based on hysteresis loops data, the energy storage properties of these materials have been studied. Moreover, we analyze the temperature-dependent scaling relations for the coercive field (EC) and remnant polarization (Pr) for all samples to study thermal stability. In addition, we undertake a more in-depth investigation of these cycles by examining the energy conversion process (pyroelectric harvesting energy) based on the Olsen cycle.

2. Experimental

Ceramics were prepared using the formula Ba1−yxPr2y/3y/3(Na0.5Bi0.5)xTi0.9Zr0.1O3 (y = 0.5%; x = 0.05, 0.1 and 0.15) (abbreviated as BZT-Pr–xNBT) by the conventional solid-state reaction method. Stoichiometric amounts of BaCO3, ZrO2, TiO2, Pr6O11, Bi2O3 and Na2O3 powders were weighed, mixed, calcined at 1150 °C for 12 h, and then pressed into discs which were sintered at 1300 °C for 4 h. X-ray diffraction (XRD) patterns were recorded at room temperature (RT) on an XPert Pro MPD diffractometer using CuKα1 and CuKα2 radiation (λ1 = 1.540598 Å, λ2 = 1.544426 Å). The results were processed by a Rietveld refinement using the software Fullprof in order to determinate the structure, cell parameters, atomic positions, bond angles and bond distances. The microstructure of the samples was examined using scanning electron microscopy (SEM) (Hitachi FlexSEM 1000) in secondary electrons (SE) mode at an accelerating voltage of 5 kV. UV-Vis diffuse reflectance spectra were obtained using a PerkinElmer Lambda 950 model. The optical band gaps of different compositions were estimated by extrapolation from Tauc plot of the Kubelka–Munk function.34 Ferroelectric (PE) loops were determined at different temperature (from 30 °C to 140 °C) by applying different electric field (1 kV mm−1, 2 kV mm−1 and 3 kV mm−1) at a frequency of 1 Hz. The current and electric field were measured during the application of cyclic electric fields, using a Keithley 428 current amplifier and a TREK Model 20/20C high-voltage amplifier.

3. Results and discussion

3.1 Structural and morphological properties

Using the software “Fullprof”,35 Rietveld refinements of the XRD patterns, at room temperature, were performed in the tetragonal structure (P4mm) for all compositions. The cell-parameters refinements, agreement factors, atomic positions, bond distances, and bond angles for the BZT-Pr–xNBT (x = 0.05, 0.1 and 0.15) ceramics are shown in Table 1. The refined XRD pattern of x = 0.1 composition, as an example, is shown in Fig. 1.
Table 1 Crystallographic data and agreement factors for BZT-Pr–xNBT compounds
Cell parameters Compositions Ba1−yPr2y/3y/3Ti0.9Zr0.1O3 (y = 0.5%)
x = 0.05 x = 0.1 x = 0.15
a = b (Å) 4.021(0) 4.018(3) 4.017(1)
c (Å) 4.029(4) 4.024(1) 4.025(7)
V3) 65.147(8) 64.975(1) 64.96(1)
Rp 7.38 5.70 5.69
Rwp 9.80 7.52 7.36
Rexp 9.05 6.76 6.11
χ2 1.17 1.24 1.45
[thin space (1/6-em)]
Atomic positions (x; y; z)
Ba/Na/Bi/Pr 0 0 0
0 0 0
−0.059(9) 0.001(5) 0.020(5)
Ti/Zr 0.5 0.5 0.5
0.5 0.5 0.5
0.4632(1) 0.526(8) 0.541(9)
O1 0.5 0.5 0.5
0 0 0
0.4179(6) 0.473(3) 0.489(0)
O2 0.5 0.5 0.5
0.5 0.5 0.5
−0.042(0) 0.003(5) −0.003(7)
[thin space (1/6-em)]
Bond angles (°)
α [O1–Ti–O2] 84.818(3) 83.88(7) 83.953(4)
[thin space (1/6-em)]
Bond distances (Å)
O1–Ti 2.0187(5) 2.020(6) 2.019(7)
Ti–O2 2.0355(9) 2.105(6) 2.196(5)



image file: d5ra00245a-f1.tif
Fig. 1 Refined X-ray diffraction pattern of BZT-Pr–xNBT (x = 0.1) ceramic.

In addition, the crystallite size, determined by both the Williamson–Hall36 and Debye–Scherrer37 methods, are shown in Table 2. There is an increase of particle size, along with an increase of crystallinity observed by XRD, when the substitution of barium (Ba) by sodium (Na) and bismuth (Bi) increases from x = 0.05 to x = 0.15.

Table 2 Crystallite size determined by Debye–Scherrer formula (DSC) and Williamson–Hall formula (DW–H) for BZT-Pr–xNBT compounds
Compositions DSC (nm) DW–H (nm)
x = 0.05 26.52 48.29
x = 0.1 29.76 55.88
x = 0.15 35.27 63.57


As shown in Fig. 2a–c, the surface morphology of the compounds, examined using scanning electron microscopy (SEM) at a magnification of 2000×, reveals a relatively homogeneous microstructure with some porosity. There is a few small voids, suggesting the formation of highly densified homogeneous samples. Additionally, the analysis of grains on fractured fragments at a magnification of 4000×, as shown in Fig. 2d–f, indicates an average grain size between 0.5 μm and a few μm. We observe a slight increase in grain size with the NBT content, which is in good agreement with the increase in crystallite size as the NBT content increases.


image file: d5ra00245a-f2.tif
Fig. 2 Surface morphology and grain size analysis of BZT-Pr–xNBT samples: (a–c) surface morphology for x = 0.05, 0.10, and 0.15; (d–f) grain observation on fractured samples for x = 0.05, 0.10, and 0.15.

3.2 Optical properties

For estimating the optical band gap from the ultraviolet-visible (UV-Vis) diffuse reflectance spectrum, the plot of (F(R))2 vs. () is shown in Fig. 3 for BZT-Pr–xNBT (x = 0.05) sample as an example. The values of the direct band gap obtained for the BZT-Pr–xNBT (x = 0.05, 0.1 and 0.15) ceramics are summarized in Table 3. These values are lower than those found in BaTiO3 (BT) ceramics38 and BZT ceramics.39,40 In addition, with the increase of xNBT content, the value of the gap energy (Eg) decreased from 2.90 eV for x = 0.05 to 2.78 eV for x = 0.15. Indeed, the substitution of Ba2+ with Na2+ and Bi3+ can promote lower energy levels, allowing easier electronic transitions at lower energies.
image file: d5ra00245a-f3.tif
Fig. 3 UV-Vis spectrum for BZT-Pr–xNBT ceramics: (a) x = 0.05, (b) x = 0.10, and (c) x = 0.15.
Table 3 Gap energies (Eg) of BZT-Pr–xNBT samples
Compositions Eg (eV)
x = 0.05 2.90
x = 0.1 2.81
x = 0.15 2.78


3.3 Ferroelectric properties

Fig. 4a–c depict the PE loops for BZT-Pr–xNBT ceramics with NBT contents of 5%, 10%, and 15%, respectively, measured under various electric fields at 30 °C and a frequency of 1 Hz. The results reveal that with increasing electric field, minor unsaturated PE loops evolve into fully saturated loops. This transformation, which occurs above a certain threshold electric field, is a consistent phenomenon within finite frequency and temperature ranges.41 The behavior is attributed to increases in the coercive field (EC), maximum polarization (Pmax), and remnant polarization (Pr), driven by domain growth and alignment under the applied electric field.
image file: d5ra00245a-f4.tif
Fig. 4 Polarization versus electric field (PE) hysteresis loops for BZT-Pr–xNBT samples at constant temperature (30 °C) for (a) x = 0.05, (b) x = 0.10, and (c) x = 0.15, measured under varying electric field magnitudes.

The data further indicate that increasing the NBT content results in a transition of the PE loops from a wide profile, characterized by a high EC at 5% NBT, to a slimmer profile. This suggests that NBT substitution is a viable strategy for optimizing materials for dielectric energy storage applications. The BZT-Pr–0.15NBT composition demonstrates a slim PE loop with a nearly linear profile, indicative of high energy storage efficiency.42 However, Pmax decreases with increasing NBT content, a behavior attributed to the core–shell structure of the ceramics. The limited diffusivity of NBT into the BZT-Pr lattice, caused by chemical inhomogeneity, hinders the uniform alignment of domain walls and reduces the material's ferroelectric properties.33,43 Previous studies have linked this reduction in ferroelectric properties to material losses resulting from local compositional fluctuations and random charge effects, which disrupt tetragonal symmetry and weaken ferroelectric order.44

Fig. 5a–c illustrate the temperature dependence of PE hysteresis loops for BZT-Pr–xNBT ceramics (x = 0.05, 0.10, and 0.15), measured at 1 Hz under a constant electric field of 30 kV cm−1. The hysteresis loops, recorded from room temperature to 140 °C, progressively become slimmer with rising temperature, tending toward a linear profile. This trend is confirmed by the ferroelectric parameters (Pmax, Pr, (ΔP = PmaxPr), and EC) extracted under the same conditions.


image file: d5ra00245a-f5.tif
Fig. 5 Polarization versus electric field (PE) hysteresis loops for BZT-Pr–xNBT samples as a function of temperature under a constant electric field of 30 kV cm−1 for (a) x = 0.05, (b) x = 0.10 and (c) x = 0.15.

Fig. 6a–c show the thermal variation of Pmax, Pr, and ΔP for x = 0.05, 0.10, and 0.15, respectively, while Fig. 6d compares EC for all compositions. The results demonstrate a consistent decrease in Pmax, Pr, and EC with increasing temperature which is attributed to the degradation of ferroelectric performance above the Curie temperature (TC) during the ferroelectric-to-paraelectric phase transition. This decline is linked to the breakdown of long-range domain order due to enhanced thermal fluctuations at elevated temperatures.45


image file: d5ra00245a-f6.tif
Fig. 6 Thermal variation of Pmax, Pr and ΔP = (PmaxPr), measured under an electric field of 30 kV cm−1 for (BZT-Pr–xNBT) samples with (a) x = 0.05, (b) x = 0.10 and (c) x = 0.15. In addition, panel (d) presents a comparison of the coercive field (EC).

Moreover, the comparison of EC values highlights the impact of NBT substitution on ferroelectric properties. The composition with the lowest NBT content (x = 0.05) exhibits the highest EC, reflecting a robust ferroelectric nature. Conversely, the composition with the highest NBT content (x = 0.15) shows the lowest EC, indicative of a soft ferroelectric behavior which requires minimal electric field for domain switching.

3.4 Energy storage properties

Ferroelectric materials, renowned for their reversible spontaneous polarization, are highly efficient in energy storage and release, rendering them ideal candidates for energy storage applications.46 Among these, relaxor ferroelectrics exhibit particularly high energy densities, making them increasingly prominent in advanced capacitors designed for high energy requirements and rapid charge–discharge cycles.47 Building on this context, and considering the enhanced ferroelectric properties observed with NBT substitution, the subsequent section focuses on assessing the energy storage capabilities of BZT-Pr–xNBT ceramics with varying NBT content. This evaluation includes key parameters such as recoverable energy storage density (Wrec), total energy storage density (Wtot), and energy storage efficiency (η), derived from the analysis of their ferroelectric hysteresis loops.

The recoverable energy density, representing the energy released during the discharge process, can be evaluated from PE hysteresis loop. This is determined by integrating the area between the discharge curve (the upper branch of the PE loop) and the polarization axis, highlighted in green in Fig. 7a. Similarly, the total energy density, representing the energy stored during the charging process, corresponds to the area between the charge curve (the lower branch of the PE loop) and the polarization axis. The relationship between these quantities is that Wtot is the sum of Wrec and the energy loss density (Wloss), with Wloss being the energy dissipated during the charging and discharging process. Wloss is calculated by integrating the area enclosed by the hysteresis loop, marked in orange in Fig. 7a.


image file: d5ra00245a-f7.tif
Fig. 7 (a) Schematic diagram illustrating the energy density of the BZT-Pr–0.15NBT sample during charging and discharging. Thermal variation of energy storage parameters measured under an electric field of 30 kV cm−1 for (BZT-Pr–xNBT) samples with (b) x = 0.05, (c) x = 0.10 and (d) x = 0.15.

These parameters are quantified using the following equations derived from the PE loops, where P represents polarization, Pmax is the maximum polarization, Pr is the remnant polarization, and E is the applied electric field:48

 
image file: d5ra00245a-t1.tif(1)
 
image file: d5ra00245a-t2.tif(2)

For practical energy storage applications, energy storage efficiency (η) is a critical metric, defined as the ratio of Wrec to Wtot. This relationship is expressed as:

 
image file: d5ra00245a-t3.tif(3)

Higher efficiency (η) indicates minimal energy loss during charge and discharge cycles, which is crucial for reducing heat generation and temperature rise within dielectric materials, thereby preserving the reliability and lifespan of capacitors.49 Optimizing Wrec and η requires a larger ΔP and high dielectric breakdown strength, as these factors enhance energy recoverability and efficiency.50

The temperature-dependent energy storage parameters (Wrec, Wloss, and η) for BZT-Pr–xNBT ceramics with NBT contents of 0.05, 0.10, and 0.15, measured under a constant electric field of 30 kV cm−1 over a temperature range of 30 °C to 140 °C, are presented in Fig. 7b–d. At room temperature (30 °C), the BZT-Pr–5% NBT sample exhibits a coercive field (EC) of 10.068 kV cm−1, a recoverable energy density (Wrec) of 37.23 mJ cm−3, and an energy storage efficiency (η) of 26.18%. Increasing the NBT content to 15% significantly reduces the coercive field to 3.28 kV cm−1, while improving the ΔP value, resulting in a recoverable energy density of 65.86 mJ cm−3 and an efficiency of 69.25%. This enhancement in energy storage properties is attributed to the effects of NBT substitution, which modifies the ferroelectric characteristics of the material. This can be explained by the increase in grain size with increasing NBT content, as previously observed in the SEM analysis. In fact, a larger grain size can enhance energy storage properties by reducing the number of grain boundaries, which can trap charges.

As the temperature increases, the PE loops become slimmer, reducing energy loss during charge–discharge cycles and yielding higher efficiency values. The maximum energy storage efficiency (η) at 140 °C is observed to be 39.66%, 63.05%, and 70.95% for NBT contents of 0.05, 0.10, and 0.15, respectively. These results clearly demonstrate that incorporating NBT significantly enhances the energy storage performance of BZT-Pr–xNBT ceramics, particularly at higher substitution levels. The slim-loop behavior of the BZT-Pr–15% NBT sample contributes to its relatively improved energy storage properties, most notably its thermal stability. This composition maintains an efficiency of approximately 70% over a wide temperature range, highlighting its potential for practical applications requiring consistent performance across wide temperature span.

3.5 The scaling behavior of hysteresis parameters

A comprehensive analysis of the scaling behavior of hysteresis parameters is crucial for understanding the interplay between temperature, electric field, and ferroelectric switching dynamics in BZT-Pr–xNBT. The coercive field and remnant polarization were evaluated as functions of temperature. Logarithmic plots of EC versus ln(T) at various electric fields, shown in Fig. 8a–c, and logarithmic plots of Pr versus ln(T), presented in Fig. 9a–c, revealed a distinct dependence on the NBT content. Using the least squares linear fitting method, the relationship between EC and T can be expressed as follows:28,31
 
ln[thin space (1/6-em)]EC = α[thin space (1/6-em)]ln[thin space (1/6-em)]T + YEC (4)
Here, YEC represents the y-intercept, and α is the slope, both dimensionless constants. The intercept YEC becomes significant as T approaches absolute zero (T = 0). At T = 0, ferroelectric domains freeze, leading to behaviors distinct from those at elevated temperatures. While the impact of EC at T = 0 is often disregarded in practical studies, YEC is retained for maintaining a linear fit. By rewriting eqn (4) in exponential form, the scaling relation can be expressed as:
 
ECTα (5)

image file: d5ra00245a-f8.tif
Fig. 8 Logarithmic plots of coercive field (EC) versus temperature (T) under varying electric fields for BZT-Pr–xNBT ceramics with (a) x = 0.05, (b) x = 0.10 and (c) x = 0.15 and the insets show the slope variations with electric field.

image file: d5ra00245a-f9.tif
Fig. 9 Logarithmic plots of remnant polarization (Pr) versus temperature (T) under varying electric fields for BZT-Pr–xNBT ceramics with (a) x = 0.05, (b) x = 0.10 and (c) x = 0.15 and the insets show the slope variations with electric field.

Similarly, the temperature dependence of the remnant polarization Pr can be described as:

 
ln[thin space (1/6-em)]Pr = γ[thin space (1/6-em)]ln[thin space (1/6-em)]T + YPr (6)
where γ represents the slope and YPr denotes the y-intercept. This equation can also be transformed into its exponential form:
 
PrTγ (7)

For x = 0.05, the scaling relations were determined to be ECT−0.243 and PrT−0.474, indicating a relatively high temperature sensitivity for both parameters. As the NBT content increased to x = 0.1, the exponents reduced to ECT−0.190 and PrT−0.385, reflecting a moderate decrease in temperature dependence. At x = 0.15, the exponents further decreased to ECT−0.089 and PrT−0.233, indicating improved stability of both the coercive field and remnant polarization with temperature.

The evolution of the magnitudes of α and γ with increasing NBT content reflects the variation in domain dynamics with NBT substitution. At low NBT content (x = 0.05), the higher temperature sensitivity, as indicated by larger absolute values of α and γ, suggests a system where domain wall motion is more strongly influenced by thermal fluctuations, leading to increased mobility and switching activity. Conversely, at higher NBT content (x = 0.15), the reduced sensitivity implies a more stabilized domain structure, likely due to stronger domain wall pinning or reduced dipole rotation caused by NBT substitution.28,31,51

The results imply that BZT-Pr–xNBT with higher NBT composition (x = 0.15) is particularly suitable for environments requiring thermal stability. Furthermore, these findings align with energy storage investigations, underscoring the material's potential in practical applications where temperature resilience is critical.

3.6 Pyroelectric energy harvesting properties

Ferroelectric materials exhibit a characteristic hysteresis in their polarization–electric field behavior, which represents the energy required to polarize and depolarize the material. Typically, we observe this hysteresis in a counterclockwise loop, which enables the conversion of electrical energy into heat. However, by switching the cycle's orientation to a clockwise direction, this technique may be used for thermal energy harvesting. The Olsen cycle is designed to convert thermal energy into electricity. This process involves two isothermal and two isoelectric field stages.52,53

The Olsen cycle consists of four key phases31,32 as depicted in Fig. 10:


image file: d5ra00245a-f10.tif
Fig. 10 Schematic representation of electrical energy conversion for BZT-Pr–0.10NBT sample as an example utilizing the Olsen cycle.

In this cycle, the process begins by charging the material at a constant lower temperature (TL) as the electric field increases from EL to EH, which results in an increase in polarization (phase 1–2). In the second stage (phase 2–3), heat (QS) is supplied to the material, raising its temperature from TL to a higher value (TH) while maintaining the electric field at EH. As a result of this heat, the material depolarizes by producing a significant current that can be harvested. The third stage (phase 3–4) involves reducing the electric field from EH to EL while maintaining the material at TH, which results in a weaker depolarization current. Finally, in the fourth stage (phase 4–1), heat is extracted at a constant electric field (EL) which brings the material to its initial state at TL and consequently completing the cycle.

The energy density produced in the cycle is determined by the area enclosed by the PE loop, which can be quantified using the following equation:28,31,32

 
ND = ∮EdP (8)
where E is the applied electric field, and dP is the change in polarization.

The thermal energy conversion potential of BZT-Pr–xNBT ceramics was investigated using the Olsen cycle under varying electric fields and temperature ranges. For compositions with NBT contents of 0.05, 0.10, and 0.15, the cycle was conducted with EL = 0 kV cm−1 and TL = 30 °C. As illustrated in Fig. 11a–d, energy density increases with greater electric field spans (EHEL) and higher operating temperatures (TH). The optimal composition, BZT-Pr–0.05NBT, achieved a maximum energy conversion of 80.77 kJ m−3 when operated between 30–140 °C and 0–30 kV cm−1. In comparison, BZT-Pr–0.10NBT and BZT-Pr–0.15NBT exhibited lower energy densities of 64.42 kJ m−3 and 37.61 kJ m−3, respectively. This reduction is attributed to the effects of NBT substitution, which altered the polarization response. Notably, a significant decrease in maximum polarization was observed for BZT-Pr–0.05NBT, explaining its superior energy conversion performance.31 To highlight the performance of BZT-Pr–0.05NBT ceramics relative to other ferroelectrics reported in the literature, Table 4 provides a comprehensive comparison of the energy harvesting capabilities of various ferroelectric materials using the Olsen cycle under different operating conditions. Table 4 shows that the response of BZT-Pr–0.05NBT ceramics is comparable to that of some materials. However, enhanced energy conversion capabilities have been observed in other modified ferroelectric compositions achieved through A-site or B-site substitution. For instance, adding 0.5% Fe to BCT–BZT–Fe ceramics improved the energy density and a maximum of 305 kJ m−3 was reached under conditions of 30–110 °C and 0–30 kV cm−1.54 Similarly, Sr substitution in BCT–BZT (BCT–BZT–Sr) resulted in a 200% increase in energy density, with 15% Sr yielding a maximum of 300 kJ m−3 over a 30–90 °C temperature range.28 For BSTZS ceramics, substituting (Ti, Zr) with Sn significantly improved energy conversion, with a 5% Sn content achieving the highest energy density of 175 kJ m−3 under 25–110 °C and 0–30 kV cm−1.31 These findings highlight the critical role of composition tailoring in optimizing the thermal energy harvesting performance of ferroelectric materials.


image file: d5ra00245a-f11.tif
Fig. 11 Energy conversion versus electric field plots for (BZT-Pr–xNBT) samples with NBT content of (a) x = 0.05, (b) x = 0.10 and (c) x = 0.15 at EL = 0 kV cm−1. In addition, panel (d) compares the energy conversion capacity across the studied compositions.
Table 4 Energy density (ND) of various ferroelectric materials
Materials ND (kJ m−3) TLTH (°C) ΔT (°C) ELEH (kV cm−1) ΔE (kV cm−1) Ref.
BZT-Pr–0.05NBT 80.77 30–140 110 0–30 30 Current work
PZT 92 25–160 135 1–20 19 55
PZST 100 146–159 13 0–29 29 52
Ba0.85Ca0.15Ti0.895Fe0.005Zr0.1O3 101 30–110 80 0–10 10 54
Ba0.85Ca0.15Ti0.895Fe0.005Zr0.1O3 305 30–110 80 0–30 30 54
Ba0.85Ca0.15Ti0.9Zr0.1O3 108 30–90 60 0–20 20 28
Ba0.85Sr0.15Ti0.9Zr0.1O3 300 30–90 60 0–30 30 28
(Ba0.95Sr0.05)(Ti0.95Zr0.05)O3 32 25–60 35 0–30 30 31
(Ba0.95Sr0.05)(Ti0.95Zr0.05)0.975Sn0.05O3 175 25–110 85 0–30 30 31
0.5Ba(Zr0.2Ti0.8)O3–0.5(Ba0.7Ca0.3)TiO3 87 20–120 100 1–10 9 56


4. Conclusion

To conclude, the structural, optical and ferroelectric properties of BZT-Pr–xNBT (x = 0.05, 0.10 and 0.15) ceramics, successfully prepared by a conventional solid-state reaction technique, were investigated. The atomic positions, bond distances, and bond angles of these compositions, which crystallized in the tetragonal structure, have been determined. Also, the crystallite size of each composition was calculated by two methods. Direct band gap values have been estimated using UV-Vis diffuse reflectance spectrum for these ceramics. Moreover, in this study, we systematically examined the ferroelectric and energy storage properties of BZT-Pr–xNBT ceramics with varying NBT contents (x = 0.05, 0.10 and 0.15) to evaluate their potential for energy storage and pyroelectric energy harvesting. Our findings reveal that increasing the NBT content significantly alters the ferroelectric characteristics of the material. Among the compositions, BZT-Pr–0.15NBT exhibited the highest energy storage properties at room temperature, with Wrec of 65.86 mJ cm−3 and η of 69.25%. This composition also demonstrated remarkable thermal stability across a wide temperature range (30–140 °C). Furthermore, the scaling relations for EC and Pr as functions of temperature were derived. For EC versus T, the scaling relations were determined as ECT−0.243, ECT−0.190 and ECT−0.089 for samples with NBT contents of x = 0.05, x = 0.10 and x = 0.15 respectively. Similarly, the scaling relations for Pr versus T were PrT−0.474, PrT−0.385, and PrT−0.233. These results indicate that BZT-Pr–0.15NBT exhibits improved thermal stability, with weaker temperature dependence of both EC and Pr which justifies the high temperature stability of the sample. Additionally, the pyroelectric energy harvesting performance was assessed under the Olsen cycle. Among the compositions, BZT-Pr–0.05NBT demonstrated the highest energy conversion efficiency, making it a promising candidate for thermal energy harvesting applications. These results underscore the critical role of NBT substitution in tailoring the ferroelectric and pyroelectric properties of BZT-Pr–xNBT ceramics.

Data availability

The data that support the findings of this study are available from the corresponding author upon a reasonable request.

Conflicts of interest

There are no conflicts of interest to declare.

Acknowledgements

The authors acknowledge the support received from the Tunisian Ministry of Higher Education and Scientific Research.

References

  1. T. Mikolajick, U. Schroeder and S. Slesazeck, IEEE Trans. Electron Devices, 2020, 67, 1434–1443 CAS.
  2. S. Gomasu, S. Saha, S. Ghosh, R. Bhowmik and D. Das, ACS Appl. Mater. Interfaces, 2024, 16, 3654–3664 CrossRef CAS.
  3. P. K. Panda and B. Sahoo, Ferroelectrics, 2015, 474, 128–143 CAS.
  4. Y. Liu, Q. Li, L. Qiao, Z. Xu and F. Li, Adv. Sci., 2022, 9, 2204631 CAS.
  5. S. Zhang, B. Malič, J. F. Li and J. Rödel, J. Mater. Res., 2021, 36, 985–995 CAS.
  6. W. Ding, J. Lu, X. Tang, L. Kou and L. Liu, ACS Omega, 2023, 8, 6164–6174 CAS.
  7. J. Hoffman, X. Hong and C. H. Ahn, Nanotechnology, 2011, 22, 254014 CAS.
  8. Q. Yuan, F. Yao, Y. Wang, R. Ma and H. Wang, J. Mater. Chem. C, 2017, 5, 9552–9558 RSC.
  9. R. Ahmed, R. Si, S. ur Rehman, Y. Yu, Q. Li and C. Wang, Results Phys., 2021, 20, 103623 CrossRef.
  10. H. Chen, C. Yang, C. Fu, J. Shi, J. Zhang and W. Leng, J. Mater. Sci.: Mater. Electron., 2008, 19, 379–382 CrossRef CAS.
  11. Y. Yuan, S. R. Zhang, X. H. Zhou, B. Tang and B. Li, J. Electron. Mater., 2009, 38, 706–710 CAS.
  12. H. I. Hsiang, L. T. Mei and Y. J. Chun, J. Am. Ceram. Soc., 2009, 92, 2768–2771 CAS.
  13. V. Batra, S. Kotru, M. Varagas and C. V. Ramana, Opt. Mater., 2015, 49, 123–128 CrossRef CAS.
  14. Y. Qi, L. Zhang, G. Jin, Y. Wan, Y. Tang, D. Xu and D. Sun, Ferroelectrics, 2014, 458, 64–69 CAS.
  15. Y. Bai, T. Siponkoski, J. Peräntie, H. Jantunen and J. Juuti, Appl. Phys. Lett., 2017, 110, 063903 Search PubMed.
  16. L. Tang, S. Han, Y. Ma, Y. Liu, L. Hua, H. Xu and J. Luo, Chem. Mater., 2022, 34, 8898–8904 CAS.
  17. V. Segouin, M. Domenjoud, Y. Bernard and L. Daniel, Acta Mater., 2021, 211, 116870 CAS.
  18. V. Athikesavan, S. Bhuvana and G. Thilakavathi, Ferroelectr. Lett. Sect., 2022, 49, 104–110 CrossRef CAS.
  19. D. Thuau, K. Kallitsis, F. D. D. Santos and G. Hadziioannou, J. Mater. Chem. C, 2017, 5, 9963–9966 RSC.
  20. M. Valant, Prog. Mater. Sci., 2012, 57, 980–1009 CAS.
  21. L. F. Zhu, B. P. Zhang, L. Zhao and J. F. Li, J. Mater. Chem. C, 2014, 2, 4764–4771 CAS.
  22. S. Khardazi, H. Zaitouni, S. Belkhadir, D. Mezzane, M. Amjoud, Y. Gagou and S. Terenchuk, J. Solid State Chem., 2022, 311, 123112 CAS.
  23. M. Maraj, W. Wei, B. Peng and W. Sun, Materials, 2019, 12, 3641 CAS.
  24. S. K. Upadhyay, V. R. Reddy, P. Bag, R. Rawat, S. M. Gupta and A. Gupta, Appl. Phys. Lett., 2014, 105, 112907 CrossRef.
  25. G. Liu, W. Yu, Y. Wang, H. Feng, M. Hao, G. Wu and Y. Yan, Ceram. Int., 2023, 49, 34387–34396 CrossRef CAS.
  26. A. P. Sharma, D. K. Pradhan, S. K. Pradhan and M. Bahoura, Sci. Rep., 2019, 9, 16809 CrossRef CAS PubMed.
  27. A. Mandal, D. Yadav, S. K. Mittal, U. Jamwal, D. Kaneria, A. Khokhar and K. L. Yadav, J. Alloys Compd., 2024, 995, 174817 CAS.
  28. S. Patel, D. Sharma, A. Singh and R. Vaish, J. Materiomics, 2016, 2, 75–86 CrossRef.
  29. C. Zhou, Q. Zhang, W. Cai, R. Yang, S. Chen, R. Gao and J. Sun, Mater. Sci. Semicond. Process., 2022, 150, 106950 CAS.
  30. W. Cai, Q. Zhang, C. Zhou, R. Gao, S. Zhang, Z. Li and C. Fu, J. Mater. Sci.: Mater. Electron., 2020, 31, 9167–9175 CrossRef CAS.
  31. I. Djemel, M. H. Khedhri, I. Kriaa, Z. Sassi, L. Seveyrat, N. Abdelmoula and H. Khemakhem, Appl. Phys. A, 2023, 129, 784 CrossRef CAS.
  32. I. Kriaa and A. Maalej, J. Energy Storage, 2024, 102, 114044 Search PubMed.
  33. I. Zouari, Z. Sassi, L. Seveyrat, N. Abdelmoula, L. Lebrun and H. Khemakhem, J. Alloys Compd., 2020, 825, 153859 CAS.
  34. J. Tauc, Optical properties and electronic structure of amorphous semiconductors, in Optical Properties of Solids: Papers from the NATO Advanced Study Institute on Optical Properties of Solids Held, Springer US, Boston, MA, at Freiburg Germany, 1966, pp. 123–136 Search PubMed.
  35. J. Rodriguez-Carvajal, A Rietveld and Pattern Matching Analysis Program (Version April 2008), Laboratoire Léon Brillouin (CEA-CNRS), France, 2008.
  36. G. K. Williamson and W. H. Hall, Acta Metall., 1953, 1, 22–31 CAS.
  37. J. Kanamori, J. Appl. Phys, 1960, 31, S14–S23 Search PubMed.
  38. H. Bai and X. Liu, Mater. Lett., 2013, 100, 1–3 Search PubMed.
  39. T. Badapanda, S. Sarangi, B. Behera, P. K. Sahoo, S. Anwar, T. P. Sinha and L. S. Cavalcante, Curr. Appl. Phys., 2014, 14, 708–715 CrossRef.
  40. C. C. Parhi, A. B. Thirumalasetty, A. R. James and M. Wuppulluri, ACS Omega, 2023, 8, 37752–37768 CrossRef CAS PubMed.
  41. S. Patel, A. Chauhan and R. Vaish, Int. J. Appl. Ceram. Technol., 2015, 12, 899–907 CrossRef CAS.
  42. M. K. Bilal, J. Wang, R. Bashir, H. Liu, S. U. Asif, J. Xie and W. Hu, J. Am. Ceram. Soc., 2021, 104, 3982–3991 CrossRef CAS.
  43. Y. Sun, H. Liu, H. Hao, S. Zhang, L. Guo and Z. Yu, Ceram. Int., 2012, 38, S41–S44 CrossRef CAS.
  44. H. I. Humburg, M. Acosta, W. Jo, K. G. Webber and J. Rödel, J. Eur. Ceram. Soc., 2015, 35, 1209–1217 CrossRef CAS.
  45. K. Chen, J. Ma, J. Wu, X. Wang, F. Miao, Y. Huang, C. Shi, W. Wu and B. Wu, J. Mater. Sci.: Mater. Electron., 2020, 31, 12292–12300 CAS.
  46. N. Li, K. Xie and H. Huang, ACS Energy Lett., 2023, 8, 4357–4370 CAS.
  47. A. R. Jayakrishnan, J. P. B. Silva, K. Kamakshi, D. Dastan, V. Annapureddy, M. Pereira and K. C. Sekhar, Prog. Mater. Sci., 2023, 132, 101046 CAS.
  48. Q. Li, W. Zhang, C. Wang, L. Ning, Y. Wen, B. Hu and H. Fan, J. Alloys Compd., 2019, 775, 116–123 CAS.
  49. F. Z. Yao, Q. Yuan, Q. Wang and H. Wang, Nanoscale, 2020, 12, 17165–17184 RSC.
  50. M. Zhou, R. Liang, Z. Zhou, S. Yan and X. Dong, ACS Sustain. Chem. Eng., 2018, 6, 12755–12765 CrossRef CAS.
  51. S. Patel, A. Chauhan, A. Chauhan and R. Vaish, Mater. Res. Express, 2015, 2, 035501 CrossRef.
  52. R. B. Olsen, D. A. Bruno and J. M. Briscoe, J. Appl. Phys., 1985, 58, 4709–4716 CAS.
  53. R. B. Olsen and D. A. Bruno, Pyroelectric conversion materials, in IECEC 1986, Proceedings of the Twenty-First Intersociety Energy Conversion Engineering Conference, 1986, vol. 1, pp. 89–93 Search PubMed.
  54. D. Sharma, S. Patel, A. Singh and R. Vaish, J. Asian Ceram. Soc., 2016, 4, 102–111 Search PubMed.
  55. S. Patel, A. Chauhan and R. Vaish, Mater. Res. Express, 2014, 1, 025504 CAS.
  56. G. Vats, A. Chauhan and R. Vaish, Int. J. Appl. Ceram. Technol., 2015, 12, E49–E54 CAS.

This journal is © The Royal Society of Chemistry 2025
Click here to see how this site uses Cookies. View our privacy policy here.