1. Introduction
Recently, the electrocaloric effect (ECE) has attracted great interest, as it could be applied for designing active cooling elements, especially in micro and power electronics, and environmentally friendly air-conditioning and cooling systems. As a consequence, a large number of electrocaloric effect studies have been published (see a recent review by Barman et al. [
1]). They bring a certain optimism regarding the prospect of implementation of electrocaloric cooling. First of all, it concerns ECE obtained in thin films. The promising results are obtained at high electric fields not suitable for bulk materials [
2,
3]. At the same time, the extremely small heat capacity of thin films significantly reduces practical application opportunities. Nevertheless, the possibility of such application is already demonstrated in PVDF thin film [
4].
Regarding ECE in bulk ferroelectrics, two cases should be considered separately: (1) moderate electric fields and (2) high electric fields. In the easily accessible moderate electric field range (20–30 kV/cm), only a few compositions have been found with values of ECE temperature change Δ
T exceeding 1 °C. All of these are related to entropy change at a first-order phase transition, which is induced by an electric field. One of these compositions is Pb
0.99Nb
0.02(Zr
0.70Sn
0.20Ti
0.05)
0.98O
3, which belongs to Pb(Sn,Zr,Ti)O
3 solid solutions and has the largest values of Δ
T found until now in this electric field range (2.6 °C at 30 kV/cm) [
5]. It is established that the reason for such high Δ
T values is the electric field-induced phase transition between ferroelectric and paraelectric states around 160 °C. Values of Δ
T at the first-order phase transition between ferroelectric and antiferroelectric states, observed for this composition at lower temperatures, as well as for other Pb(Sn,Zr,Ti)O
3 compositions, are much lower [
6,
7]. Unfortunately, the temperature region of the most expressed ECE in this case is not appropriate for cooling devices. PbSc
1/2Ta
1/2O
3 (PST)-based compositions are much more promising from this point of view. It is found that, at the first-order phase transition between ferroelectric and paraelectric states in partly B-site-ordered PST, observed in the temperature range from 0 to 10 °C, Δ
T exceeds 1 °C at 20 kV/cm [
8]. Later, increased B-site ordering and modification of PST with Sb allowed a Δ
T increase of up to 1.5–1.8 °C at 25 kV/cm [
9,
10]. Still, this value is below the level necessary for implementing a promising prototype, which goes beyond mere demonstration of the working principle. At the moment, it seems that there are no suggestions for seeking out new materials with higher values of Δ
T in this electric field range.
Another way to increase Δ
T is to apply higher electric fields. Indeed, in PST, at 50 kV/cm, Δ
T reached 2.3 °C [
9]. Promising values of
ΔT even at higher electric fields, exceeding 100 kV/cm, were reported in recent years. This was achieved by reducing the sample thickness and creating multilayers. In such a case, the reduced thickness allows one to increase the applied electric field, while the multilayer structure compensates for reduced heat capacity of a thin single layer. In 0.9PbMg
1/3Nb
2/3O
3-0.1PbTiO
3 ceramics at 160 kV/cm, Δ
T reaches 3.5 °C [
11]. In Ba(Ti
0.8Zr
0.2)O
3 ceramics at 145 kV/cm, Δ
T is 4 °C [
12], while in a multilayer structure of the same composition, even Δ
T of 6 °C was achieved at 150 kV/cm [
13]. Taking into account the large passive heat capacity in multilayer structures, which does not contribute to ECE, but consumes part of the generated heat [
14], such a high value of Δ
T is surprising. Δ
T of 4 °C is reported also in a BaTiO
3 multilayer structure, but at two very different applied electric field values—176 and 350 kV/cm [
15,
16]. In a PST multilayer structure, Δ
T of 5.5 °C was observed at 290 kV/cm [
14]. Lower values are obtained in PbMg
1/3Nb
2/3O
3—at 90 kV/cm, Δ
T of 2.4 °C was reported by Rožič et al. [
17], and only 1.5 °C in Peräntie et al. [
18]. Such results bring optimism about increasing Δ
T by increasing of applied electric field. Even though the obtained values of Δ
T are remarkably lower compared with thin films, they are usually within a range useful for development of the first prototype coolers with real output power, assuming problems related with breakdown fields are resolved. Unfortunately, measurements of Δ
T in bulk ceramics at high electric fields are published rarely. Moreover, sometimes they are not mutually consistent, which can be seen also from the results mentioned in the previous paragraph.
The experimental results described above are obtained via direct measurement of the temperature change created by ECE (heat flow measurements by adapted DSC are also used). As follows from a large number of published results, another method—the indirect approach—is even more widely used for the study of ECE. In this case, Δ
T values are calculated using the Maxwell’s relations from electric field-temperature dependences of polarization
P(E,T), which are extracted from experimentally measured polarization hysteresis loops. The great popularity of such an approach is understandable—the experimental part is much easier to implement. It is especially valuable for studies in thin films, where performing direct ECE measurements is extremely complicated. Promising values of Δ
T, obtained by the indirect method in PbZr
0.95Ti
0.05O
3 thin film [
19], alone started the “renaissance” in ECE research.
In Na
1/2Bi
1/2TiO
3 (NBT)-based materials, Δ
T values above 1 °C, as calculated by an indirect method, are reported repeatedly, paying more attention to the morphotropic phase boundary of NBT–BaTiO
3 solid solutions. Controversial results are obtained in a pure 0.94NBT–0.06BaTiO
3 composition: in [
20], the minimum of temperature dependence of Δ
T (∆
T = −2.3 °C at 70 kV/cm) is found at 50 °C; starting at room temperature, only positive, monotonously increasing values of Δ
T at higher temperatures were measured by Li et al. [
21]. Modification of this composition by 0.5 wt% of La [
22] facilitated improved Δ
T (∆
T = −2.6 °C at 50 kV/cm, 65 °C). A similar value in the same temperature range was obtained in 0.9NBT–0.1K
1/2Bi
1/2TiO
3 [
23]. ECE in NBT, doped with various lanthanides, is surprisingly sensitive to the choice of doping element [
24], in particular for doping with 2 mol% Gd—reaching Δ
T = 1.8 °C (90 kV/cm, 150 °C). Positive and negative ECE was frequently reported in the same composition across different temperature ranges [
20,
21,
22].
Δ
T values obtained by a direct method are lower. Le Goupil and Alford [
25] have observed Δ
T = 1.5 °C at 70 kV/cm, 125 °C in a 0.94NBT–0.06KNbO
3 solid solution. A similar value was obtained at 100 °C in 0.94NBT–0.06BT, where 1.5 mol% Bi was replaced by Sm [
26].
The correspondence between the direct and the indirect methods is discussed repeatedly. Sometimes both methods give comparable results [
27,
28,
29], although sometimes obtained results are clearly different [
29,
30,
31]; even the sign of Δ
T may differ [
30,
32]. Nonegodicity [
31] and the presence of polar nanoregions (PNRs) [
29,
33] are used to explain such nonconformance. Violation of conditions necessary to apply the Maxwell’s relations due to polydomain state existing in the studied materials is also mentioned [
30,
34]. Usability of polarization data should be treated with extreme caution if they are extracted from polarization hysteresis loops in bipolar electric fields, because reorientation of domains is essential for this kind of loops. Polarization hysteresis loops in the ferroelectric state should be at least with well-expressed saturation of polarization, while in the nonpolar phase, they should have low contribution of dielectric losses. Ignoring these requirements yields disputable results, such as negative ECE, sometimes obtained in Na
1/2Bi
1/2TiO
3-based compositions [
22,
23,
35,
36], or physically unrealistically high values of Δ
T [
37].
Taking into account the above considerations, extending the range of ECE results obtained by the direct method across a higher range of electric fields is of great importance. The presented research reports on detailed studies of ECE temperature change Δ
T in (1−x)(0.8Na
1/2Bi
1/2TiO
3-0.2BaTiO
3)−xCaTiO
3 ((1−x)(0.8NBT-0.2BT)−xCT) solid solutions in electric field range up to 100 kV/cm. The initial composition of this group—0.8NBT-0.2BT (x = 0.00)—is located in the tetragonal side from the morphotropic phase boundary in the phase diagram of NBT-BT solid solutions. This choice is motivated by well-expressed ferroelectric–nonferroelectric phase transition in this composition with the largest polarization and the highest jump of lattice parameter comparing to the whole NBT-BT concentration range [
38]. CaTiO
3 was added in order to reduce the depolarization temperature and increase Δ
T around and below room temperature, which is a relevant condition considering the application of ECE in cooling devices.
2. Materials and Methods
(1−x)(0.8Na0.5Bi0.5TiO3-0.2BaTiO3)−xCaTiO3 ((1−x)(0.8NBT-0.2BT)−xCT) ceramics with various CaTiO3 concentrations (x = 0.050, 0.075, 0.100, and 0.125) were prepared by solid-state reaction from chemical-grade oxides and carbonates: Bi2O3 (purity 99.9%), TiO2 (99.8%), Na2CO3 (99%), BaCO3 (≥99%), and CaCO3 (≥99%). The two-stage calcination was performed for 2 h at temperatures 850 °C and 1000 °C. Sintering was carried out for 3 h at 1140–1200 °C depending on composition.
The crystal structures of the crushed ceramic samples were determined using a PANalytical X’Pert PRO X-ray diffractometer (Malvern Panalytical Ltd., Malvern, UK).
Weak-field dielectric permittivity measurements were performed in a temperature range from room temperature to above 400 °C in the frequency range from 130 Hz to 1 MHz, using an impedance analyzer HP precision LCR meter 4284A. Samples of the compositions x = 0.050 and x = 0.075 were poled at room temperature before measurements.
Direct measurements of ECE (
Figure S1, Supplementary file) were taken in samples of two different thicknesses for each (1−x)(0.8NBT-0.2BT)−xCT composition: up to 70 kV/cm for samples 0.3 mm thick and up to 100 kV/cm for samples 0.2 mm thick. Silver electrodes were fired on the sample surfaces at temperatures 500–600 °C. Then, 15 s long electric field pulses of various amplitudes up to 2 kV were applied via high-voltage amplifier TREK 609E-6. Temperature change of the measured sample and polarization current were simultaneously measured by Keithley nanovoltmeter 2182A and Keithley picoampermeter 6485, accordingly. In order to perform the measurements of temperature change, a thin copper–constantan thermocouple was glued with silver paste onto the sample surface, which was electrically grounded. For application of voltage to the sample, thin copper wires (0.08 mm diameter) were also glued onto the sample electrodes. The edges of each sample were covered with a breakdown-resistant acryl coating. This coating, together with the glued thermocouple and the wires used to apply an electric field, created passive thermal capacity. The impact of this thermal capacity was evaluated in two ways: (1) by calculating it from weight of the glue, wires, and lacquer, and (2) by comparing the measured temperature change Δ
T values of the samples of different thicknesses (from 0.5 mm, with only minor influence of passive thermal capacity to 0.2 mm) with the same composition. It was found that the measured Δ
T value of 0.3 mm thick samples is 23% lower, and 41% lower in for 0.2 mm samples—compared with the real temperature change caused by ECE. The shape of time dependence of the temperature change Δ
T(t) was carefully followed for absence of narrow pulses induced by change of electric field. The measured Δ
T(t) curves reflected slow thermal decay of temperature difference (
Figure S2, Supplementary File), related to incomplete thermal isolation between the sample and the environment, without any additional features. During the measurements, electric field dependence of Δ
T was measured at sequential electric field pulses at chosen temperatures. The temperature-dependent measurements were taken in a Delta 9023 chamber up to 150 °C.
Polarization hysteresis loops were measured by the Sawyer–Tower method at triangular pulses, with voltage U(t) slope of 100 V/s. Temperature dependence of remnant polarization was determined by measuring static pyroelectric effect on heating for the samples previously poled at room temperature at E = 70 kV/cm.
3. Results and Discussion
The X-ray diffraction patterns of the studied (1−x)(0.8NBT-0.2BT)−xCT compositions confirm pure perovskite structure. Splitting of pseudocubic [200] maximum for the composition x = 0.050 at 2θ = 46° is seen, while at higher CaTiO
3 concentrations, no deviations from the cubic pattern were detected (
Figure 1).
Temperature dependences of the dielectric permittivity (
ε) of dielectric losses (
tgδ) for all studied compositions are presented in
Supplementary file, Figure S3, and reveal characteristics for NBT-based composition features: diffused frequency-independent maximum (for studied compositions located around 250 °C) and frequency-dependent shoulder at lower temperatures. Well-expressed, frequency-independent maximum in
tgδ dependence on temperature of composition x = 0.05 slightly above 100 °C is related to depolarization temperature (
Td). For composition x = 0.075, corresponding anomaly is transformed into small jump close to 60 °C.
Polarization hysteresis loops for all studied compositions at bipolar electric field pulses at various temperatures are shown in
Figure 2. Well-expressed double hysteresis loops for the compositions x = 0.050 and x = 0.075 are observed in the region of
Td. Signs of a double hysteresis loop at room temperature region are also observed for the composition x = 0.100, while
P(E) for the composition x = 0.125 is almost linear.
Temperature dependences of remnant polarization
Prem(T) for the compositions which are in the ferroelectric phase at room temperature (x = 0.050, x = 0.075) are presented in
Figure 3. The dependences were determined in two ways: (1) calculated from measurements of static pyroelectric effect, and (2) extracted from bipolar hysteresis loops at different temperatures. In both cases,
Prem(T) decreases if temperature is increased, but, in the case of static pyroelectric effect, this dependence is more expressed. The abrupt drop of
Prem(T) around 100 °C for the composition with x = 0.050 and around 40 °C for the composition with x = 0.075 indicates the region of depolarization temperature
Td.
Double hysteresis loops for (1−x)(0.8NBT-0.2BT)−xCT compositions x = 0.050 and x = 0.075 allow us to determine the temperature dependence of the critical electric fields in these compositions (insets in
Figure 3). For both compositions, temperature dependence of the critical electric field
Ec2, corresponding to the phase transition from ferroelectric state to nonpolar state upon reducing the electric field, crosses the
x-axis (
Ec2 = 0). This temperature indicates the temperature at which the poled state remains after electric field is switched off in the direction of low temperatures. Therefore, this is another way to determine the depolarization temperature
Td. For the composition with x = 0.050, T
d is 100 °C, while for the composition with x = 0.075, it is 50 °C. Temperature dependence of another critical electric field,
Ec1—characterizing the phase transition from nonpolar to ferroelectric state upon increasing of field—is less expressed.
Temperature dependences of Δ
T of all studied compositions at 100 kV/cm pulses are presented in
Figure 4. The presented experimental data are measured upon switching off the electric field pulse in order to avoid contribution of the Joule heat in the measured Δ
T, which appears at the highest electric fields and temperatures. The composition with x = 0.050 at room temperature is in ferroelectric state, which is the reason for Δ
T being low. Significant increasing of Δ
T(T) is observed around 100 °C, corresponding to the depolarization temperature, where a drop in remnant polarization occurs (
Figure 2). The diffused jump of Δ
T(T) in the temperature range between 70 and 130 °C is apparently related to a mixed ferroelectric–nonpolar state at
E = 0 below 130 °C. With increasing temperature, concentration of the ferroelectric phase at
E = 0 decreases. In the case of the composition with x = 0.075, the jump is shifted towards lower temperatures and is more pronounced. The broad maximum of Δ
T(T) indicates that the applied electric fields are sufficient to induce ferroelectric state in a wide range of temperatures above
Td. In the compositions with higher concentrations of CaTiO
3,
Td is no longer found in the temperature range above room temperature, and Δ
T only weakly depends on temperature, while Δ
T values are significantly reduced.
Electric field dependences of Δ
T for all studied compositions at different temperatures are shown in
Figure 5. For the composition with x = 0.050 (
Figure 5a), at measurement temperatures up to 90 °C, where the composition is in the ferroelectric state, Δ
T is low across the whole electric field
E range, with weakly negative curvature of Δ
T(E). Above this temperature range, Δ
T(E) contains significant contribution from the jump, which is observed in the vicinity of the critical electric field
Ec and corresponds to the entropy jump
(∆
SPT) characteristic for the first-order phase transition:
where
cEc is heat capacity at the critical electric field. The Δ
TPT value, roughly evaluated from Δ
T(E), is 0.65 °C. On the other hand, Δ
TPT can be evaluated indirectly from the Clausius–Clapeyron equation:
Considering that Equation (2) is derived from the condition that chemical potentials of both phases are equal, it does not take into account the actually-observed hysteresis of the critical electric field values upon increasing and decreasing the electric field (
Ec1 or
Ec2—accordingly) due to metastability of both states. For this reason, temperature dependence of the critical electric field
Ec(T), which could be used in Equation (1), should have values between
Ec1(T) and
Ec2(T). In order to compare directly and indirectly determined values of Δ
TPT for the composition with x = 0.050, we have used the polarization jump at phase transition Δ
PPT and slopes of
Ec1(T) and
Ec1(T), as well as heat capacity of 0.92NBT–0.08BT at
E = 0 [
39], which should be close to the heat capacity of the measured compositions. Near the depolarization temperature,
dEc1/dT = 475 V/cm and
dEc2/dT = 775 V/cm. Δ
TPT calculated from these values is 1.3 and 2.2 °C accordingly—which is significantly higher than the directly measured value. The reason for such a difference could be the nature of the relaxor state existing in this case [
38], which can be considered as a mixture of polar (polar nanoregions) and nonpolar phases instead of a pure nonpolar (paraelectric) phase.
For the composition with x = 0.075, transformation of Δ
T(E) with increasing or decreasing temperature (
Figure 5b) follows the same trend as for the composition with x = 0.050. Obviously, due to lower
Td, related to the phase transition induced by the electric field, the jump of Δ
T(E) is shifted towards lower temperatures, while the value of the jump is similar.
At higher temperatures, ΔT(E) for the compositions with x = 0.050 and x = 0.075 becomes diffuse, losing signs of electric field-induced first-order phase transition, and its values at maximal electric field pulses decrease significantly. At the same time, moderate increasing of ΔT due to the Joule heat, upon application of an electric field pulse, is observed.
Δ
T(E) for the composition with x = 0.075 (
Figure 5c) maintains the inflection point, presumably reflecting the diffuse transfer to the partly ferroelectric state, which almost disappears with increasing temperature. In the composition with x = 0.125, the inflection point disappears completely at higher temperatures (
Figure 5d).
The difference between Δ
T values with the electric field switched on and off is not related to Joule heat alone. This behavior is illustrated in the case of the composition with x = 0.075 at
T = 40 °C (
Figure 6a), where a significant difference between Δ
T values appears already at moderate temperatures and electric fields. The reason for such a difference is slow relaxation of polarization observed in the unipolar
P(E) cycle (
Figure 6b), which is responsible for
P(E) hysteresis and persists even after the cycle is finished and electric field reaches
E = 0. Relaxation is complete within a few minutes and can be followed by measurement of polarization. However, the corresponding contribution to Δ
T is outside the peak of Δ
T(t), which is used to measure Δ
T and is determined by thermal relaxation between the sample and thermocouple. This slow relaxation of polarization is very expressed in the temperature range where the electric field induced phase transition occurs.
Experimentally measured polarization data are frequently used to determine Δ
T by the Maxwell relation. In the ferroelectric state, data from the upper branch of bipolar hysteresis loops in the electric field ranging from the maximal value
Emax to 0 are used for this purpose. Earlier, we noted that polarization, especially in the case of bipolar electric field pulses, depends significantly on ferroelectric domains switched by an electric field. Obviously, a switched domain pattern at different temperatures can be different, while polarization at fixed values of
E (extracted from polarization hysteresis loops) upon change of temperature does not describe the continuous development of the same domain state. Especially in the case where
Emax does not significantly exceed the coercive electric field, temperature dependence of concentration of switched domains is remarkably reduced comparing to the fully oriented case. As the coercive field usually decreases upon approaching the phase transition, concentration of the domains oriented in direction of electric field also increases. It can lead even to increased polarization upon increasing of temperature, which results in negative ECE, sometimes determined in NBT-based compositions by the indirect method [
28,
29,
30]. At first sight, the problem with dependence of remnant polarization on domain switching can be avoided using measurements of static pyroelectric coefficient. However, unfortunately also in this case, contribution of domain reorientation or, more presumably, irreversible shift of domain boundaries contributes to polarization. These problems are illustrated in
Figure 3. Temperature dependence of remnant polarization extracted from polarization hysteresis loops is more weakly expressed, compared with remnant polarization determined from measurements of static pyroelectric effect, and the corresponding
Td is shifted towards higher temperatures. At the same time, temporary change of direction of temperature for the composition with x = 0.050 around 50 °C (inset in
Figure 2) clearly reveals irreversibility, related to irreversible domain processes. On the other hand, incomplete switching of domains, which should be quite a common behavior in polydomain ferroelectrics, reduces the value of Δ
T compared with the fully aligned domain pattern. Therefore, both directly measured values—∆
T and
—can be far from thermodynamic conditions where direct and indirect values can be compared. Their coincidence in the ferroelectric state within the electric field range usually accessible for bulk ferroelectrics (where contribution of domain processes is expected) could be rather accidental.
Outside of a stable ferroelectric state, polarization consists of solely electric field-induced polarization, and evaluation of Δ
T from the Maxwell relation seems more appropriate. As unipolar pulses are used for measurements of Δ
T,
P(E) dependence could be also more adequate if measured at unipolar pulses. In the present case, the Maxwell relation will be applied to temperature dependence of polarization at different electric field values, obtained from measurements of polarization current at the same electric field pulses which were used in the measurements of Δ
T (
Figure 7).
From our point of view, comparison of values of Δ
T obtained by the direct and the indirect methods is not the best way, because integration of
does not allow one to extract correspondence between both methods at different ranges of electric field. For that reason, differential layout of the Maxwell formula (Equation (1)) was transformed, allowing us to compare
dT/dE and
.
As with the calculation of the jump of Δ
TPT, dependence of
cE on the electric field is neglected and values of heat capacity are borrowed from [
39].
Examples of comparison of the results obtained by the direct and the indirect methods are presented in
Figure 8. In the case where Δ
T(E) has a clear enough inflection point, the shapes of curves representing the left and the right side of Equation (3) mutually correlate. However, the differences between these values can be significant. If the inflection point is not as pronounced (maximum of
dT/dE as a function of E is low), even the shapes of the two curves do not correlate. This might partly be explained by rather flat
P(T)E dependence in this case (low values of
accordingly) when any additional contribution to polarization can seriously influence the values of
. From another point of view, as it was assumed above, the expressed inflection point could reflect at least a partial phase transition between nonpolar and ferroelectric phases, which could present a more appropriate case for using the Maxwell equation rather than rearrangement of polarization in the nonpolar phase clearly representing the relaxor state.