3.1 Compositional heterogeneity and microstructure
To analyze the lattice structures of BNT-BZT-xGaN samples, the XRD patterns and XRD Rietveld refinement are presented in Fig. 3(a)-(d). As shown in Fig. 3(a), the pure perovskite structures can be observed for all BNT-BZT-xGaN samples, suggesting that GaN has diffused into BNT-BZT lattice (as shown in Fig. 2(d)), being consistent with the analysis of EDS. The perovskite phase structure of BNT-BZT-xGaN ceramics can be confirmed by assessing the splitting peaks at around 40.0o and 46.5o as magnified in Fig. 3(b). For the rhombohedral phase, only one peak can be observed at 46.5o, while the (200) peak turns into two peaks of tetragonal phase [25–26]. As shown in Fig. 3(b), both the (100) peak at 40.0o and the (200) peak divide into two peaks, implying the coexistence of rhombohedral and tetragonal phase in BNT-BZT-xGaN system [27]. The rhombohedral phase dominates in the BNT-BZT-xGaN ceramic matrix because of the agglomeration of BaTiO3 as shown in Fig. 2. The main crystal phase with R3c space group (as shown in Fig. 3(d)) can be identified by the XRD Rietveld refinement as shown in Fig. S5(a)-(e). The reliability factors of weighted patterns (Rwp), patterns (Rp) and the goodness-of-fit indexes (χ2) are under 7.5%, 5.5% and 5.5, respectively, indicating that the crystal structure mode is valid and the refinement result is reliable. The lattice parameters (a = b = c, α = β = γ and volume) and the distance of Ti-O are obtained as shown in Fig. 3(c). The angles (α = β = γ) tend to decrease, while both the lattice parameters (a, b, c) and volumes first increase and then decrease with addition of GaN. In addition, the similar tendency is also observed in the Ti-O bond distances of BNT-BZT-xGaN samples. As the GaN content increases to 0.1 wt%, the Ti-O (1) and Ti-O (2) bond lengths reach the maximum value of 1.894 Å and 2.044 Å.
3.2 Lattice vibrations and phonon structure
In BNT-based dielectric solids, lattice vibrations dominate the heat transport and a phonon refers to the quantum of atomic vibrational energy. The thermal transport in solids is regarded as the diffusion of phonons actuated by a heat gradient [28]. To analyze the relation between lattice vibration and phonon structure of BNT-BZT-xGaN samples, the calculated phonon-dispersion curves along the high-symmetry lines of Brillouin zone of BNT-BZT and BNT-BZT-xGaN with x = 0.1 wt% are presented in Fig. S6. There are three low-frequency acoustic modes observed near the Г-point and a lack of imaginary line in phonon band structures proves the dynamical stability for both BNT-BZT and BNT-BZT-xGaN with x = 0.1 wt% samples [29–30]. Meanwhile, the dispersion curves show some optical modes above 80 cm− 1 and a band gap between 600 ~ 680 cm− 1. The optical modes below 800 cm− 1 and the acoustic modes determine the lattice thermal conductivity, playing an important role in heat transport. To find out intrinsic nature of lattice thermal conduction, the phonon density of state (PDOS) of BNT-BZT and BNT-BZT-xGaN with x = 0.1 wt% is presented in Fig. 4(a)-(b). As shown in Fig. 4(a)-(b), there are different shadowed areas representing Bi, Na, Ti, O, Ga and N atom phonon DOS, respectively. The atoms of BNT-based matrix mainly possess phonon distribution within the low-frequency part blew 800 cm− 1. This means that atoms for BNT-based matrix vibrate at the phonon modes below 800 cm− 1 and conduct heat, corresponding to the phonon band dispersion as shown in Fig. S6(b). In addition, the phonon distribution of Ga atom is located between 10 and 375 cm− 1, which has overlapped regions with the atom phonon DOS of Na, Ti, and O. These overlapped phonon distributions stand for the resonance vibration between Ga and Na, Ti, O atoms, which contributes to the enhancement of the lattice vibration of BNT-BZT. Thus, the enhancement of the lattice vibration of BNT-BZT not only gives rise to the enhancement of the lattice heat conduction, but also improves the spontaneous polarization of the samples. Because the spontaneous polarization of BNT-based ceramics is prevailingly derived from the movement of Ti and O atoms, the resonance vibration among Ga and Ti, O atoms (as shown in Fig. 4(b)) improve the vibration amplitude of Ti-O (as shown in Fig. 4(f) and supplementary Video V1), which can be confirmed by the group velocity of BNT-BZT and BNT-BZT-xGaN with x = 0.1 wt% as shown in Fig. 4(c)-(d). The phonon group velocity of BNT-BZT-xGaN with x = 0.1% at about 300 cm− 1, roughly corresponding to the Ti-O vibration (shown in Fig. 4(e) and Fig. S7), is 684.1 m s− 1. This value is higher than 651.0 m s− 1 of pristine BNT-BZT sample as shown in Fig. 4(c). Meanwhile, as shown in Fig. 4(e) and Fig. S7, the intensity of Raman spectrum between 10 and 375 cm− 1 for BNT-BZT-xGaN with x = 0.1% is significantly stronger than that of pure BNT-BZT sample, further confirming the resonance vibration among Ga and Ti, O atoms. These results show that the introduction of GaN achieves the simultaneous improvement of thermal conductivity and spontaneous polarization, which can be verified by the measurement of heat conductivity and electrical properties as shown in Fig. 5 and Fig. 6, respectively.
3.3 Thermal conductivity and heat transport
Based on the kinetic theory of gases and Debye’s specific heat theory, the lattice thermal conductivity of a crystalline solid can be expressed as: \(\kappa =\frac{1}{3}{C}_{v}{v}_{g}^{2}\tau\), where \({C}_{v}\) is the specific heat per unit volume, \({v}_{g}\) is the phonon group velocity, \(\tau\)is the relaxation time [31]. According to the equation, the thermal conductivity is proportional to the square of the phonon group velocity. As discussed in Fig. 4(c)-(d), with the modification of GaN, the phonon group velocity of BNT-BZT-xGaN with x = 0.1 wt% is higher than that of pristine BNT-BZT sample. Theoretically, the thermal conductivity of BNT-BZT-xGaN with x = 0.1 wt% is correspondingly higher than that of pristine BNT-BZT sample. To evaluate the actual heat conductivity, Fig. 5(a) and Fig. S8 give the temperature-dependent thermal parameters of BNT-BZT-xGaN ceramics with different amount of GaN. As shown in Fig. 5(a), while the thermal conductivity dose not vary much with an increase in temperature, the introduction of GaN contributes to an increase in the thermal conductivity from 1.48 W m− 1 K− 1 to 1.61 W m− 1 K− 1 when the GaN content increases from 0 to 0.1 wt% at room temperature. This may be contributed to the resonance vibration and enhancement in the phonon group velocity resulted from introduction of GaN. In addition, the temperature distributions of the pristine and BNT-BZT-xGaN sample with x = 0.1 wt% calculated by COMSOL Multiphysics software at the same time duration are presented in Fig. 5(b). It can be observed that BNT-BZT-xGaN sample with x = 0.1 wt% is able to conduct heat faster for the same thermal input in comparison with the pure ceramic. Namely, the construction of BNT-BZT-xGaN hybrid ceramic can be an effective method to improve the heat transfer performance compared to the pure BNT-BZT ceramic samples.
3.4 Dielectric and ferroelectric properties
Figure 6(a) and Fig. S9 present the dielectric constant (εr) and loss (tanδ) with temperature range from 25 to 400 oC for polarized BNT-BZT-xGaN samples. As shown in Fig. 6(a) and Fig.S9, three dielectric peaks can be observed in dielectric constant curves for all samples, which correspond to ferroelectric to paraelectric phase transition (Tm) at around 300 oC, ergodic relaxor to ferroelectric state (Ts) at about 150 oC, and non-ergodic to ergodic relaxor state (TFR) near ambient temperature, respectively [32]. The Tm decreases slightly while the Ts and TFR increase appropriately with the increase in GaN content as presented in Fig. 5(b). In particular, the temperature of non-ergodic to ergodic relaxor state (TFR) for BNT-BZT-xGaN samples shifts from 40 to 65 oC with increasing GaN, which is corresponding to the peak pyroelectric temperature as shown in Fig. 7(a). In addition, the dielectric constant (εr) and loss (tanδ) at corresponding ambient temperature range decrease with the increase in GaN content, especially in room temperature from ~ 1750 to ~ 1400 and from ~ 0.074 to ~ 0.05, respectively. Because the pyroelectric energy harvesting figure of merit is inversely proportional to the dielectric constant, the decrease in dielectric constant (εr) facilitates the improvement of pyroelectric energy harvesting performance. The polarization-electric field (P-E) hysteresis loops of BNT-BZT-xGaN samples with various content of GaN at room temperature are presented in Fig. 6(c). All the samples exhibit typical ferroelectric hysteresis loops. Both the saturation polarization (Ps) and remnant polarization (Pr) first increase and then decrease with increasing GaN content. As the GaN content increases to 0.1 wt%, both the Ps and Pr gradually increase to the maximum value of 40 µC cm− 2 and 23 µC cm− 2, respectively, which is attributed to the enhancement of Ti-O vibration as discussed in Fig. 4.
3.5 Pyroelectric effect and mechanism of pyroelectric energy harvesting
Figure 7(a) presents the temperature dependent pyroelectric coefficient of BNT-BZT-xGaN ceramics with different content of GaN in the temperature range of 20 ~ 90 oC. It can be observed from Fig. 7(a) that the peak pyroelectric coefficient of BNT-BZT-xGaN samples shows a similar trend as that of saturation polarization and remnant polarization. The peak pyroelectric coefficient firstly increases and then decreases with the addition of GaN, and reaches the maximum value of 850×10− 4 C m− 2 K− 1 when the content of GaN increases to 0.1 wt%. Meanwhile, the peak position gradually shifts to higher temperature with the increase in GaN, in accordance with the temperature of non-ergodic to ergodic relaxor state as shown in Fig. 7(b). The mechanisms of pyroelectric energy harvesting for BNT-BZT and BNT-BZT-xGaN with x = 0.1 wt% are illustrated in Fig. 7(c) and Fig. 7(d), respectively. At a steady state (dT/dt = 0), for BNT-BZT samples, the polarized dipoles reorient along the applied electric field and oscillate randomly to reach an equilibrium condition, concomitantly obtaining the attracted charges on both surfaces of samples as shown in Fig. 7(c). For BNT-BZT-xGaN with x = 0.1 wt%, the poled dipoles are easier to reorient and absorb a greater quantity of opposite electrical charge on both sides of samples (as presented in Fig. 7(d)) due to the enhancement of Ti-O vibration, which provides the potential for giant pyroelectric effect. When the sample is heated, the dipoles in BNT-BZT ceramics oscillate within a large degree of alignment, leading to the decrease in spontaneous polarization and absorbed electrical charges, and thus the reduced currents flow in the external circuit. With the introduction of GaN, heat transfer is faster in BNT-BZT-xGaN with x = 0.1 wt% owing to the improved thermal conductivity, which contributes to dipoles’ oscillation in even larger degree and larger currents in the circuits as illustrated in Fig. 7(d).
3.6 Pyroelectric energy harvesting
Figure 8 presents the energy harvesting performance of BNT-BZT-xGaN samples in response to temperature change between 25 to 50 oC. The corresponding time dependent temperature variation curves are exhibited in Fig. 8(a). It can be seen that the thermal distribution is a combination of continuous rising trend with periodic temperature fluctuation of about 2 oC every 10 s. The short-circuit current and open-circuit voltage of BNT-BZT-xGaN ceramic induced by the low temperature gradient of 2 oC are shown in Fig. 8(b)-(c). As shown in Fig. 8(b)-(c), both short-circuit current and open-circuit voltage increase first and then decrease with the increase in GaN content. As the content of GaN reaches 0.1 wt%, the BNT-BZT-xGaN sample with x = 0.1 wt% can generate a short-circuit peak current of 0.12 µA, or an open-circuit peak voltage of 58 V. The improvements in pyroelectric short-circuit current and open-circuit voltage are mainly owing to the thermoelectrical coupling effect in BNT-BZT-xGaN sample with x = 0.1 wt% modulated by GaN. Meanwhile, the output energy density of BNT-BZT-xGaN ceramics with x = 0.1 wt% has been investigated by measuring the pyroelectric current and voltage with a range of load resistance and calculation based on the integral formula \(\text{w}={\int }_{t=0}^{t=T}\frac{I\times U}{S\times d}dt\), where I, U, T, S and d are pyroelectric current, voltage, cycle time, sample area and thickness, respectively. The pyroelectric current and voltage obtained at different load resistance are presented in Fig. 7(d)-(e). Figure 7(f) shows the load resistance dependent energy density curves of BNT-BZT-xGaN ceramics with x=0.1 wt% induced by the low temperature gradient of 2 oC. As shown in Fig. 7(f), the output energy density firstly increases and then decreases with the increase in load resistance and reaches the maximum value of 80 µJ cm−3 at the load resistance of 600 MΩ. In addition, Fig. 7(g) presents the electrical circuit diagram of the pyroelectric energy harvesting device where a 4.8 µF load capacitor and BNT-BZT-0.1 wt%GaN sample are connected by a phase-controlled rectifier. The pyroelectric samples generate alternating current driven by a low thermal gradient, which is converted by the phase-controlled rectifier to direct current and stores the electricity in the capacitor. When the voltage of the capacitor reaches 2.2 V, the capacitor is triggered to discharge and light up the LED bulb. The charging-discharging voltage variation curves of the capacitor are presented in Fig. S10. Also, the practical applicability of the thermal energy harvesting device induced by low temperature gradient is illustrated in supplementary Video V2. The BNT-BZT-xGaN ceramics with x = 0.1 wt% exhibit superior pyroelectric energy density driven by unit temperature gradient in comparison with state-of-the-art pyroelectric energy harvesting systems (as shown in Table 1), which have great potential application in low temperature driven thermal energy harvesting devices [16, 33–45].
Table 1
Comparison of pyroelectric energy harvesting performance of BNT-BZT-xGaN ceramic with that of other state-of-the-art pyroelectric energy harvesting systems driven by unit temperature gradient.
Materials
|
ΔT
(K)
|
p
(10− 4 C m− 2 K− 1)
|
Uo
(V K− 1)
|
Io
(µA K− 1)
|
Energy
density for volume
(µJ cm− 3 K− 1)
|
Power
Density for area
(µw cm− 2 K− 1)
|
Ref.
|
ITO/BNT-BZT/Ag
|
30.5
|
5.0
|
1.67
|
0.003
|
~
|
~
|
[33]
|
BZT/BCT/STO
|
11.8
|
34.3
|
~
|
0.11
|
~
|
~
|
[34]
|
I−/I3−/MC/KCl (TGC)
|
15.0
|
~
|
0.01
|
~
|
5.34
|
~
|
[35]
|
PPGO
|
110.0
|
~
|
0.02
|
0.10
|
~
|
0.17
|
[36]
|
PMN-PMS-PZT: xCNT
|
20.0
|
43.3
|
0.67
|
~
|
~
|
~
|
[16]
|
HfO2/Hf0.5Zr0.5O2:La
|
160.0
|
0.7
|
~
|
~
|
~
|
~
|
[37]
|
PyEHs
|
4.0
|
~
|
15
|
~
|
~
|
~
|
[38]
|
PEG
|
70.0
|
~
|
0.43
|
0.27
|
~
|
0.09
|
[39]
|
SnS:Na
|
3.9
|
~
|
0.001
|
24.51
|
~
|
~
|
[40]
|
PEDOT: Tos
|
13.8
|
~
|
5.41
|
0.001
|
~
|
0.03
|
[41]
|
WKF-based PTM
|
80.0
|
~
|
~
|
~
|
~
|
0.002
|
[42]
|
S-PENG
|
29.0
|
~
|
4.14
|
0.072
|
~
|
0.009
|
[43]
|
Al/PVDF/Al
|
5.0
|
17.9
|
9.60
|
0.005
|
0.23
|
~
|
[44]
|
PPNG
|
38.0
|
~
|
0.001
|
~
|
~
|
0.003
|
[45]
|
BNT-BZT-GaN
|
2.0
|
850.0
|
29.0
|
0.06
|
40
|
~
|
This work
|